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Perturbative Uberpriifung der Wirkungs- und 
Energie- Gittersummenregeln 

Z usammenfassung 

Gittersummenregeln werden mittels Storungstheorie auf dem Gitter iiberpriift. 

Die Wirkungssummenregel liefert eine Beziehung zwischen dem Quark-Anti- 
quark-Potential, dessen logarithmischer Ableitung nach dem Abstand und dem 
Erwartungswert der Wirkung; die Energiesummenregel driickt das Potential als 
Summe der Energie in den Gluonenfeldern und eines anomalen Terms aus. Zwei 
verschiedene unabhangigc Bcrechnungen des Quark- Antiquark-Potentials wer- 
den vorgestellt, und die Transversalitat der gluonischen Vakuumpolarisation 
auf dem Gitter wird bewiesen. Der wesentliche Teil der Wirkungssummen- 
regel ist eine Identitat, deren explizite Uberpriifung Methoden und Ergebnisse 
zur Verfiigung stcllt, die bci der Bchandhmg der Energiesummenregel niitzlich 
werden. Zusatzlich wird die Eichinvarianz des Erwartungswertes des Wilson- 
Loops bis zur nachstfiihrenden Ordnung bewiesen. Die Mogliclikeit, den Er- 
wartungswert der Wirkung auf den Erwartungswert der Summe der Plaketten 
zu cincr festen Zeit einzuschranken, wird diskuticrt. Die Energiesummenregel 
wird storungsthcoretisch bis zur nachstfiihrenden Ordnung iiberpriift, und cs 
wird gezeigt, dass sie mit guter numerischer Genauigkeit erfiillt ist. Die einzel- 
nen Beitrage zum Quark-Antiquark- Potential werden analysiert, und die Ein- 
schrankung des Erwartungswertes der Summe iiber alle raumlichen Plaketten 
(die Energie in den magnetischen Fcldcrn) auf den Erwartungswert der Summe 
der raumlichen Plaketten zu einer festen Zeit wird untersucht. 

Perturbative Check of the Action and Energy Lattice Sum Rules 

Abstract 

Lattice sum rules are checked using lattice perturbation theory. The action 
sum rule gives a relation between the quark-antiquark potential, its logarithmic 
derivative with respect to distance and the expectation value of the action; 
the energy sum rule expresses the potential as the sum of the energy in the 
gluon fields and of an anomalous term. Two different independent calculations 
of the quark-antiquark potential are presented, and the transversality of the 
gluonic vacuum polarization on the lattice is proven. The crucial part of the 
action sum rule is an identity whose explicit check using perturbation theory 
provides methods and results which are useful for checking the energy sum 
rule. Additionally, the gauge invariance of the expectation value of the Wilson 
loop up to next-to-leading order is proven. The possibility of restricting the 
expectation value of the action to one fixed time slice is discussed. The energy 
sum rule is checked perturbatively up to next-to-leading order and shown to 
be satisfied with good numerical accuracy. The various contributions to the 
quark-antiquark potential are analyzed, and the restriction of the expectation 
value of the sum over all spatial plaquettes (the energy in the magnetic fields) 
to one fixed time slice is examined. 



There are very beautiful and elegant ways of getting these things 
these days; but suppose that you were inventing it, what would you 
do to find [it]? You fiddle around. All the elegant stuff is found 
later; the way to learn is not to learn elegant things, it's to fiddle 
around blind and stupid. Later you see how it works; polish it up; 
remove the scaffolding and publish the result for other students to 
be amazed at your ingenuity. 



Richard Feynman 
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Chapter 1 

Introduction 



1.1 Quantum Chromodynamics 

Only 40 years ago, protons and neutrons were still assumed to be elementary 
particles. But already in the fifties of the 20th century, more and more new 
so-called "elementary" particles had been found at high energy colliders (see, 
for example, eventually dozens were known, most of them with very short 
life times. 

The turning point came 1964 with the invention of the quark model [H 
in], which postulates that all of these particles (nowadays known as hadrons) 
are composed of sub-particles, the quarks. With this model, the variety of 
hadrons could be sorted in a systematic way: the particles with integer spin, 
the mesons, consist of a quark-antiquark-pair, the particles with half-integer 
spin, the baryons, are bound states of three quarks or three antiquarks. Today 
hundreds of hadrons are known, and most of them fit well into the quark model. 
The few exceptions can be attributed to so-called gluehalls and other exotic 
states. 

Early experimental support came from high-energy scattering experiments, 
which showed that there are sub-particles, called partons, in the proton [3]. The 
scale invariance which was observed in these experiments was well explained by 
the quark model 0- Based on these and many other experimental results, 
nowadays the quark model is universally accepted; a good introduction can be 
found in jH], for instance. 

The theoretical description of the quarks and their interactions is provided 
by Quantum Chromodynamics (QCD) [TIIH]. In this quantum field theory, the 
quarks are described by Dirac spinors with an additional degree of freedom 
called color; they are coupled to SU{3) gauge fields. The quantization of the 
bosonic gauge fields leads to the eight spin-one particles called gluons which 
mediate the interactions of the quarks. A prediction of QCD is the existence 
of bound states from these gluon fields, the already above mentioned glueballs. 
An early theoretical review of QCD can be found in 0; for a review of the 
experimental support see, for instance, ^0] and PT] . 

A large problem remains, however: no free quarks or gluons ever have been 
observed, despite all efforts taken in the years following the invention of the 
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quark model ^Ij. The standard explanation for this is that quarks and gluons 
can only exist in bound states; this phenomenon is called confinement. 

It can be explained if one assumes that the force between the quark-anti- 
quark pair does not decrease with increasing distance as in electrodynamics, 
but instead stays constant — a flux tube or string of gluons forms between the 
quarks. Therefore, if the distance is increased, a point will be reached at which 
the gluon string contains enough energy for the generation of a new quark- 
antiquark pair out of the vacuum, breaking the string. Hence instead of sepa- 
rating the original quark-antiquark-pair, one has produced two new pairs (this 
effect is called screening of the color charge). Recent Monte Carlo simulations 
on the lattice (see below) confirm this picture |13j . 

1.2 Lattice Gauge Theory 

The most commonly used method in Quantum Field Theory is perturbation 
theory. But for many important problems, unfortunately including quark con- 
finement, perturbative expansions fail: In QCD, the coupling constant is small 
at high energies, corresponding to small separations of the quarks (this is known 
as asymptotic freedom and is based on the non-abelian nature of the gauge 
group), but large at the low energies and relatively large distances occuring in 
hadronic bound states. Therefore many non-perturbative methods have been 
developed in the past decades in order to treat confinement. 

Most promising is the formulation of QCD on a space-time lattice: space 
and time are discretized, so that the number of degrees of freedom is reduced 
drastically. Usually one works in euclidean space-time. If additionally one only 
looks at a finite volume and a finite time interval, there is only a finite number 
of degrees of freedom left, and the generating functional for this discretized 
version of QCD can in principle be calculated on a computer. In the limit of 
infinitely small lattice spacing, the results one obtains should reduce to the 
continuum results. 

Lattice formulations of QCD where first suggested 1974 by Wilson jTl], 
Kogut and Susskind ^^l- In the three decades since these proposals, lattice 
gauge theory has become a branch of particle physics in its own right. Addi- 
tionally, it is closely connected to statistical mechanics and therefore of interest 
not only to particle physicist, but also to physicists working on solid state 
physics and other many-particle problems. Today, there are lots of encouraging 
results from lattice calculations, including the quark-antiquark potential, the 
string tension, hadron masses, the QCD phase transition from the hadronic 
bound states to a quark-gluon plasma and so on. An introduction to all of 
these concepts and a summary of the results can be found in JHI) for example. 

But one still has to use a very large number of degrees of freedom in order to 
get a reasonable approximation to continuum, infinite volume QCD. Because of 
this, very high-dimensional integrals have to be computed. For example, if one 
uses a lattice with an extension of only 10 in each of the four directions and for 
simplicity considers only the gauge fields, there are already 320.000 degrees of 
freedom (10^- 4 (polarization) • 8 (color)). Hence a 320.000-dimensional integral 
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has to be performed. Using only 10 points in every direction in this 320.000- 
dimensional space in order to evaluate this integral numerically, one would have 
to compute and add 10'^^'^'''^'^ terms! 

Obviously this is not feasible; hence the numerical integrals are usually done 
by Monte Carlo methods. Most gauge field configurations have a large action 
and thus only contribute very little to the generating functional. Therefore the 
functional integral can be approximated by generating a sequence of gauge field 
configurations with a probability distribution given by the Boltzmann factor; 
this technique is called importance sampling. Then the expectation values of 
operators can be approximated by their ensemble average over these represen- 
tative configurations. 

Simulating the gauge fields is relatively simple; the quarks pose bigger 
problems — for example, fermion doubling [17| and the sign of the fermionic 
determinant. Therefore one often works in the so-called quenched approxima- 
tion: one neglects the effects of dynamical quarks. Obviously the screening 
and the string breaking due to the creation of quark-antiquark pairs can not 
be observed in such simulations, but they still can give valuable results on the 
short distance behaviour of bound states. 

For a description of the commonly used Monte Carlo algorithms and a review 
of the most important results, see e. g. |16j . 



1.3 Sum Rules for the Quark- Antiquark-Potential 

Although the lattice results reproduce the expected behaviour of the quark- 
antiquark potential, including the confinement, quite nicely and give good esti- 
mates for the string tension, an analytic description would be preferable instead 
of these solely numerical results. In order to gain some insight into the quark- 
antiquark potential and the corresponding energy densities of the gluon fields, 
in 1987 C. Michael invented lattice sum rules These connect the quark- 

antiquark potential with the expectation values of the energy in the gluon fields; 
even a separation of the electric and magnetic contributions is possible. 

However, Michael's derivation was partly incorrect (see below); this was 
pointed out by Dosch, Nachtmann and Rueter in 1995 Shortly afterwards 
H. J. Rothe published a corrected version [53], and when Michael extended 
his sum rules in 1996 (2^, he took this correction into account. Additionally, 
Rothe discovered that the trace anomaly of the energy momentum tensor gives 
an important contribution to the quark-antiquark potential |24j . 

The main goal of this dissertation is a perturbative check of these sum rules. 
In the following subsections, their derivations will be presented, based on |23j . 
[24| and then I will outline the problems with the sum rules which will be 
investigated using perturbation theory and give an overview of my approach to 
solve these problems. Some results will be presented for a general gauge group 
G, but in most cases I will restrict myself to the commonly used group SU{N). 
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1.3.1 The Wilson Loop and Wilson's Action 

The starting point for the derivation of both sum rules is the standard formula 
for the quark-antiquark potential V: 

V{R) = - lim i In < W{R, f) > . (1.1) 

f^oo T 

The Wilson loop W{R, T) is the path ordered product of the link variables Ui 
along a closed rectangle C with spatial and temporal size R and T, respectively: 

The link variables Ui are elements of the gauge group, and d{R) is the dimension 
of the representation R of the group in which the sources of the gauge field are. 
Usually the fundamental representation F is used; then d{R) = d[F) = N for 
the group SU{N). 

Wilson introduced this loop in 1974 as an order parameter for Quantum 
Chromodynamics ^1]; a similar structure was already used in 1971 by F. J. 
Wegner for the Ising model jl8j . The derivation of the formula (jLlj) which 
connects the quark-antiquark-potential to the expectation value of the Wilson 
loop can be found in JH] and in the review article [20], for example. 

The potential and the size parameters of the Wilson loop are measured in 
lattice spacings a: 



R 

a 
T 



R 

f 

a 

V = V-a. (1.3) 



The expectation value in (|l.lj) is calculated with Wilson's action (the sum over 
all plaquettes) in the quenched approximation (no dynamical fermions): 

^ = E (l - ^Tr([/^4n) + C/^,(n))) =: /3P (1.4) 
with the lattice coupling constant 

^ = ^ (1.5) 

and the plaquette variables 

U^.H = U^{n)U,{n + fi)Ul{n + ml{n). (1-6) 

where jl resp. v are unit vectors pointing in the corresponding direction. The 
sum in the action runs over all lattice sites n and all possible orientations of 
the plaquettes. 
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1.3.2 Action Sum Rule 

Taking the derivative of with respect to the logarithm of the couphng 

constant, one arrives at: 

r = /?— r = hm ^ < 5 >gg-0, (1.7) 

where for a general operator O, the expectation value in the quark-antiquark 
state with respect to the vacuum state is defined as 

< OW > 

< o >,,_o= <o>- (1.8) 

Now the scaling behaviour ()1.3|) of the potential can be used. In the continuum 
limit a ^ 0, the physical potential V obviously should be independent of the 
lattice spacing. This gives: 



for a ^ 0. (1.9) 



Carrying out the differentation, using the scaling behaviour of R, the definition 
H1.5|) and introducing the lattice beta function 

PUgo) = -a^, (1.10) 

one arrives at 

^J^fZ = V + R?i. (1.11) 

go dp dR 

Using this in (|1.7|) . the result is: 




y + i?^ = ^ hm i<5>,,_0 

dR go f-^ooT 



(1.12) 



This is essentially the action sum rule for the quark-antiquark potential which 
was derived by H. J. Rothe in 1995 |23j . 

As already pointed out, Michael derived this sum rule first (in 1987, HJ), 
but instead of (|l.llj) . he used a wrong scaling behaviour for the potential, and 
therefore in his final formula, the logarithmic derivative of the potential with 
respect to R was missing. For a confining potential, V = aR, this is crucial: 
using the correct sum rule, the left hand side gives 2V, whereas Michael's wrong 
sum rule only had V and therefore was wrong by a factor of 2. 

On the other hand, in his 1987 paper, Michael claimed that the expectation 
value < S >qq-o on the right hand side can be further simplified in the limit 
T ^ oo (for a more detailed derivation, using transfer matrix methods, see 
[29j): plaquettes on different time slices should give the same contribution to 
the expectation value, so that 

< S >gq-0^ f < L{t) >gg_0, (1.13) 
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where L(t) is the sum over all plaquettes at a fixed time slice t: 

Obviously this is the Lagrangian at time t. 
Therefore the action sum rule becomes: 

V + R?Xr = ^ lim < L{t) >,,-_o . (1.15) 
dR 90 f^oo 

Rothe noticed 1995 [23 that this can be rewritten using the trace anomaly 
of the energy momentum tensor. This was motivated by Ji's observation that 
this trace anomaly contributes 1/4 of the hadron masses j26j . Already in 1977 
it had been shown by Collins, Dunkan and Joglekar |27j that the trace of the 
energy momentum tensor in the continuum formulation of SU{N) gauge theory 
is given by 

^ 2g g ^ ^ 

in the quenched approximation, where C is the Lagrangian density, which is 
exactly the combination appearing in the formula above. In 1992, Caracciolo, 
Menotti and Pelisetto showed that this is also true in lattice perturbation theory 
|28j . Using this one finally gets for the action sum rule: 

- dV 

V + R-^= hm ^ < T,,^{x, t) >,,-_o . (1.17) 
1.3.3 Energy Sum Rule 

In order to derive the energy sum rule, one has to use an anisotropic lattice 
with lattice spacing at in the temporal direction. The anisotropy parameter is 
defined as 

i ■■= at/a. 
Then the action has to be rewritten: 

=: psVs + PtVt, (1.18) 

where one has to use different coupling constants for the spatial and the tem- 
poral plaquettes. 

In the continuum limit, the expectation value of the Wilson loop should 
be independent of the asymmetry, as long as the extension in the temporal 
direction is the same: 

< W{R, f) >5=i=< W{R, if ) >^ for a ^ 0, 
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and accordingly 
with 



V{R,(3s{i),m)) = - lim -.\n<W{R,if)> . 

5T^oo ST 



It follows that in the continuum limit 

d 



iv{R,m).m)) 







(1.19) 



should be satisfied. Carrying out the differentation and then returning to the 
isotropic lattice by setting = 1, one arrives at: 



lim — 

f^oo T 



V{R,(3)-- 
With the abbreviations 



-Q^ < Pt >qq-0 +~Q^ < >qq-Q 



(1.20) 



J 5=1 



'-=2 



(1.21) 



this can be rewritten as 



V{R, f3) = lim ^ [7?_ < -Vt + Vs >qq-0 <Vt+Vs >qq-0] ■ 



In j42j . Karsch had shown that 



4 go 



(1.22) 



by requiring that the string tension calculated from space-time-like Wilson loops 
should be identical to the one calculated from purely spatial ones. Hence the 
potential becomes: 



V{R, p) = lim i 

f^oo T 



V- < -T^t + Vs >qq-0 +- /3 <Vt + Vs >qq-0 



4 go 

where obviously the second term contains the expectation value of the action: 



V{RJ) = lim i 



V- < -Vt + Vs >qg-0 + 



1 2/3L(ffo) 
4 go 



<S> 



qq-O 



[1.23) 



This is essentially the energy sum rule which Rothe derived in |23 . Using the 
relation between the action and the trace of the energy momentum tensor, it 
can be rewritten as 



V{R, /3) = lim 



1 1 

V^J< -Vt + Vs >qq-0 +4 < t) >qq_o 
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On the other hand, Karsch [17 showed that in the continuum Hmit (go 0): 

and if one again assumes that the expectation value for plaquettes on different 
time shces is identical for large T, one can rewrite the first expectation value: 

lim 1<-Vt + Vs >gg-o= lim < -V^it) + V',{t) >q,-_o, 

T->oo T T->oo 

where V'g f{t) is the sum over all spatial respectively temporal plaquettes on the 
fixed time slice t. Therefore in the continuum limit, the first term on the right 
hand side of ()1.23|) reduces to: 

lim r/„i < -Vt + Vs >qg-o^ 1™ < -^^i^) + ^^(^) >gg-o, 

T^oo T i ^ '^^'^ 

(1.24) 

which is just the euclidean version of the energy in the gluon fields. Hence the 
energy sum rule tells us that the potential energy of the quark-antiquark pair 
is given by: 

energy in the gluon fields + -trace anomaly of the energy momentum tensor. 

An especially interesting case to look at is a purely confining potential: 

V{R) = ak. 
From the action sum rule, one can deduce that 

- dV 

hm ^ < r^^(f , t) >gg_o= V + R-^ = 2V; 
hence the energy sum rule gives: 

IV = - lim r?_i < -Vt + Vs >qq-0 ■ 
^ T^oo T 

Therefore for a confining potential, the contributions from the energy in the 
gluon fields and from the trace anomaly of the energy-momentum-tensor have 
the same magnitude: both account for exactly one half of the total potential 
energy. 

1.4 Problems with the Sum Rules 

The main goal of this work is to check the lattice sum rules. This is necessary 
because of several reasons: first, neither for the argument leading to the restric- 
tion to one fixed time slice, nor for the scaling arguments on the anisotropic 
lattice, it is known if they are justified. Hence their validity should be checked. 

But beyond these validity doubts, there are also other reasons to take a 
closer look at the sum rules. Because of the opposite signs of the magnetic 
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and electric field energies in the euclidean formulation, cancellations between 
these two contributions can occur when their expectation values are calculated 
in Monte Carlo simulations. These cancellations are a possible source of errors; 
a more thorough investigation of what happens there exactly should be helpful. 
Additionally, it would be interesting to see how do the contributions of the trace 
anomaly of the energy-momentum tensor look like. 

There are in principle two ways to check the sum rules and to examine 
these interesting details. On the one hand, one could perform Monte Carlo 
simulations. Using these, one could investigate the physically interesting region 
where the coupling constant is large. But there are some drawbacks, too: much 
computer power and time is needed for these simulations, and the results are 
not very enlightening in general; numerical results do not tell as much about 
the underlying physics than analytical calculations. 

Because of these reasons, in this work lattice perturbation theory is used to 
study the sum rules. Obviously this is only possible for weak coupling, hence 
the physically interesting and relevant region is missed. But this is countered 
by the fact that much less computer power and time is needed (only some 
low-dimensional numerical integrals in momentum space have to be done). Ad- 
ditionally, one gets analytical results, which are expected to shed some light on 
the problems mentioned above. 

The Feynman rules for lattice perturbation theory are well-known and can 
be taken, for example, from JHI; there the necessary weak coupling expansion 
of the action is also explained in detail. The consistency of these perturbative 
calculations, a general power counting theorem, and considerations about the 
continuum limit can be found in 38 . The last necessary ingredient for the 
check of the sum rules in the weak coupling approximation is an expression for 
the potential up to next-to-leading order, which is provided in '3^1 . 

The outline of the check is as follows: first, I will explain how one obtains a 
weak coupling expansion of the quark-antiquark potential in lattice perturba- 
tion theory; I will outline an approach by Kovacs to this problem |3U| as well as 
the methods used by Heller and Karsch, already mentioned above 1391. In the 
same chapter, I will prove additionally the transversality of the gluonic vacuum 
polarization, which appears in the calculation of the potential. 

Then I will turn to the check of the action sum rule. In principle the validity 
of this sum rule is clear - after all, it is more or less an identity. Nevertheless, 
it seemed to be an easier task to check it rather than the energy sum rule, 
and it provides an opportunity to demonstrate some of the necessary methods. 
Additionally, it will turn out that this sum rule is closely connected with the 
gauge invariance of the Wilson loop — only some small additional calculations 
will be needed to check this gauge invariance, too. And there is one part of the 
action sum rule for which a check is really needed: the validity of the restriction 
to one time slice, which was already mentioned above, is not clear at all even 
for the action sum rule. 

Finally the energy sum rule also has to be checked. Many results of the check 
of the action sum rule will become useful there (a further reason to check that 
sum rule first), hence this check will be easier than the check of the action sum 
rule, although it looks more complicated at first sight. After the check, I will 



10 



CHAPTER 1. INTRODUCTION 



explain how one can derive results about the sizes of the various contributions to 
the potential (electric, magnetic, anomal) from my calculations and the energy 
sum rule. I will close that chapter again with a section on the restriction to one 
time slice. 

The last chapter contains a summary and the conclusions. In the appendix, 
useful formulas are summarized, and some basic calculations are presented. A 
detailed account of the Fourier transformation needed in chapter 2 for calculat- 
ing the potential can also be found there. 



Chapter 2 

The Quark- Antiquark 
Potential 



In this chapter, I will line out how one can calculate the static quark-antiquark 
potential on the lattice, which is needed in order to check the sum rules. First I 
will give a general derivation, based on the running of the coupling constant and 
renormalization group arguments (for more details, see ISHI)) then I will discuss 
how one can derive these results directly by a perturbative lattice calculation. 
Additionally, I will prove that the gluonic vacuum polarization is transversal in 
leading order of the lattice perturbation theory. 



2.1 The running coupling constant 



In leading order, the quark-antiquark potential is given by: 

^2 



V{R) 



9 



AttR 



C2{F), 



(2.1) 



where g is the coupling constant of the strong interaction and C2{F) the 
quadratic Casimir operator of fundamental representation F; for SU{N): 



C2{F) 



1 



2N 



(2.2) 



If one takes into account that the coupling constant depends also on the 
distance between the quark and the anti-quark ("running coupling constant"), 
then the formula above needs corrections: 



V{R) 



9HR) 
4ttR 



C2{F) 
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AirR 



C2{F) + 0{g' 



(2.3) 



There are various methods to calculate these corrections. In .'{5 . they were 
determined in momentum space, using dimensional regularization and the MS 
renormalization scheme (modified minimial subtraction: not only the pole in e 
is subtracted, but also terms of ln47r). The result is: 



V{t) 



-C2{F) 



7 + 



31 

33 



(2.4) 
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Here the Feynman gauge (gauge parameter a = 1) and the quenched approx- 
imation (number of dynamic fermions Uf = 0) was used. is the couphng 
constant defined in the MS scheme, fi a mass parameter which has to be intro- 
duced during renormahzation, 7 ~ 0.577216 is the Euler-Mascheroni constant 
and /?o the first coefficient in the expansion of the beta function: 



/3o 



11C2(G') 
487r2 - 



(2.5) 



where C2{G) = N is the quadratic Casimir operator of SU{N) in the adjoint 
representation. 

From this one can get an expression for the potential which should result 
from lattice calculations by using the relation between g^jg and the coupling 
constant go{a) appearing in perturbative calculations on the lattice. This is 
done in two steps; first the following relation can be used: 



1 



A{a,nf,N) 2 ( 2^ , r>( 4^' 
9MOMij^ )+0{g ) 



Air 



where gMS and gMOM are the coupling constants defined in the MS and 
MOM (momentum-space subtraction) renormahzation scheme, respectively. 
A(a,nf, N) was determined in [5T| for various values of the gauge parameter a 
and numbers of dynamical fermions nf, the value A{1,0,N) which is relevant 
here is given by: 



A{1,0,N) 



C2{G) 
2ti 



11 11 ^ 23 1 

7 H ln47r H \ / 

6 ' 6 6 36 



with 



Inx 



■c"^ — X + 1 



2.3439072. 



But here gj^g is needed; for this one obtains: 
with 



1 5mom(/^ ) + 0(5 



, , 11C2(G) C2(G) 
^(1,0,A^) = A(1,0,A^) r4^1n47r- ^ 



11 23 1 

7 + — + 



6 



6 36 



(2.6) 



(2.7) 



127r 2tt 

The second step is to connect gMOuifJ''^) with go{a) in the continuum limit 
a ^ (31031: 

^2 



a 



TT 



1 + g'oia) /3o In ^ + R{N) + 0(5^ 



(2.8) 



with 



R{N) 



C2{G) 
C2{G) 



1 



487r2 
1 

487r2 



23 - 221n7r-F-/j +2P 

22,- 22\mT + -I] +2P -- 
6 



1 

w 



+ -C2(F), (2.9) 
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where / is defined as above and P ~ 0.0849780. 

Putting everything together, the result for the potential in momentum space 

is: 

ffo(«), 



-C2{F) 



l + 9o(«) 



/3o In 



vr^ 31- 
aV ^ ^ 33 



A{l,^,N) 
4tt 



+ RiN) 



In order to get the potential in coordinate space, one has to do a Fourier trans- 
formation (for details, see appendix ID)) : 



ViR) 



alia) 



C2{F) 



AttR 
l + 5o(a) 



/3o In 



a2 +^+33 



3l\ A(1,0,A^) 



47r 



+ R{N) 



hence the potential, measured in lattice spacings, is: 
5o(a), 



V{R) 



AttR 



■C2{F) 



(2.10) 



l+5o(«) 



2/3o In f7re^/2+3i/66^A _ M^Ml + r(^n) 



Inserting the explicit expressions for A{1,0,N) and R{N) given above, one 
finally arrives at: 



V{R) 



9oia) 
AttR 



C2{F) 



(2.11) 



l + 5o(a) 



or 



go(a) 
A-kR 

^2 



C2(F) 



2/3oln (e^+3i/66+48.^P/ii-3.Vii^) + lc'2(F) 

C72(G)ln (7.501i?) + ^C2(F) 



l+5o(«) 



11 
24^ 



y(i?) 



AttR 

alio) 

AttR 

alia) 

AttR 



J-C2iF) [l + 2gl{a)Po In (7.501e^^(^)/s^oi?) 



C2(F) 1 + 2(7^(0)^0 In (7.501e' 



l + 25g(a)/3oln (56.47i?J 
l + 2c/§(a)^oln (82.07-RJ 



„37r2(Ar2-l)/llAr2 

5f/(2) 
SUi3) 



R 



(2.12) 



Alternatively, one can express these results by using the QCD scale param- 
eters A of the different renormalization schemes; the following relations hold 



^(l,0,iV) 



2/3o In ■ 



A 



MOM 



4tt A 
R{N) = 2/3oln 



MS 
J^MOM 



(2.13) 
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With this, one can write for the potential 
9o{a) 



V{R) 



AttR 

9o{a) 
AttR 



C2{F) 
C2{F) 



l + 25^(a)/3oln 



,7/2+31/66^ 



MS 



R 



,A 



1 + 2g^{a)po In ( 2.135^^ i? 



The ratio of the scale parameters is: 



A 



MS 



Al 



26.45 SUi2) 
38.45 SU{3) 



(2.14) 



(2.15) 



These expressions for the potential agree with results from Monte Carlo 
simulations j3Uj . as well as with a perturbative calculation on the lattice |39j . 
Additionally, the potential shows the right scaling behaviour, compared with 



In some older publications, one can find the following simpler expression for 
the potential instead of the one derived above 



V{R) 



AttR 



C2{F) l + 2gl{aWnR 



(2.16) 



the (big) factor 2.135 J^'^ in the logarithm is missing there completely! There 
are two reasons for this: 

1. The contribution from the tadpole graph in the gluonic vacuum polar- 
ization was left out by these authors; therefore the term proportional to 
C2{Ff is missing. That term contributes a factor 7.529 {SU{2)) or 10.942 
{SU{2>)), respectively, in the logarithm. 

2. For the ultraviolet and infrared cutoffs, they used only rough approxima- 
tions, choosing them idential to a and R, respectively; hence the ratio 
of the cutoffs is identical to R in their results. In contrast, the correct 
calculation yields an additional factor of 7.501 in the logarithm. 



2.2 Weak coupling expansion on the lattice 

A weak coupling expansion of the quark-antiquark potential on the lattice was 
done first by Heller and Karsch in 1985 [SHI- Their methods and results are 
summarized in this section for future reference. 

As already mentioned in the introduction, the quark-antiquark potential 
can be derived from the expectation value of the Wilson loop: 

V{R, /3) = - lim i In < W{R, f) >subtr, (2.17) 

with 

W = -T^Tr Yl Ui. (path ordered) (2.18) 
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The subscript "subtr" means that in order to get the correct potential from the 
formula above, one has to subtract the self energy contributions to the quarks. 
This subscript will be omitted in most places in the following. 

The rectangular curve C is chosen to have the corners no, no + Rfi, uq + 
Rjl + Ti>, no + Tz), where /i and z) are the unit vectors in one spatial and the 
temporal direction of the lattice, respectively. The link variables are connected 
with the vector potential by the following relation: 

^^ = eWA,_ (2.19) 

Using the Baker-Hausdorf formula repeatedly, one then arrives at the fol- 
lowing expression for W: 



W= ^Trexp 



igoJ^Ai - ^ {Ai^,Ai^\ - -^igl ^ [[A,,,^,,], ^,3] 

I h<l2 h<l2<l3 



{h,l2)<l3 h<l2 

where Ai is given by: 

Ai = Afj_{n) for / = (n, n + /i) 

Ai = — A^(n — /x) for I = {n,n — fi). 

Expanding the exponential and taking the trace, one obtains the perturbative 
expansion for the Wilson loop: 



with 

,(2) 



W = l- 5ga;P) - 5o^a;(^) - ^o^a;^ + 0(^0^) (2.20) 



E^n (2-21) 



M{F) 



TMEA +:r7^Tr E^^EIA.A,] (2.22) 



Qd{F) VV / 2(i(F) 



I h<l2 



■ZTTinTrjE^' E [[Ah,Ai,],Ai,] 
^ ^lE^' E [Ah,[Ai2,Ai,]] 



12d{F) 
1 



I (Il,l2)<l3 



12d(^)^ [E^^E [[Ai„Ai,],Ai,] 
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For doing explicit calculations with these expressions, it is convenient to split 
uj^^^ up into its parts in the following way: 

^{4) = _^(4A)__^^_J4F)^ (2.24) 

Additionally, in order to calculate the expectation value of the Wilson loop, 
one needs a perturbative expansion for the action. Expectation values are 
calculated in the usual way: 

where S is Wilson's action (|1.4)) and the integral runs over all elements U of the 
gauge group. Carrying out the Faddeev-Popov gauge fixing and rewriting the 
integration measure using the vector potential A^, one arrives at the following 
perturbative expansion (see, for example, |16j): 

S,ff = + 505(1) ^ ^2^(2) ^ g2sp^ ^ g2sgl,^ ^ Q^gZ^y (2.26) 

5(1) 

gives the three-gluon, S^"^^ the four-gluon vertex; S^^^^ comes from the 
Fadeev-Popov determinant, and S^"^^""^^ from the transformation of the inte- 
gration measure. S^'^^ is the term of order g^, which is quadratic in the gauge 
fields: 

BZ 

Here, a is the usual gauge parameter; for all the results summarized below, the 
Feynman gauge a = 1 had been used. The subscript BZ denotes the region of 
integration: the first Brillouin zone, all components of p run from — vr to +-7r. 
Additionally, the following standard abbreviations were introduced above: 

= 2sin(p^/2) 

= jZPl- (2.28) 

Note also that the sum over the Lorentz indices /i and v was explicitly written 
out, i. e. the usual sum convention is not used in this work: Repeated indices 
only have to he summed if this is explicitly denoted. 
Now the expectation value of an operator becomes: 

jDAOe-^^ff jDAOe-^'-°\^ 

=: <0>o+0(5o), (2.29) 
and the free gluon propagator is given by 

6 - (l - -] ^'f" 
<A^{p)A^{q)>0 = S^^{27r)'5{p + q) ^ °^ ^ (2.30) 



p^ 



S"^^ {2tt)'^6{p + q)-^ (Feynman gauge) 
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The explicit expressions for the other coefficients in the expansion of the 
action can also be found partly in j39j (but notice that the action is defined 
there with just the opposite sign in comparison to the convention here!); for the 
three- and four-gluon vertices, see ISH] or |37j. A more modern version, better 
suited for calculations, is given in 16_ • 

2.2.1 Leading order 

Using the formula and the expansion of the Wilson loop outlined above, 
one gets in leading order: 

V = g^ Urn l< tj(2) >o ^^fej^, (2.31) 

T^oo T 



with 



< >o (2.32) 



BZ 



(2vr)4 p2 lsin2(p^/2) sin2(p^/2) 



For a detailed derivation of this expression, see appendix El As mentioned 
above, subtr means that the self-energy contributions have to be subtracted, 
and I' represent here the spatial respectively temporal direction of the Wilson 
loop; these are different, fixed indizes and are not summed over. Because of 
the symmetry of the lattice, these two indices can be chosen arbitrarily without 
changing the result. 

In this order, the potential can be calculated exactly. First, perform the 
limit T oo. The following relations hold: 

1-^ 1 sm\p^R/2)sm\p,f/2) ^ ^ 
f^oo f sin2(p^/2) 



1_ sin\p^R/2)sin\p ^T/2) 
f-^^o f siv?{pu/2) 



lim — „ — — = 2Tr6{pu 

With this one gets: 



f = 2C.(F)sg/ 1^322!%^. (2.33) 



where 



hence 



BZ 



^ \ " -2 



p 



V 



2C2{F)gl f d^p sm\p^R/2) 



BZ 



(27r)3 



2C2(F)5, 



P 

2 , dV sin^(p;,i?/2) 



TV /a ^ 

° (2vr)3-^^ sin2(pW2)' 

-TT/a 
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In the continuum limit a —> 0, the integral gets its main contributions for 
p')^a ~ 0. Hence one can expand the sines in the denominator: sin^ (p^a/2)) ~ 
^p'^a^. The integration limits ^ir/a go to ±00 in the limit a — > 0. Therefore 
one gets: 

—00 

Now use the relation 

/■oc 

dte 



L 



00 1 

tx , 



X 

with this one can write: 

00 „ ^ 

V = 2C,{F)gl j Ap dtsin'{p'^R/2)e-'^^^^''^ 



12 



00 



2C,{F)glj^dt J ^sin\p'^R/2)e-^^'.' J ^e"-^ 



2 



-00 

00 



—00 
—00 

After completing the square in the exponent, the integral over p'^j^ can be done: 

Only the second term depends on R; the first gives an infinite contribution 
which is independent of R — hence it is the self energy contribution (this can 
also be seen by a close inspection of the graphs contributing to < w^^^ >o). 
After subtracting it, what remains is: 

V = ^Co(F) Ae-«V4t 



Att^/'^R 



Hence indeed one recovers the correct formula for the potential in leading 
order from lattice perturbation theory. 

2.2.2 Next-to-leading order 

In next-to-leading order, there are several different contributions to the expec- 
tation value of the Wilson loop: 



<W>=l-g^o... + gl< a;(3)5(i) >eo„„ -g^ < a;^^) >o -g^Wyp + 0{gf 



0^ 

(2.36) 
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(note the different sign for the w^^^-term, compared to [331) where the subscript 
"conn" means that one has to consider only the connected graphs contributing 
to the expectation value, and Wyp stands for the contributions to < > 
coming from the vacuum polarization graphs: 



Wvp 



(2) ^(2) 



^ ^ '-'meas -^conn 



< q(2) 

^ 'Jpp ^conn 

(2.37) 



2C2{F) 



d^p sin2(p^i?/2)sin2(p^f/2) 



BZ 



(^2)2 



+ 



Again, the signs of three of the terms are opposite to the ones in 
of the different sign convention for the action. 

The explicit contributions are, given here for future reference: 

<^ UJ^ ^ ^ ^ conn 

= -C2{G)C2{F) 



(2.38) 
because 



d^p 



d^k 



1 



d Zj D Zj 



+4 



' sin2(p^i?/2) 
sin(p^/2) 

sin(p^f /2) sin((2p + fc)^/2) 

sin(pi./2) sin(fc^/2) sin((p + k)yl2) sin(A;,,/2) sin((p + k)yl2) ^ 
sin((A: - p) J2) cos((p + A;)^/2) 



sin(p^f /2) cos(p^/2) sin((2p + A:)^/2) 

sin((2p + k)yfl2) 



sin(p^/2) sin(fc^/2) sin((p + /c)^/2) 



• siii{p^R/2) sm{p^f/2) sm{{p + k)^f/2) sm{k^R/2) 
■ cos((p + k)^R/2) cos{k^f/2) \ + \{fi,R)^ {u, f) \ 



(2.39) 



Here again the same comment with respect to the sign as above applies. 

For the calculation of < w^^^ >, use the decomposition (|2.24j) and determine 
each of the six contributions separately j4Uj : 



2C2{Ff - ^C2{G)C2{F) 



d^p s\i\^{p^R/2)sin^{kyf/2) 



(2.40) 
\ n 2 



■BZ 



{2nY 



sin2(p^; 



+ 



sin^(pi.) 



< UJ 



d^k 



BZ BZ 
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' sm'^{Pf,R/2)sm'^{kuf/2) 
sin2(p^/2)sin2(fc^/2) 

(2^ ^ (27r)4p2^2 



+ {(//,i?)^(l/,f)} 

lsin^((p^ + fc^)^/2) 
2 sin^(p^/2) 



sin2((p, + /,v)r/2) 



sm((jv - Av,)i?/2) sin(0v-Av,)/2) "- 
sin((p^ + A;^)i2/2) sin((p^ + A;^)/2) _ 



+2 



sin2(A;^/2) 

sm2((A;^ - p^)f /2) sm^p^R/2) sin^( fc^^/2) j 



sin2(p^/2)sin2(A;^/2) 
+ {(//,i?)^(i/,f)} 



(2.41) 



2 ; .V ..u y ^27r)4 J (27r)4^2^2 

sin(p,/2) 2sm(p.r/2)(S,-S2) 

sin(p^^/2) sin2(p^f /2) sin(A;^^/2) cos(A;^f ) 
sm(p^/2) sin(A;^/2) 

, sin2(p^i?/2)sin2(p^r/2) 



sm^(p^/2) 

sin(p^i?/2) sin2(p^f /2) sin(fc^i?/2) cos{k^f ) 

sin(p^/2) sin(/c^/2) 
sin^(y^i?/2)sin2(ff,r/2) 



Si?(/c/„-fc^)l + {(/x,i?)^(z/,f)} 



sin2(p^/2) 

sm'^{p^ R/2) sin^(p ^T/2) 
sin2(p^/2) 

sin^(p^.R/2) sin^(y^f /2) sin^(fc^^/2) cos(fc^f) 
sin2(p^/2) sin2(A;^/2) 



(2.42) 



BZ BZ ^ 

[ sin(p^/2) 

(Sr(0, -Pi,) + El - S2 - SR(fc^ - p^, -A;^)) 
sin(p^i?/2) . 2 



sin^(p^T/2) 



sin(;j^/2) 

(Sr(-p^, 0) + El - E2 - Sr(-A^, A^ - Pi,)) 
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sm'^iKf/2) 



sm{p^R/2) . 2 



sin2(p^/2) 



sin2(p^T/2) 



sin(A:i//2) 

sin^(fc^fi/2) 
sin2(A;^/2) 



(l + cos{kyf 



sm(p^/2) 



sin2(p^r/2) 



sin(j3^i?/2) 
sin(p^/2) 

sm{plR/2) 
' sin2(p^/2) 



sin(A:^/2) 
sin(fc^i?/2) 



cos{Kf) 



sin2(p^r/2) 



sin(A;^/2) 



cos(/s,,r)Efl(-p^,fe^) 



sin^(p,,/2) I >- J 



sin2(p^^/2) 



sm'^{kf,R/2) 



sm^{pj2) 



sm\k^/2) 



_^ sin^(p,fi/2) ^i^2(^^^/2)ST(fc., -fc.) 

sm^(p^/2) 



sin2(p^i?/2) . 2 



sin2(p^/2) 



sin2(p^r/2) 



sm\k^R/2) 
sm^{k^/2) 



cos(k^f) 



sm'^(puR/2) .2. r^, ,sm^(kuf/2) 



(2.43) 



C2iG)C2{F)g^o 



(fp 



d^k 1 



(27r)4 J (27r)4p2^2 



1 sin(p^/2) 
sin(p^^/2) 



sin2(p,f/2)E«(-p^,0) 



sm 



sin^(p,.f /2)Si?(A;^, -p^ - k^ 



sin(p,i?/2), ^.^2(^^^/2)5,^(0, -p,) 



sin(p^/2) 



sin(p^i2/2) . 2. ^/oxv / . . N 

sin {p,,T/2)^r{-p^ - K/i) 



sin(p^/2 



sin2(p^i2/2)sin2(p,f/2) sin2(p^i2/2)sin2(j9,f/2) ^ 



sin2(p^/2) 
sm(p^^/2) ^.^2/ 



sin^(p^/2) 



I ^ . , , 2.„ ^ sm(fc^J^/2) sm((p^ + fc^)i^/2) 
^ sin(p^/2) ^^'^ / ^ sin(A;^/2) sin((p^ + A;^)/2) 



• cos{kuf) 
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\ if^, R) - f) I + 2^^^^!^f^ sm\p^f/2)f 
[ J sm (p^/2) 



(2.44) 



I I 

and finally, due to symmetry properties in the colour indices 

< >o= 0. 

the functions S are the parts of the following functions S which are even 

^ Pill f^U ^ kn. 



(2.45) 



Here 

in Pf, -Pf„kf, 



Si 



S2 = e 



R-3 R-2 

gipp(R-l)/2 j2 J2 

xi=0 X2=xi- 
R-3 R-2 

>^(i?-i)/2 j2 j2 

xi=0 X2=x\- 



R-2 R-1 

^-iPtiXi ^ikfj.{x2-X3) 



R-2 R-1 

E E 

=xr+l X3=X2+l 



XI- 

R-2 R-1 

^-iPtiX2^ikf^{xi-X3) 



R-2 R-1 

E E 

=xr+l xs=X2+l 



xi=0 X2=xi+1 

The explicit expressions are for R> 2 j4Uj : 

,si n(p^i?/2) 1 sin(p^(i? - 2)/2) 
/OA . sin(p^/2^ 



El = 



-(^-2; 



+ 4 



1 sin((fc^ - p^){R - 2)/2) cos(fc^(j? + l)/2 - p^) 
4 sin((A;^-p^)/2)sin2(A:^/2) 
1 sin(p^(i? - 2)/2) cos(p^ - /c, ' 
4 sin(p^/2)sin2(A:^/2) 



S2 



sin(p^/2) sin^ 
lsin(p^(i^-2)/2) cos(fc^-p^/2) 

^ '"^ sin(/c^/2) sin((p^ + k^)/2) 

cos((p^ + k^)R/2 + fc^/2 - p^) 



(2.46) 



n(PM/2) 



1 sin(/c^|^« — z 
"^4 sin2(A;^/2 



,XR-2)/2) cos((p^ + fc^), 
i2(A:^/2) sin((p^ 



+ ^/.)/2) 



sin 

+ k^)R/2 + fc^^ 
-1 sm vt^/^j sM{Pfj, + k^,)/2j 

1 sin((A;^ - p^)(j? - 2)/2) cos(fc^(j? + 2)/2 - p^/ 
4 sin(p^/2) sin(A;^/2) sin((A;^ - p^)/2) 

1 sin(fc^(^ - 2)/2) cos((p^ + k^)R/2 + k^- p^) 
sin(p^/2)sin2(A;^/2) 



For R < 2, both of these functions v 
sums. Additionally for R > 1 one has: 

Sij(p^, kf,) = Re tuip^, k^^ 

l sin(p^(j?-l)/2)sin(fc^j;/2 
2 sin(p„/2)sin(fc„/2) 



sin(p^/2)sin- ^ ^, ^ 

vanish because of the constraints 



(2.47) 
in the 



.{k ^R/2-p^/2) 
fc^)(i?-l)/2)sin( 



sin(p^/2) sin 
jin((p^ + k^,){R - l)/2) sin(p^/2) 



l sin((p^ 

'2 sin((p^ + /j^)/2) sin(A;^/2 

1 sin((p^ + k^)R/2) 
-hM0\ ^^-^^l 



1 sin(p^i?/2) sm{kf,R/2) 1 sin((p^ + k 

2 sin(p^/2) sin(A;^/2) 2 sin((p^ + a,^ 



k^.)/2) 
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For R = 1, this function also vanishes, again because of the constraints in the 
sums. Obviously the function is also even under ^ k^j^; therefore terms like 

in the expressions ()2.4H) to 1)2.44(1 vanish. 

Additionally, the part of S^j which is odd under —Pfn k^ —kfj_ 

appears also several times in the expressions above (in < w^^) >o and in < 
^(3)5(1) ii ig given by: 



(2.49) 



sin((ff^ + k^){R - l)/2) cos(p^/2) 
sin((p^ + kf,)/2) sin(A:^/2) 

cos{k^,R/2 - p^/2) sm{pf,{R - l)/2) 



sin(/c^/2) 

1 sin((p^ - fc^)/2) 
4 sin(p^/2) sin(/c^/2) 



sin(p^/2) 

sin((p^ + k^)R/2) sm{{p^ - k^)R/2) 
sin((p^ + k^)/2) sin((p^ - k^)/2) 



Obviously, OR{Pf_i, k^) is also odd under Pf_i k^^. 

Finally, the vacuum polarization is given by the sum of the following five con- 
tributions, where the fourth and the fifth are unique to lattice perturbation 
theory: 



jj^gluonloop 



yrgluontadpole 



Y^ghostloop 



jrghosttadpole 



jrmeasure 



-C2{G) 



d'^k 



1 



BZ 



(27r)« 



p^ [p + k 



25^^ cos2(A:^/2) (2p^A;)' 



+ {2k+p)) i^k + p)^ E cos2(p^/2) 



+2 (2f+ k) [k^p) ^ cos{{p^ + k^)/2) cos{k^/2) 



(2.50) 



- 2C2(F)--C2(G) vr-p^p, 



2d 



+ gC2(G)V 



{d + 3)Ao + I- - + 3 A0-— cos(p^ 



d 



2d 



- 7An - 



2d 



-C2(G) 



p 
d'^k 



(2.51) 



1 



BZ 



fp' [p + k 



2{p + k] kyCos{k^/2)cos{{pu + k^)/2) 



1 

1 

"12 



(2.52) 
(2.53) 
(2.54) 
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For later convenience, the results are given here for an arbitrary numbers of 
dimensions d. Additonally, the abbreviation 

Ao = / (2.55) 

BZ 

has been introduced (compare appendix EJ. 



It is not obvious how one could get an explicit expression for the potential 
from these eight-dimensional integrals in momentum space, reproducing the 
result (|2.12|) . One way is to evaluate the integrals numerically; Heller and 
Karsch did this already in their original paper ^39J and fitted their results to 
the following formula: 

Here the factor | comes from the quadratic Casimir operator C2{F) for SU{?i) in 
the fundamental representation, the factor comes from the first coefficient 
of the beta function, and Vseif denotes the contributions from the self energy. 
Comparing with (|2.12|) . one should get: 



M = 82.07 



respectively 

InM = 4.408. 

The results Heller and Karsch obtained for In M for i? ^ 1 are consistent with 
this prediction - a further verification that l|2.12|l indeed gives the right result 
for the quark-antiquark potential. 



2.2.3 Transversality of the gluonic vacuum polarization 

Using the results Heller and Karsch obtained for the vacuum polarization |39j . 
one can show that it is transversal even for finite lattice spacing. This is not 
crucial for this work, but an interesting result in itself. So far, the transversality 
of the vacuum polarization, calculated in lattice perturbation theory, only had 
been checked in the continuum limit |37j . 

The contribution proportional to C2{F) of the vacuum polarization is ob- 
viously transversal. For the contributions proportional to C2{G) it has to be 
shown that 

E»n^^ = 0, (2.56) 

which proves the transversality. In order to show this, look at the various 
contributions separately: 
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+ 



{Pi 



-I +P-- Ao-- 



1 



n Y\9luontadpole 



C2(0)9? 

if) 



/((l-d)A„ + l-jL) 



+ 



En ^ Tighosttadpole _ 
PlJ.Pl^^'-jXl/ 

ET) -n Timeasure _ 
PuPui-i-fxi, 

H,v 

Here the abbreviation 
I : 



v2 1 

4 



Ar 



1 

2d 



An 



p 



1 

2d 
2 1 . 



<:?2(G)5o' 



1 



2d 
(2.57) 



(2.58) 
(2.59) 

(2.60) 

(2.61) 



d^q S('^\p + k + q) 



BZ 



{2ttY J i2T:Y q^k'^ 
BZ 



was introduced. In the calculation, the following Ward identity for the three- 
gluon vertex was used: 

E^MrM<^A(p,g,A;) = igo{2-aY5^'^\p + q + k) (2.62) 

• cos{jpx/2){q^5^x - qAx) - cos(jp^/2){k^5j,x - Kk\ 

Another useful formula which was needed here is: 

pk cos((p + k)/2) = i (ip + kf -p^ - P) . (2.63) 

For the evaluation of the other integrals appearing in the derivation, see ap- 
pendix C. 

Adding up all these contributions, everything cancels, and one arrives indeed 
at the result that the vacuum polarization is transversal. Hence it can be written 
in the following form: 



^nu{p) = (s^ivP^ -p^,Pu) n(p). 



(2.64) 



Additionally, one should show that the function Hip) is regular at p = 

or alternatively that n^jy(O) = 0. For the contribution proportional to C2{F) 
this is again obvious; and again the contributions proportional to C2{G) are 
determined separately: 



j^gluontadpole^Q^ = C2{G)gl5^, 



1 



1 



1 



(c!-l)(-A„ + ^ 



(2.65) 
(2.66) 
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^(0) = C2(G)<7o'<5, 



yimeasure i 



1 
1 

12 



(2.67) 
(2.68) 
(2.69) 



Again, some of the integrals listed in appendix C were used. Adding everything, 
one indeed arrives at the claimed result 



n^.(o) = 0. 



(2.70) 



Chapter 3 

Action Sum Rule 



In order to derive the action sum rule 

V + = ^ lim < L{t) >gg_o, (3.1) 

oR go T^oo 

essentially three steps were needed: 

1. Taking the logarithmic derivative of this led to the identity (|1.7() . 

2. Using the scaling behaviour of the potential; the result was p.l2p . 

3. Taking the limit of large T and thereby restricting the sum over all pla- 
quettes to one fixed time slice. 

The explicit form of the potential up to next-to-leading order was already 
derived in chapter 2; its easy to see that it has indeed the right scaling behaviour. 
Hence only the first and the third step remain to be checked. As already pointed 
out, the first step gives essentially an identity, so that in principle a check for 
this is not needed. 

Nevertheless, it is instructive to see how this identity looks like on the level of 
Feynman graphs, and the results derived here will be helpful later. Additionally, 
it will turn out that this identity is closely connected with the gauge invariance 
of the Wilson loop, so that this invariance can also be checked perturbatively. 

The first section will give an introduction and contains the necessary ex- 
pansions; the second section will deal with the perturbative check of (|1.7|) up to 
next-to-leading order. Then I will point out the connections to the gauge invari- 
ance of the Wilson loop, and how this can be checked. Finally the restriction 
to one fixed time slice will be examined. 



3.1 Preliminaries 

The identity (|1.7)) . when inserting the definition (|1.8j) and the formula 
becomes 
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Obviously, this identity has to be satisfied even for all finite T: 

d - - < SW > 

- P^ln<WiR,T) >= ' -<S> . (3.3) 
dp <W> 

Carrying out the derivative and multiplying by < W > gives: 

- P-ir <W >=< SW > - < S ><W > . (3.4) 
dp 

Now insert the relation H1.5|) between the lattice coupling constant P and the 
coupling constant go{a) appearing in perturbative lattice calculations; this fi- 
nally leads to: 



5o7^ < W > = < SW > - < S ><W > = < SW >conn ■ (3.5) 



1. 

That is the formula which will be checked perturbatively in the next section. 



3.2 Perturbative check 

3.2.1 Leading order 

Using the expansions 

S = 5(°)+O(5o). (3.6) 

and 

iy = l-5oV2)+0(5|]), (3.7) 

(|3.5j) becomes 

<a;(2) >0 = <C^(2)5(°) >eonn. (3.8) 

The expectation value on the left hand side was already given in section Ti. 2. II 

< w^^) >o (3.9) 
d^p s\T?{\p^R)s{ii^{\p^f) 



^'^'^^^ J (2^ ^2 1 — 27irT + .^.2/1 

BZ 



(2^)4 K^^^^ihp^.) sin2(ip,)^ 



this result was derived using Feynman gauge a = 1. In the following, if it is 
not explicitly indicated otherwise, all expectation values will be calculated in 
that gauge. 

The easiest way to calculate the correlator on the right hand side of 1)3. 8|) is 
to compute the correlator between two gauge fields and the action first; this is 
equivalent to the insertion of the action into a gluon line: 

5 - ^1^" 

< A^{p)A^{q)S >eo„n= 5^^^^^—:^{2T:f5{p + q). (3.10) 

This results hold for every gauge parameter a — hence inserting the action into a 
gluon line with arbitrary gauge parameter simply gives a propagator in Landau 
gauge a = oo! Therefore the calculation of the correlator < w^^^S^'^^ >conn gives 
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exactly the same result as if one would calculate the expectation value < w^^^ >o 
using the propagator in Landau gauge. For convenience, this expectation value 
will now be calculated in an arbitrary gauge. 
First, a;(^) is given by: 



UJ 



(2) = -J— i \- 



E^i" ' (3-11) 



Ad{F) 

with the sum along the Wilson loop: 

R-l f-i 
Y^Af = E^^K + ZA) + E^^K + /ii? + Z^) (3.12) 

I 1=0 1=0 

R-l f-1 

- E A^ino + fLR + 0f-lfL)-J2 ^uino + i>T - W). 

1=0 1=0 



Expressing this in Fourier space, one gets: 

-——sm{p^R/2)sm{p^^,.,^ . 
(2vr)4 sm[pa/2) 



= I |^sin(p,i;/2).i,.(p.f/2).--. -^^'^);^;^-^-' , (3.13) 

' BZ 

where 



nc = no + (iR/2 + i'T/2 

is the center of the Wilson loop. 

Hence calculating the expectation value, using the propagator (|2.3UI) . the 
result is: 



I ^ f d^p f d 



i{p+q)nc 



d(F)^^y^ (2^)4^yj2^)4 sin(W2) 



^fil3 — ^ul3 
sin{q0/2^ 



2C2(F)E / 



sm{p^R/2) sin(p^r/2) sin(g^fl/2) sin(g^r/2) 
d'^p sm^{pf,R/2)s\n^{pyf/2) 



(27r)4 p 



2 



'//a - ^va (5^/3 - ^vP ( ^ - (l- 



sin{pa/2) sin{p/3/2) y^^" \ a) p^ 
d^p sin2(ip^i?)sin2(ip^f) 



2C2(F) 

BZ 
( 1 



(27r)4 



,2 



+ 



ysin^(ip^) sin2(ipi,)y 
= < a;(2) >p^^^^ . (3.14) 

Therefore the expectation value < u;*^^) >q does not depend on the gauge pa- 
rameter a, i. e., it is gauge invariant! Naively one could have expected this 
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result, because the Wilson loop is defined in a gauge- invariant way. But gauge 
transformations on the lattice are defined with respect to the link variables, 
and it is not entirely clear what effects this can have in a weak-coupling expan- 
sion. Hence it is nice to see that even in the limit of small coupling, the gauge 
invariance is preserved. The close connection between the action sum rule and 
the gauge invariance of the Wilson loop which becomes apparent here will be 
discussed further in section ESI 

Additionally, using this gauge invariance, it is obvious that indeed 

holds — hence the action sum rule is valid in leading order. 

This validity as well as the gauge invariance of the Wilson loop in leading 
order can also be expressed in diagrammatical form in the following way: 




In these diagrams, a sum of the end points of the gluon line along the Wilson 
loop as well as along the plaquette, and a sum over all possible positions and 
orientations of the plaquette is to be understood. In the third graph, the gauge 
parameter a can be arbitrary. Gluon lines without an explicit index a are in 
Feynman gauge. 

Another way to express this more concisely is to introduce the operator 

G:=IJ A^{p)p^p,A^{-p); (3.15) 

BZ 

then one can write: 

< A^{p)A''{q) >o,a=< A^{p)A''{q) >o - (^1 - < yl^(p)A^(g)G >,o„„ . 

(3.16) 

If one denotes the insertion of this operator into a gluon line by a cross, one 
can also express this diagrammatically: 

a f i\ 
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where again a gluon line without an exphcit index represents a propagator in 
Feynman gauge. 

Now, using this notation, the gauge invariance of the Wilson loop in leading 
order can be expressed in the following concise way: 



0, 



or again diagrammatically: 



(3.17) 




0. 



3.2.2 Next-to-leading order 

Taking the higher terms in the expansions of the Wilson loop 

1^ = 1 - 5oV^) - 5oV3) - <7oV^) + O(5o^) 
into account, one gets: 



2 d 



<w> 



(3.18) 



+25o^<a;(3)5(i) >g +o{gl). 



2g^ < u;^v 



(3.19) 



For the right side of ()3.5() . the higher terms in the expansion of the action 
have also to be used: 



S = S^'^+goSW+giS^'^ + 0{gl). 



The result is: 



<WS > - <W >< S > 
-gl < +^4 < 5(0)^(2)^(2) 



+50 < conn 

-5^ < 5(0)c^(4) >,o„„ +0{gl). (3.20) 




Inserting this into the formula H3.5|) and using the vahdity of the action sum 
rule in leading order, one arrives at the following formula which has to be shown 
to be true: 



< S'-^'S'-'iO'-' >conn "2 < y'^") ' >conn + < ^'f'^UJ'^^' > 



conn 



^FP > 

conn 



+3 < cj(2)5(2) >^^^„ _2 < cj(2) (5(1)) >^^„,„ +2 < J^^sPp 
+2 < J^^Sil, >,onn +3 < J^^^S^^^ >o -2 < c^(4) >o . (3.21) 

Now one can use that the insertion of S^^^ into a gluon line transforms the 
propagator into Landau gauge. For every graph, one has to count the numbers 
of gluons lines into which S^'^'^ can be inserted; then one gets the following 
simple results: 

< S-^O^OJ^^) q(2) ^ _ o ^ ^(2) r,(2) 

^ ;j uj L> conn,a=oo — ^ uv '^conn,a=oo 
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o 







Figure 3.2: Contributions to the r.h.s. of (|3.21 



< "^LO^ -'(S'^ >conn,a=oo- 

/ ^(0), ,(2) c(2) ^ 

^ tj Lu '-^pp ^ conn,a=oo~ 

< S^^^J'^KS^^'^ > - 2 < tj(2) c{2) > 

^ tj LU '-'meas conn,a=cx> — ^ ^ ^meas conn,a=oo 



4 < cj(2)fc(i))2 > 

^ ^ V'-' / conn,a.=oo 

9 ^ , ,(2) . 
z, ^ uy ^pp conn, a=oo 



^ tj LU ij ^ conn,a=oo 
r.n'n.n..n'=nci — 

2 < LO'^' > 



3 < 



> 



conn,a=oo 



(3.22) 



Inserting these expressions into ()3.21|) , one sees that indeed the left hand side 
is equal to the right hand side. Hence if one calculates every graph in Landau 
gauge, the action sum rule is obviously valid up to next-to-leading order. But 
this is not entirely satisfying, because the results for the potential in chapter 
2, taken from 35], were calculated in Feynman gauge instead of Landau gauge; 
and it will turn out later that checking the energy sum rule can also be done 
much more easily in Feynman gauge. 

Hence what is needed is either a check that the expectation value of the 
Wilson loop is gauge invariant up to next-to-leading order, or an explicit check 



34 



CHAPTER 3. ACTION SUM RULE 



of H3.21() in the Feynman gauge. As will be explained in section 1531 these two 
are closely related; checking ()3.21|) will provide more than half of the check of 
the gauge invariance already. 

Therefore now 1)3.21(1 has to be checked. The relevant Feynman diagrams 
contributing to the left and right hand side, respectively, are given in the figures 
13.11 and 13.21 There are essentially three types of graphs: vacuum polarization 
graphs, one with a three-gluon vertex and one with two independent gluon 
lines. For the last two, the end points of the gluon lines are not summed over 
the entire Wilson loop; there are constraints: see the explicit form of to^^^ and 
(^(4) given in section 0^ In contrast, in the vacuum polarization graphs both 
gluon end points are summed over the complete Wilson loop. Additionally, 
in all of the graphs, a sum over all possible positions and orientations of the 
plaquette is to be understood. 

The vacuum polarization graphs can be grouped into two types again: the 
ones where the action is inserted into an internal gluon line (the second and 
fourth graph in figure 13.11) , and the ones where it is inserted into an external 
line (the first, third, fifth, sixth and seventh graph). The second class includes 
a gluon loop, a gluon tadpole, a ghost loop, a ghost tadpole and a contribution 
coming from the integration measure; the latter two are special for lattice per- 
turbation theory. These graphs can be dealt with easily if one uses the following 
observation: 

^ < AA;^(X)SW >conn= E < Al^^i^) >0,a ■ (3.23) 
I I 

This means that even if only one end of the gluon line into which the action is 
inserted is summed over the whole Wilson loop, while the other end is entirely 
arbitrary, the result will be the same as if one would have used a gluon line 
in an arbitrary gauge, but without the insertion. This can also be expressed 
diagrammatically: 



Using this, it becomes obvious that 

< S'^'^^tJ^^^S'*-^-' >conn = 2 < uj^"^^ S^'^^ >conn + insertions into internal lines 
< S^^^J^\s'''^^f >conn = 2 < a;(2)(5(i))2 >^^^^ + ins. into internal lines 



^ q{0), ,(2) 0(2) ^ - 9 <r , ,(2) ^(2) ^ 



<S^^^J^^S(^L>,onn = 2<J^)S^±,>,onn (3.24) 

or, again with diagrams: 
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and the same for the graphs with ghost hnes and the graph with the insertion 

of the contribution of the integration measure. 

Then there are the graphs where the action is inserted into internal gluon 
lines of the vacuum polarization. Taking them also into account, one can write: 

< 5(0)a;(2)5(2) >conn= 3 < S^'^u^'^ >,onn +A < S^^)u^^) >, (3.25) 



where the A-terms come from the operator G contained in the action. The 
explicit results are: 

A < (5(^)2a;(2) > 

(27r) 



sin(pa/2) sin(p^/2) 



+ 5o^cos(fcQ/2)cos(%/2)((2]r+ A;)^)^ 

+ kakp cos{{k + p)a/2) cos((p + k)f3/2) [k^pf 

~ XI ^t(^ ~p)i cos(p^/2)[(p + 2A;)afc^cos((p + k)p/2) + a (5] 

7 

- fc(2p + A;)[(p + 2A;)a^;3Cos(%/2)cos(p/3/2) + a ^ f3] 
+ k{2-f+k)[{k^p)oc cos{ka/2)kp cos((p + k)p/2) + a ^ p] 



(3.26) 
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and 



A < u;(2)^(2) > 

-C2{G)C2iF)J2 ^ '^"^ sm^{p^R/2)sm^{p,T/2)6, 



Sua S/i/B - Sul3 



a^BZ 



(2vr) 



sin(pa/2) sm(p_a/2) 



BZ 

2 



(fk 1 



Sapi^Y. Pp cos\k^/2) -icos\k^/2) J2 sm''ip,/2)k 



s\n\kal2)p^ + 2ki cos\k^/2) sin^(p„/2) 



[4 



+ P<,p^|-sin^(A:„/2) + - sm^(A;/3/2) - 2sin^(A:„/2) sin^(A;^/2) 



(3.27) 



For the graph with the three-gluon vertex (the so-cahed "spider graph"), 
one can make use of the fact that in w*-^^ , one of the end points of the three gluon 
lines is summed over the entire gluon loop, without any constraints. Therefore, 
according to (|3.23|) . an insertion of the action into this line just reproduces the 
normal spider graph without the insertion, and one has to consider only the 
effect of inserting the action into the other two lines. The result is, the A-term 
again coming from the operator G contained in the action: 



< 



5(o)^(3)_5(i) 



(3.28) 



with 



A < 0.(3)^(1) 



> 



-2C2{G)C2{F) 



d^p f d^q sm{pf,R/2)sm{p^T/2) 



(27r)4 J (27r)4 p^q^p'^qf 
BZ BZ 

1 ( cos{q^/2) ^^^^ „ 

sin((p, + q,)/2) l,sin(p,/2) + ^^^'^^ 

^ g^(g'+»^cos(g^/2)cos(p^/2) ^ cos(g^/ 2) cos((gj. + Pv)/2)qyp^ 



sm.{py/2) 
sin(p^r/2)sin((p^ + 2g^)/2) 



sin(p,,/2) 
sin((p^ + 2q^)R/2) sm{p^R/2) 



sin((p^ + 2q^)/2) sin(p^/2) 
+2 sin(g^r/2) sin((p^ + q^,)R/2) cos((p^ + qu)T/2) cos{q^R/2) ] \ 



+ {if,,R)^{i^,T) 



(3.29) 



For calculating the effect of inserting the action into the graph with two 
independent gluon lines, it is convenient to split up a;*-'*^ as shown already in 
1)2.24(1 and treat the various terms separately. In w^^"^), all end points of gluon 
lines are summed over the entire Wilson loop, hence using H3.23() . one gets: 



(3.30) 
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The other contributions are much more compHcated. In u;^^*^^ to uj^'^^\ at least 
one of the end points of the gluon lines is summed over the entire Wilson loop, 
and therefore (|3.23|) can be used to simplify the calculation: insertions of the 
action in that gluon line simply reproduce the same result as for the graph 
without any insertions. In uj^^^\ even this simplification is not possible. 
The end result is, where again the A-terms come from the insertion of G 



(3.31) 



respectively 



with: 



< -'S''- >conn 

r.nnn. 



2 < J'^^ >o +A < J^^^ >; 



(3.32) 



A < 0.(4^) > 

C2{G)C2{F) 4 

Z 9o 



BZ 



(2^ J (27r)4p2^2 
BZ 



^ A^\^i^sm\p^T/2)cos\k^T/2)+sm\k^T/2)cos\p^T/2) 
sm\p^R/2) sm^{k^R/2) 



+ 



sin2(p^/2) sin2(A:^/2) 

sin^{k^R/2) sm\p^T/2 ) ^ lsm'^{{p,, + k^)T/2)sm'^{{p^ - k^)/2) 
sm\k^/2) sm^{pj2) 



+ 



sin((p^ + k^)R/2) sin((p^ - k^)R/2) 
sin((p^ + k^)/2) sin((p^ - k^)/2) 

kfX ( sm^{pf,R/2) sm\kf,R/2) sm\k^T/2) 
fc2 1^ sin2(p^/2) sin(A:^/2) sin(A;^/2) 

1 sin((p^ + k^)T/2) sm{{pf, - k^,)/2) 

2 sin(p^/2)sin(fc^,/2) 

sin((p^ + k^,)R/2) sin((p^ - k^)R/2) 
sin((p^ + k^)/2) sin((p^ - k^)/2) 



sin2(p^/2)sin2(A:^/2) 

2 \ 



cos2(p^r/2) 




{^i,R) ^ {u,T) 



(3.33) 



A < > 

C2{G)C2{F) 4 
'a fo 



d'^p 



d^k 



(27r)4 J (27r)4p2fc2 
BZ BZ ^ 



' kl sm{p^R/2) 
fe2 sin(p^/2) 



sin2(p^r/2)(Si -S2) 



sm\p^R/2) sm^ip^T/2) 
sin'^{pa/2) 



f sm\k^R/2) 
'^2 I sin2(A:^/2) 



R 
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7^ — sm [pyl 1) — — — 

/ sin(p^(i?-l)/2) cos(p^/2) 
y sm(p^/2) sm(A;^/2) 
cos{k^R/2 + sin((p^ + k^){R - l)/2) 

sin(A:^/2) sin((p,, + k^)/2) 

sinikfAR - l)/2) siu(A,,(i? + l)/2) sin(p^i2/2) 



+2 



sin(A;^/2) sin(p^/2) 
cos(fe^(i? + l)/2) sm{{p^ + A;^)(i? - l)/2) \ 1 



sin(p^/2) 



sin((p^ + kf,)/2) 



+ {ili,R)^iiy,T) 



^0-2/ ^,,, sin2(p.T/2) fcg sm^(fc^i?/2) 
+2sin^(p^.R/2) 2/,„ /ON „:,^2^r /o\ cos(fc^T) 



+2siii2(p^r/2) 



sin^(p^/2) A:2 sm^(fc^/2) 
sm^{Pf,R/2) kl ( sin^(fc^r/2) 
sin2(p^/2) P \ sin2(fc^/2) 



-2 



sin^(p^i?/2) . 2. rr^in^ h^^u sm^{kf,R/2) sm^(fc^T/2) 
sin /2; ^.^^^^^^^2) sm{kj2) 



sin (p^/2 

+2sin2(pui?/2)^ , . 

^ ^ sin(p^/2) sin(A;/,/2) 

/ sin(p^r/2) sin(A:^r/2) _ sin((p^ + fe^)r/2) \ 
' V sin(p^/2) sin(A;^/2) sin((p^ + A;^)/2) J 

A < a;(4^) > 



A;2 



C2(G)C2(F) 4 



12 



(27r)^ 



BZ 



k^ Hm{p^R/2) 



4^ 



fe2 sm{p^/2) 



(27r)4p2^2 

^ _^ l sin(p^i?/2) l sm{k^R/2)sm{{p^ + k^)R/2) \ 
' ' 2 sin(p^/2) 2 sin(fc^/2) sin((p^ + A;^)/2) J 



(3.34) 



sm^{puT/2) 



sm^{p^R/2) sin^(p^r/2) 
sin2(p„/2) 



^2 (sm\k^R/2) 



fc2 \^ sin^(A;^/2) 



.2,^ D/o^sin^(P^^/2) sin2(A;^i?/2) kl 



+2sm'{p^R/2)- 



+2sin2(p^r/2) 



-2sin2(p^r/2) 



+2sin2(p^r/2) 



sin2(p^/2) sin2(A;^/2) p 
sin^(p^i?/2) sin^(fc^i?/2)fc^ 



sin^(p^/2) sin^(A;^/2) k^ 



(1 + cos{kuT)) 



sin(p^i?/2) sin(fc^i?/2) k^ 
sin(p^/2) sin(A:^/2) fc2 

sin^(p^i;/2) sin^(fc^i?/2) sin2(fc^r/2) k^K 
sin2(p^/2) sin(fe^/2) sin(fe^/2) k"^ 
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-sin2(p^r/2) 



sin(p^fi/2) sin^(fc^r/2) k^K ( ^sm{p^,{R - I) /2) cos(j>^/2) 



sin(p/,/2) sm(fc^/2) fc2 



sin(p^/2) sin(A;^/2) 



^ cos{k^R/2+p^/2) sm((p^ + - l)/2) 

sm(/c^/2) sin((p/, + /c^)/2) 

s\n{k^{R - l)/2) sin(fc^(i^ + l)/2) sin(p^i?/2) 
sin(A;^/2) sin(p^/2) 
cos(fc/,(i2 + l)/2) sm((p^ + A;^)(i2 - l)/2)' ~ 



sin(p^/2) sin((p^ + k^)/2) 

.2,^ j,i^.^HPvT/2)sin^{k^R/2) 



-2s\n\p^R/2) — 



sin(fc,.T/2) 



kiikj. 



sin(p^/2) sin(A:;,/2) — ' ' p 

' sm(p^r/2) sin(/g^T/2) _ sin((p^ + /c^)T/2) \ 
^ sin(p^/2) sm{k„/2) sin((p^ + A:^)/2) J 



-2sin2(p^r/2) 



sm^(p^i?/2) fcg / sin^(fc,r/2) 
sm2(p^/2) ^2 1^ sin2(fe^/2) 



, ^ sin^(p.T/2)sin^(p^i^/2) ^^ 
sin^(p„/2) 



Sau) 



sin2(A:^i?/2) sin2(A;j,r/2) fe^fcj, 
sin(A;^/2) sin(A;^/2) p~ 



.2,^ D/oNsin'(p.T/2) sm\k^R/2) kf, 



+4sin^(p^i?/2)- 



sin2(p^/2) sin2(fe^/2) k^ 



+4sin2(p^r/2) 



-2sin^(p^i?/2)^ 



sin^(p^i?/2) sm^(fc,T/2) 
siii2(p^/2) sin2(fe^/2) 

sin2(p^r/2) sin2(A:/,i?/2) 



cos(/c/,i?)^ 
cos(fc.T)-^ 



+2sin2(p^r/2) 



sin^(p^/2) sin^(fc^/2) 

siu'{p,M/2) siu-'(A-,,i?/2) siu2(Avr/2) t^.A,, 



sin2(p^/2) sin(A;^/2) sin(/e^/2) p 



(3.35) 



A < ^(4^) > 

C2(G)C2(F) 4 

6 



' sin(p^i?/2) 
sin(p^/2) 



sin2(p^r/2) 



^l sin(p^i^/2) fc^ sm{k^R/2) sm((p^ + fc^)fi/2) 

^4 sin(p^/2) p sin(A;^/2) sin((p^ + A:^)/2) 



•(cos(A;^T) - 1) + 



k^h sm2(fc.T/2) . ^^^^ sin((p^ + A:^)i?/2) ^ 
fe2 sin(/c./2) ^^"^^A^^/^^ sin((p^ + k,)/2) } 



1^^^ ^ 7 2 ^ sin2(p^/2) k^ sm(j>,/2) 



■ sin\p,R/2f^^^^^J^ sin(fc.r/2) 



sin(A;^/2) 



sin((p^ + fc^)r/2) 
sin((p^ + fe^)/2) 



(3.36) 



and finally 



A < >= 0, (3.37) 

again because of symmetry properties in the colour indices. The functions Ei, 
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T,2 and S/j can be found in section r2.2.21 
Putting everything together, the result is: 



/ C2{G)C2{F) 
I 6 

BZ 

l sin(p^(i?-2)/2) \ 
2 sin(p^/2) i 



+C2iG)C2{F) 



d'^p f d^k 1 sin2(p^i?/2)sin2(p^r/2) 



(27r)4 7 (2^)4 ^2^2 sin2(p^/2) 



(kl sm^iKT/2) Ikl sui\k^R/2) \ 
\P sin2(/c^/2) 6^2 sin2(/c^/2) j 

2C2(G)C72(F) [ d^p f d^k 1 sin2(p^i?/2)sin2(p^r/2) 



(2^)4 7 (2^)4 p2 ^4 sin2(p^/2) 



1 A;' 



.sin2(A:^i?/2)sin2(A:,r/2) . . 

ysm (A;;^/2) 2sm (A;^/2) 

C^2(G)C2(F) /■ d^p f d^k 1 fc^ sin^(p^i?/2)sin^(fc^i?/2) 



2 7 (27r)4) 7 (27r)4) p2^2 ^2 sin2(p^/2)sin2(A:^/2) 



■s\n\k^T/2) 

C2{G)C2{F) f d^p f d^k 1 kl sm\{p^-k^)/2) 



9 7 (27r)4 7 (27r)4^2^2^2sin2(p^/2)sin2(A:^/2) 



•sin2((p^ + A;^)r/2) 



' sin((p^ + k^)R/2) sm{{p^ - k^)R/2) \ ' 



sin((p^ + k^)/2) sin((p^ - k^)/2) J 
C2{G)C2{F) f d^p f d^k 1 A:^sin2(p^r/2)sin2(A:^i?/2) 



2 7 (27r)4 7 (27r)4p2^2^2 sin2(p^/2) sin2(A;^/2) 

C2(G)C2(F) f d^p f d'k 1 H^-^^2^p^j./2)^Hp^.Rm 



(27r)4) 7 (2^)4) p2^2 ^2 ^-"^ ' ^ sin(p^/2) 
BZ BZ 

sm{k^R/2) sm{{p^ + k^)R/2) 
' sin(/c^/2) sin((p^ + kf,)/2) 

I C2(G)C2(F) r A / rf^fe \ ^'cin^(p,r/2) ^^"^^^^/^^ 



+ 



(27r)4 7 (2^)4 p2p ^2 ^^-^ ' ' sin(p^/2) 

BZ 

( sm{pf,{R - 2)/2) cos(p^ + k^) 
^ sin(p^/2)sin2(/c^/2) 
sin((ff^ + k^){R - 2)/2) cos(fc^(i? + l)/2 + p^J 
sin((p^ + A:^)/2)sin2(A;^/2) 
sin(p^(ii - 2)/2) cos(p^/2 - k^) cos{{p^ + kf,)R/2 + k^/2 - p^) 



sin(p^/2) sin(A:^/2) sin((p^ + k^)/2) sin((p^ + k^)/2) 
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sm{k^{R - 2)/2) ^ sin((p^ + k^){R - 2)/2) cos{k^{R + 2)/2 +p^/2) 



sin^(A:^/2) 



sin(p^/2) sin(A;^/2) sm((ji^ + k^)/2) 



sm{k^{R - 2)/2) cos((p^ - A;^)i?/2 - A;^ - p^) 



+C2(G)C2(F) 



sm(p^/2) sm2(A:^/2) 

d^p /■ d'^k 1 k^ky 



(27r)4 y (2^)4 p2^2 p 



sm^(p^i?/2) 
sm2(p^/2) 



cos2(p^r/2) 



s\n^{k^R/2) sm2(fc^T/2) 1 sin((p^ + k^)T/2) sm((p^ - fc/,)/2) 



sin(/c^/2) sm{k^/2) 2 
/ sin((p^ + k^)R/2) sin((pp 



y sm((p^ + k^)/2) sm((p^ - A;^)/2) y 
sm{p^R/2) sm{k^T/2) 



sin(p^/2)sin(A;^/2) 
fc u)i?/2) \ 



sin(p^/2) sm{ki,/2) 



■ cos{pyT/2) cos{k^R/2) 



+- 



C2{G)C2{F) f (fp f d^k 1 k^k^ 



6 



BZ 



{2nY 



BZ 



(27r)4p2^2 ^2 



■sm2(p^r/2) 



sm(p^i?/2) 
sm(p^/2) 



■ rA- p/gN sin^(W2) / sin(p^i?/2) sin(/c^fi/2) sm{{p^ + k^)R/2) \ 
■''''^ / ^ sin(A:,/2) [ smip^/2) sm{k^/2) + sin((p^ + A;^)/2) ^ 

C2{G)C2{F) f d'^p 



12 



(2vr)^ 



(27r)4p2^2 p 



1 k^ku . 2, ^/o\Sin(PMi?/2) 



sm(p^/2) 



-5 
-7 
+7 



sin^(fc,r/2) / sin(p^(i?-l)/2) cos(p^/2) 
sin(/c;,/2) y sin(p^/2) sin(A;^/2) 
cos(fc^/?/2 + p^/2) sm((p^ + k^){R - l)/2) 

sin(A:^/2) sin((p^ + A;/,)/2) 

sin(fc^(i; - l)/2) sm(fc^(i? + l)/2) sm(p^i?/2) 
sin(A:^/2) sin(p^/2) 

cos{k^{R + l)/2) sin((p^ + k^){R - l)/2) 



sin(p^/2) 



sm((p^ + /c^)/2) 



(3.38) 



Looking again at (|3.2H) . one sees now that 

A < 5(2)^(2) > -Ia < > +A < 5(1)^(3) > _^ ^ ^(4) >^ 

(3.39) 

has to be satisfied in order to show that (|3.5|) is true in next-to-leading order. 

The only way to show this is to numerically evaluate the various eight- 
dimensional integrals in momentum space for arbitrary R and T. Obviously 
this would take a lot of computer time; in the light of the fact that (|3.5j) is 
an identity and therefore in principle has not to be checked (as explained at 
the beginning of the chapter, these calculations are mainly done for illustrative 
purposes, and because the results will become useful later) this is not justifiable. 
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Therefore only a check for the very special case R = T = I was done. Most 
of the results above simplify significantly then (for example, the functions Si, 
E2, TiR and Or vanish); most of the integrations can even be done almost 
exactly, with only the one numerical constant Aq remaining. 

But the contributions from A < ^S"*^^^^ co^"^^ > and A < S^^^co^^^ >, which 
are now given by the following integrals 



A<(S«)" 



R=T=1 



BZ 



(27r)4 J (27r)4p2)2(p)2(^'^^)2 



BZ 
4 



{p + 2k)l J2 K cos'(Pp/2) + cos2(A;^/2) ((2|j + k)k 
+klcos\{p + k)^/2){k'^pf 

4 

+ 2k)f,k^ cos((p + k)^/2) ^ kp{k - p)p cos(pp/2) 
{2f+ k)k{{p V2k)^kf, cos(A;p/2) cos(p^/2) 



-{k - p)f, cos{kf,/2)kf, cos((p + k)p/2) 



- {p + 2k)^{p + 2k), J2 klcos\pp/2) 
\ p=i 

+k^K cos((p + k)^/2) cos((p + k)u/2){k^pf 
4 

- ^ kp{k-p)pCos{pp/2) 

■ \{p '+2k)p,K cos((p + k)^/2) + [p '+2k)^kp cos{{p + fc)p/2) 
-i2f+ k)k {p'+2k)p,ky cos{k^/2) cos{p^/2) 

+{p + 2k)ukp_ cos{kp/2) cos(p^/2) - {k - p)p cos{kf^/2)k^ cos((p + k),y/2) 



-{k - p)u cos{ku/2)ku cos((p + k)p/2) 
A < 5(1)0.(3) > 



\ppPu (3.40) 



d'^p 



R=T=1 
d^k 



BZ 



Pu cos{k^/2) ^{p + k)pkp 



(27r)4 ] (27r)4p2(^2)2(p+^)2 ^ 

+Pfj.kp{p + k),y cos{ky/2) cos(pp/2) + p^Jt^ cos{kp/2) cos{{p + k)^/2) 
■ku cos((p + k)u/2) cos{kp/2) > + < /i ^ z> 



(3.41) 



still have to evaluated completely numerically. This was done using the standard 
routine Vegas from the Numerical Recipes |41j . On the other hand, the integrals 
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which can be done by hand give: 



> 



A < a;(^) > 



R=f=i 192 ° 36 ° ^ ' ' 

- - = — Ao--A^ (3.43) 
R=T=i 192 " 6 ° ^ ^ 



Using the result of a numerical integration for Aq, one finally gets the fol- 
lowing results: 

A < 5^^)^^^) > = 0.00070077 + -0.00000005 
-^A < (^S^^^y u^^^ > = 0.00187522 + -0.00000145 

A < 5(^)cj(^) > = -0.00254364 + -0.00000265 

-A < J"^^ > = -0.00003397 + -0.00000062 (3.44) 

Add everything up: 

A < > -1a < (5(1))' ^(2) > ^ ^(1)^(3) > _^ < ^(4) > 

= -0.00000162 + -0.00000478, (3.45) 

hence in the range of the numerical error, ()3.39p is indeed satisfied — or, in 
other words, the action sum rule respectively the identity (|3.5|) is valid up to 
next-to-leading order (for R = T = 1). 



3.3 Gauge invariance of the Wilson Loop 

The gauge invariance of the expectation value of the Wilson loop in leading 
order has already been checked in section 13. 2. 11 For checking this in next-to- 
leading order, one makes use again of the relation ()3.16|) which connects the 
propagator in Feynman gauge with the propagator in an arbitrary gauge by 
using the insertion of the operator G. 

If one applies this relation to an arbitrary graph with n gluon lines in an 
arbitrary gauge a, it follows that one can write this graph as a sum of the 
following contributions: the same graph with all propagators in Feynman gauge, 
n graphs with n — 1 propagators in Feynman gauge and the insertion of G in 
the remaining line, n(n — l)/2 graphs with n — 2 propagators in Feynman gauge 
and insertions of G into the two remaining lines, and so on; the last one is a 
graph with insertions of G into all lines. The coefficients of the terms in this 
sum are (-1)™ 1-^ , where m is the number of insertions of G. 

Here is an example to illustrate this (n = 3); as usual a gluon line without 
an explicit index a represents a propagator in Feynman gauge, and a cross de- 
notes the insertion of the operator G (|3.15l) : 
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Now looking at the graphs contributing to < > in next-to-leading order 
(see figure IH^ . ones sees that there are up to four gluon lines in them. Hence 
one would expect that if one calculates < 1^ > in an arbitrary gauge, one gets 
a polynomial in (^1 — of degree four. If one can show that all terms in this 
polynomial with the exception of one of order zero (which is independent of a) 
vanish, then the gauge invariance of the expectation value of the Wilson loop 
is proven. 



Terms of order three and four 

The only contributions with a power of (^1 — greater than two come from 
the two vacuum polarization graphs (gluon tadpole and gluon loop) where G 
can be inserted into external as well as into internal lines, and from the spider 
graph. One can easily show that these contributions give zero. For this, use 
(|3.2,3j) : introducing the operator G again, that formula is equivalent to: 



E 

I 



< AiA^ix)G >, 



0; 



(3.46) 



This means that even if only one end of the gluon line into which G is inserted 
is summed over the whole Wilson loop, while the other end is entirely arbitrary, 
the contribution of the graph will vanish. Diagrammatically: 
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So, obviously the contributions from the vacuum polarization graphs where G 
is inserted into an external line, and from the spider graph where G is inserted 
into the line whose endpoints is summed over the whole loop, vanish, and only 

terms up to ( 1 — — ] remain in the polynomial which results from < W >. 



Terms of order one 

The contributions of order (^1 — were already given in section rj.2.21 - the A 
terms were just the contributions to < SW >conn which resulted from exactly 
one insertion of G into a gluon line. It also has been shown in that section that 
the condition that all of these terms add up to zero 1)3. 39() is equivalent to the 
validity of (|3.5|) in next-to-leading order. But as already stressed several times, 
(|3.5|) is an identity, hence one can conclude that the contributions of the graphs 
with exactly one G insertion indeed add up to zero. 



Terms of order two 



The only contributions which are left are the ones with two insertions of G. 
Using 1)3.46(1 . one sees that such contributions can only come from < lo^^^S^^^ >o, 

^conn- Fortunately, the explicit calculation of 
these three contributions is not too hard; the end result is: 



■sm\pf,R/2) sin2(p^f/2) 

.(kcos{{p + k)^/2)-k^cos{{p + k)j2)] . (3.47) 



Considering now that <W>is given by 

<W>= -^gl < (5(1))' >eo„„ -9t < u^^'^ >o +9^0 < ^(^^5(1) >o + . . . , 

(3.48) 

one sees that the three contributions of order ( 1 — — ) indeed add up to zero. 
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Hence it has been shown that in the polynomial in — — j of degree 4 one 
gets from calculating < 1^ > in an arbitrary gauge, only the first term (of order 

^1 — and hence independent of a) survives: The term of order 1 gives zero 
because of (|3.39p , which is equivalent to the identity H3.5() ; the term of order 2 
has been shown to vanish exactly above, and the terms of order 3 and 4 vanish 
because of (|3.46p . Therefore the proof is complete: up to next-to-leading order, 
the expectation value of the Wilson loop is indeed gauge invariant. 

3.4 Restriction to a fixed time slice 

What remains to be checked is the restriction to one fixed time slice, expressed 
by the equation 

lim i < 5 >qq-o= lim < L{t) >gg-_o . (3.49) 
Using H1.8() . this is equivalent to: 

hm ^ < SW >conn= lim < LW >conn ■ (3.50) 
f^oo T f^oo 

Alternatively, with 1)3. 5() . it can also be written as: 

d 1 

lim < LW >conn= lim <W > (3.51) 

f^oo CIQq f-*oo T 



The expansion of < LW >conn is very similar to the one for < SW >conn (!3.2n|) : 

<LW>-<L><W> 
= -gl < L(0)a;(2) +^0^ < L^'^^S^'^J'^ >,onn 

+5o' < L(0) 5(1)^(3) _g4 < ^(0)^(4) ^0(^6) (3 52) 

For the right hand side of (|3.51|1 . one can again use (|3.19|1 . 
3.4.1 Leading order 

On the right hand side of ()3.5ip . in leading order the only contribution comes 
from 



< >Q = (F) J 



BZ 
1 



d^p sin2(p3i?/2)sin2(p4f/2) 
(27r)4 p' 



n2 



+ -^7-^ ■ (3-53) 



sin^(p3/2) sin^(p4/2) 
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For convenience, the spatial direction of the Wilson loop, previously denoted 
simply by /i, has been chosen to be the X3 direction here, and euclidean time, 
previously denoted by 1^, has been identified with X4. 
Now the limit of large T has to be considered. Obviously 

^. 1 f d^p sin2(p4f/2)sin2(p3ii/2) 
lim — / - — — TT- -— = (J. 

f^ooT J (27r)4 p2 sin2(p3/2) 

BZ 

Only the second term can give a non-vanishing contribution; using 

^^^^ ^sm^(p,f/2) ^ 
f^ooT sin2(p4/2) 

(compare with section r2.2.1() . what remains is: 

r 1/ (2)^ on (j^^ [ sin^(p3i?/2) 

hm ;^ < ^ >o= 2C2{F) / ——3 ^2 (3-54) 

with 

P =2^Pr 



On the other hand, an insertion of L^'^^ leads to non-conservation of the fourth 
component of the momentum vector: 



<A^{p)A^{q)L(') >, 



i;AB 



p2q2 



-i{p+q)^/2 



(3.55) 



Making use of this, an explicit calculation yields the following result: 

= 2C2{F) / — —- / ^sm^(p3ii/2)sm(p4T/2)sm(g4r/2) 



BZ 



(2vr)4_y 2. ^2(^-?+^2 



(3.56) 



/ cos((p4 + g4)(ree,4 - t)) cos((p4 + g4)("-c,4 - t - 1/2)) \ 
1^ sin2(p3/2) sin(p4/2)sin(g4/2) j' 

where again ric is the center of the Wilson loop and hence 

= no,4 + ^f. 

In the limit of large T, the first term vanishes due to the fast oscillations of the 
two sines depending on T; the second gives only a non- vanishing contribution 
for p4 = q4 = 0: 



lim < 2C2(F) / 

BZ 



d^p sin^(p3^/2) 



(3.57) 
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which is identical to the result above. Hence in leading order the restriction to 
the fixed time slice t indeed works: 

lim i < S >qq-o= lim < L{t) >,g-_o +0{g^). (3.58) 



For the special case nc,4 = t (the fixed time slice lying in the middle of the 
Wilson loop), the result above simplifies to 

^ ±1 id ^conn 

l^^ll flflf^ '"'^'''^"^ 



(2t)' 

BZ -TT 



p. 



1 ^ cos((p, + ,,)/2)) \ 



y sin2(p3/2) sin(p4/2) sin(g4/2) ^ 

and it is obvious that the first term vanishes for all T, because the integrand 
is odd in p4 as well as in 54. For the cosine in the numerator of the second 
term, use the trigonometric relation cos((p4 + g4)/2)) = cos(p4/2) cos((74/2) — 
sin(p4/2) sin((74/2); the sines do not contribute because they would again give 
a function odd in ^4 as well as in (74. Hence what remains is: 

= 2C.(F) fp^, 1'-^^ ^jf + sin\p,R/2f^^P^^/'^ ^^^'^^^'^^ 



BZ 



(2^)4 J 2tt p2 (^^ + ^ ' ' tan(p4/2) tan(g4/2) ' 

an expression which in the limit T — 00 again gives 

lim < >_= 2C.(F) / (3.61) 



T^oo J (2vrj 

BZ ^ 

3.4.2 Next-to-leading order 

The graphs corresponding to the next-to-leading order contributions were mostly 
given in the figures ITT] and alreadv. but now the plaquettes are not summed 
over all possible positions, but only over the ones with a fixed time. Ad- 
ditionally, three extra graphs have to be considered now, corresponding to 
< >^„„„, < l(i)5(i)c^(2) >^„„„ and < L^J^) >^„„„ (see figure ESJ. 

There again the possible positions of the plaquette lie only on the fixed time 
slice. 

By inspecting the calculation of 5*^^^ and S"^^-* in JH] and modifying it to give 
explicit results for L^^^ and L^'^\ one sees that essentially the contributions of 
the three graphs above can be obtained by using the usual four-gluon vertex in 
conn and the usual three-gluon vertex in the two other expectation 
values. The only crucial differences are that one has to replace the usual four 
momentum conservation at these vertices with conservation of only the spatial 
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components of the momentum, and that some additional phase factors appear. 
This will be studied in more detail below. 

Now, as usual, consider the different types of graphs separately. There are: 

• vacuum polarization graphs 

• spider graphs 

• graphs with two independent gluon lines 

For simplicity, the Wilson loop will now be chosen to lie symmetrically to the 
fixed time slice, and only the special case t = will be treated. Thus one has 
no,4 = -^T. 



The vacuum polarization graphs 



In contrast to the situation in section there are now two additional vacuum 
polarization graphs (the first two displayed in figure 13. 3|) ; the other five are 
very similar to the ones already encountered in section but now instead of 
calculating the effect of inserting the operator G, one has to use (|3.55|) in the 
gluon lines. 

First, look at the cases where L is inserted into an external line (represented 
by the first, third, fifth, sixth and seventh graph in figure IXT]) . Because of the 
non-conservation of the time component of the momentum, the relation (|3.23|) 
can not be used any more; one has to calculate the contributions explicitly. The 
result is: 



An \- f d^P f d^q sin\p^R/2) sm{p^f/2) sin(g4T/2) 
^^'^^^ J J^^i (^ 



sin(pQ,/2) sin(q'^/2) 



1 
2 



BZ 



{P- 



-i{p+q)(t/2 



{27rf6ip + q). 



(3.62) 



In the limit T — > 00, only the terms with a = 7 = 4 give a non- vanishing 
contribution, and these only for p4 = (74 = 0. Hence the limit is, carrying out 




Figure 3.3: Additional contributions to < LW > 
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the integration over qi, q2 and (73 using the Delta- function: 
This is identical to 

2 lim —Wyp, 

where TVyp can be found in section 12.2.21 and, using (|3.24)) . therefore also 
identical to 



conn,ext 



lim U-< S('W'^S('^ >_n,e.* + < 5Wc.(2)l (5(1))' > 

f~>oo \ 2 V / 

— < "^UJ^ Spp > conn,ext — < 5*^ ^UJ^ ^ ^meas ^ conn,ext ^ ; (3.64) 

where only the insertions of S^^^ into the external gluon lines are considered. 

Next, look at the graphs where the sum over the plaquettes on the time slice 
t = is inserted into an internal line, represented by the second and fourth 
graph in figure 13.11 and at the two graphs with insertions of L^^^ respectively 
L^"^^ depicted in figure ESI They give the following contribution: 

(F\\^ f J sin^fe^/2) sin(p4f /2) sin(g4f /2) 



sin(p,/2) sin((7^/2) "'^ ^ 



nL(p,P4,g4), (3.65) 



where 11^ represents the vacuum polarization tensor with L^^^ inserted into one 
of its internal lines respectively the two contributions with L^^^ or L^'^\ In the 
limit T oo, this reduces to: 

^Cm^ I ^'-^^^^^TiiMO,0). (3.06) 

°'Pbz ^ ' [p) 

Now first consider the contribution from the second graph in figure 13.11 
where L^''^ is inserted into the internal line of the gluon tadpole. The four-gluon 
vertex appearing there is denoted by T'^^^^{k, q, r, s). Its explicit form can be 
found in jl6j . for example; it is not important here. Then the contribution 
coming from this graph is proportional to: 



7|4(2vr)^^(p + . + r + «)r^45^(p, ,,,, 
P'^BZ BZ 



(3.67) 



where the Delta-function coming from the four-gluon vertex has been extraced 
explicitly from T. Splitting this four-dimensional Delta- function up into the 
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spatial and temporal components and looking now only at the relevant compo- 
nent where = = 0, this reduces to: 



BZ BZ 



p,\ 
1 

2 



Carrying out the four s-integrations, using the second and third Delta-function, 
yields: 

E / 7|^(2^)^(P + ^Ti^^^m 0), (g, 0),r, -r) ^'^^~J'"' ■ (3-68) 



P'^BZ 



Using (|3.1fl)) , one sees that this is the same result which one would have obtained 
if one would have inserted S^^^ into the internal line of the gluon tadpole diagram 
and then looked only at the contribution for ^4 = ^4 = 0. Hence the result is: 

hm < hm I < S^'^^J'^S^'^ >,onn, (3.69) 

r^oo r^oo T 

if only the insertion into the internal line is considered. But looking at the 
results for insertions into the external lines obtained above, one sees that this 
formula is true even if all possible insertions are considered. 

Exactly the same arguments can be made for the fourth graph in KS.lL where 
L^^^ is inserted into internal lines of the gluon loop. Here, too, one obtains 

hm < (s^Y >_= lim I < fgiDy (3 70) 

using the Delta- functions from the two three-gluon vertices. Using the results 
obtained above for insertions into external lines, the end result is: 



hm I - < LW^(2)^(2) + < ^(0)^(2) 1 ^gil)y >^ 



_ / r(0)^(2)o(2) _ / r(0)^(2)o(2) . ] 

<^ 1j lu Jpp ^conn ^ UJ >-'meas ^conn I 

= hm l(-< 5Wa;(2)5(2) + < S^'^J'^]- (s^'A' 

f ^00 T\ 2 V y 

— < S"^ ■'w^ Spp >conn — <i ■'w^ ■''S'^eas ^conn ^ i (3.71) 

where now insertions into all gluon lines are allowed. 

What remains are the two additional vacuum polarization graphs which were 
depicted in figure E31 incorporating L^^^ and L^'^\ First look at the second of 
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these graphs, which involves a gluon loop where one of the two three-gluon 
vertices is replaced by L^^^ . As already mentioned above, the insertion of L^^^ 
at this place amounts to using a slightly altered version of the three-gluon vertex 
T'^^^{p, k, q). The usual form can be found again in [TB], for example; it comes 
from S^^\ Doing the same derivation as outlined there, but for L(0) instead 
for S, one gets: 



igo{2TTf6{p + q + k)fABC 



^(p+,+/c)(A+A)/2^^^ cos{k^/2){p'^q)x\ . (3.72) 



Using this, the graph 




k, a, D 



gives a contribution proportional to 

BZ BZ 

E -«)rS''''^(fc, r, (3.73) 

where the Delta-functions have been extracted from the Fs and are displayed 
explicitly. Doing the integrations over the three spatial components of s using 
the first Delta-function, this gives: 

[ 4i-J $7 . (2vr)^(g4 - .4 - r,){27rf6{k + q) 

(2vr)*_J^ 27r f.2 ('(f _ ^2 + ^2j 

E r^.T(Q, -r, (f, -^4))rit^;;''''^(fe, r, (-f, .4)). 

Now in the limit T — > cxo, again only the contribution for (74 = /c4 = is needed. 
Then the first Delta-function can be used to do the S4-integration, and one gets: 

^ r d^T {27rf6{k + q) ^ABC((^Q^ _^ . iL),DBC.,r . . 

Looking now at (|3.72|) . one sees that in this special case the phase factors give 
only factors of one, so that 

rS''''''((^> 0), r, -r) = Ci^((fc, 0), r, -r). 
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and the contribution from the graph is simply 

5: / f ; '\ 0), -r, r)rZTiCk, 0),r, -r). (3.74) 

But that is identical to the contribution which the general gluon loop graph 
would give if one would also consider only 54 = A;4 = there! Hence the final 
result is: 

hm < j.^^ \ ^ (5(i))'c^(2) . (3.75) 

Exactly the same type of argumentation can be used for the first graph, with 
the insertion of , so that one gets also: 

lim < L(2)a;(2) >^„„„= Hm 4 < S(2)a;(2) >^„„„ . (3.76) 

Summarizing: in the limit T — > 00, all vacuum polarization graphs con- 
tributing to the correlator of L(0) and W give the same contribution as the 
vacuum polarization graphs contributing to the correlator of S and W ^ divided 
by f. 

The spider graphs 

There are two spider graphs here: the usual one with the three-gluon vertex, 
where L'^^^ is inserted into a gluon line (the eighth in figure IXT|) . and the one 
where the three-gluon vertex is replaced by L*-^^ (the third in figure l3^ . For 
the normal spider graph, one has: 

lim \ < 5(^)^(2) >o= 0, (3.77) 

f ^00 T 

hence in the limit of high T, the contributions of these two spider graphs should 
also vanish. 

Unfortunately there is no such elegant argument here as in the case of the 
vacuum polarization graphs; the contributions have to be calculated explicitly. 
The second one is simpler tban the first; one gets for it: 

^ uj ±j --^ conn 

C2{F)C2{G) f d^p f d'^k } dq^ sm{p3R/2)sm{p^f/2) 



BZ 



^ cos((fc + p)4/2) cos((p + fc)3/2)(fc - p)4 g,^^^ sin((p + fc)3fi/2) 

4 ; — — cos /ca it/ 2 — — —— 

sin(p3/2) \ i I > + ^)^/2) 

sm.{kif/2) 



sin(/c4/2) 
cos{pi/2){2k+p)z 



sin(p3/2) 



cos(<74(r - l)/2) 

sm.{pi/2) cos(g4/2) cos(A;4/2) 
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sin(g4f /2) sin(fc4f /2) , . a, ^ 

sin(g4/2) sm(K;4/2) 

_^ cos(p4/2)(2fc +p)3 sin(pg^/2) sin(p4/2) cos((fc + q)^/2)OT{q4, h) 



-4 



sm(p3/2) 

cos(p4/2) cos2(A;4/2)(2irf A;)3sin((p + A;)3i?/2) sin(/c4f/2) 



sm(p4/2) 
■ cos(A:3i?/2) cos(g4(T' - l)/2 
^ ^ cos(p3/2)(g - fc)4 

sin(p4/2) 
•Oii(A;3,-P3 -^3). 



sin((p + A;)3/2) sin(fc4/2) 
cos(p4/2) cos((g + /s)4/2) sm{{q + k)if/'2-) 



(3.78) 



In the limit T — > 00, terms like 



sin(g4r/2) 
sin(g'4/2) 

give constant contributions, whereas sines and cosines which depend on T, but 
are not divided by the corresponding sines of the coordinates, give contributions 
proportional to T~^ because of the fast oscillations of these functions. Therefore 
all of the terms in the expression above go to zero in the limit: 



lim < a;(^)L(i) >, 



0. 



(3.79) 



The first spider graph is much more complicated: 



l,5,e,ipQZ BZ -TT 

2 

g-i(p+fe+q)7/2 

q=—p—k 



p + k 



p^P {p + k^^ {(p + k\ + 



Sj^ (^p + k^ - + ^) (p + k 



-i(p+k+q)^/2 



\ 



p- 



'k'^ (^p + k^ i (p + k 



hs — ^45 



sm(j)3R/2) sin(p4f /2) 



sin{ps/2) 

sin(fe3^/2) sin(g4f /2) 



2S3eSi^ 



sin(fc3/2) sin(g4/2) 
-i sin((p + A;)3-R/2 + A;4f /2)) 

+S3eS3,p(^- isin((g - k)if/2) 



cos{{p + k)3R/2 + kif/2) 



sin(A;3i?/2) sin((j; + A:)3i?/2) 
sin(A;3/2) sin((p + A;)3/2) 
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+2sin((g + k)4f/2)OR{k3, -P3 - ks)) 

sm{k4f/2)sm{(p + k)sR/2) ^ . 

^ sin(A;4/2) sin((p + A;)3/2) v ' ' 

+zsin(g4r/2 + A;3^/2)) 

, , (■ ■ n o,N A /ox 8111(^4^/2) sin(g4r/2) 
-(54e54^Usin p + 2fc 3i^/2 . \r /-/ ■ 7 ' 

sin(A;4/2) sin(g4/2) 

-2sin(p3i?/2)OT(fe4,(?4))j 

- (-TT^TW^ + — '^^((^^ + sm(g4f /2) 

Vsin((p + /c)3/2) 3111(^4/2)/ 

/ sin(fc3^/2) sin(p4T/2) ( . ^ 

• 2(53ed45^-7— — — — cos(p3i?/2 - A;4T/2) 

y sin(fc3/2) sm(p4/2) V ' ' 

+ism(p3^/2 + fe4f'/2 



^ / . .X ^ sm(fc3i?/2) sin(p3/2/2) 

+53e<^35( - ^sin p - ^4^/2 .\l > .7 ;/ 
V sm(A:3/2) sm(p3/2) 

+2sin((p + fc)4f /2)Or(A;3,P3)) 

sin(fc4T'/2) sin(p3.R/2) ^ . . 

sm(A;4/2) sm(p3/2) ^ 

+ism{p4f/2 + k3R/2 ' 



-d4eS4S[ism{{p - k)3R 2) . . . 

V sm(A;4/2) sin(p4/2} 

-2sin((p + /c)3i?/2)OT(/e4,P4) 



(3.80) 



Here the explicit expression for the three-gluon vertex Tg^^ was not inserted 
and the sums were not carried out; this would only have given an even more 
complicated result. This is not necessary because even here one can already see 
that all of the terms which appear vanish in the limit (using the same arguments 
as above). Hence the result for this spider graph is also: 

lim < L^^^J^^S^^^ >conn= 0. (3.81) 
T-^-oo 

Thus indeed both spider graphs vanish in the limit of large T. 



The graphs with two independent gluon lines 

As usual, a;W is splitted up into its parts and every contribution is calculated 
separately. First one obtains, again using the symmetry properties of the colour 
indices: 

< L(0)a;("^) >conn= 0. (3.82) 



56 



CHAPTER 3. ACTION SUM RULE 



The contribution oi u^^'^ which is proportional to (C2(-F))^ only comes from 
u^^^^ and can also be easily calculated: 



< 



--< LO 



(2) X 



+C2iG)C2{F) 



(3.83) 



In the hmit of large f, < L^^^u^'^^ > 

conn gives the same contribution as < 
5(0)^(2) >^j,„„ /T and therefore < L^'^^uj^'^^^ >conn gives the same contribution 
as < 5Ww(4A) 

The other parts are much more complicated. Because of the non-conserva- 
tion of the momentum when L^^^ is inserted, (|3.23p no longer applies, so that 
the insertion into both lines has to be considered now. This leads to a doubling 
of the number of the graphs which have to be computed. 

It is convenient to split the various contributions up again, classifying them 
according to the number of links in spatial and temporal direction. First look 
at the graphs where only spatial links appear; denote the relevant parts of a;^^^ 

by ^RR- 



< ^ ^RR >c 



(27r)4 J (27r)4 J 2tt p 



-J 



C2{G)C2{F) f d^p f d^k f dq^ I f + \pl + \qi 



BZ BZ 

sm\p3R/2) sin^(fc3j?/2) . 2 - . ^ . 

— ^Wn^^^"^ {k4T/2)sm{piT/2)sm{qiT/2) 

sm^(p3/2) sm^(A;3/2) 

^ r(o), ,(4iJ) 

ij ^RR -^conn 



C2{G)C2{F) f d^p 



BZ 



d^k 



(27r)4 J (27r) 

BZ 



Pi 



27r A;2 



p'^ (yP + qj 



sm{{k - p)if/2) sin((fc + q)Af /2) 



sm^{p3R/2) sin^iksR/i) 



sm\p3/2) sin2(A:3/2) 
4sin((p + k)^f/2) sm{{q - k)4f/2)OUps, k^) 



^ ^ ^RR 



C2{F)C2{G) f d^p f d'^k 



(27r)4 J (27r) 



dq4 1 + \pl + \ql - Pl 



■s\ii^{kiT/2) 



BZ BZ 

sm{ksR/2) 



2tt 



p2 (^p + qfj 



S2)L„_i.„ + ■ ^r{P3,-P3) + —T-, -—T.R{p3,k3] 



sin(A:3/2) 
2 (El - S2)|p3„fc3 

^ RR. -^conn 

C2{F)C2{G) 



cos(p4T/2) cos(g4T/2) 



(3.84) 



(3.85) 



(3.86) 



sin(fe3/2) 



d'^k 



sin(p3/2) 



BZ 



d'p 

(27r)4 J (27r) 
BZ 



u^p +\pI + hoi 



Pl 



27r 



+ ql 
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cos(p4T/2) cos((/4T/2) 



+2 



sin(A;3/2) 

2 (El 

sin(/c3i?/2) sm2(p3i?/2) sin(p3i?/2) 



S2)|p3^fc3 + 2 Sr(0, -fcs) - 2 S,j(p3 - ks, -P3) 



sin(/c3/2) sin2(p3/2) 
Sr(p3, -Pa) 



sin(p3/2) 



sin(fc3/2 
C2(F)C2(G) 



BZ 



(27r)4 y (27r)4 7 27r fcS 



BZ 



• cos(p4r/2) cos{q4T/2) sin^(A;4r/2) 



sm(/c3i?/2) 



sm(fe3/2) 

sin(A;3i?/2) sin(p3i?/2) sin((p + k)3R/2)' 
' sin(fc3/2) sin(p3/2) sin((p + A;)3/2) 



sm(A;3i?/2) 
sm(A;3/2) 



(3.87) 



(3.88) 



Because of the fast oscillations of the two cosines for T — > oo, the contributions 
of these integrals vanish in the limit of large T. 

The next group consists of the graphs with two spatial and two temporal 
links, where L^^^ is inserted into the line connecting the two spatial links. The 
corresponding parts of a;^^) are denoted by oo^^tv 



RTl -^conn 



C2(G)C2(F) f d'^p 



(27r) 

BZ BZ 



d^k 
(2^ 



dq4 1 P + \pI + \ql - p\ 
27r A;2 



sin2(j33i?/2) sin^(A:4r/2) 
sin2(p3/2) sin2(A;4/2) 



C2(G)C2(F) [ d^p 



sin^(/c3i?/2) sin(j94r/2) sin(Q4r/2) (3.89) 



i (2^7 (2^7 



dq^\_p +^pl + ^ql-pi 
27r jfc2' 



(2vr)^ 



BZ BZ 

sin2(p3^/2) sin2(A;4f/2) , ^ , - 

• ■ 27 /oV ■ 2tl J cos(y4r/2) cos(g4r/2) 
sm (^3/2) sm (K4/2) 

RTl ^conn 



P 



(3.90) 



C2(G)C2(F) f d^p f d^k 



BZ 



(27r)4 J (27r)4 J 27r A;2 



dq4 1 P + \pI + \ql - pI 



BZ 



(p+qI) 



sin^(fc4r/2) 
sin2(fc4/2) 



sin^(/c3^/2) cos(p4r/2) cos(g4r/2)Sij(p3, -^3) 
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sin2(A;4r/2) 
sin2(fe4/2) 



sm\k3R/2) cos(j)iT/2) cos{qiT/2) — 



sin2(p3i?/2) 



sin^(p3/2) 



(3.91) 



/ 7-(0), 



dqA I f + y>\ + jqj -Pj 



6 



sin2(fc4r/2) 



+2 



sin2(fc4/2) 
sin2(A;4f/2) . 2 



sm^{k3R/2) 008(^4^/2) cos (94^/2) Ei^(p3, -pa) 



sin2(A;4/2) 



sin^(fc3i?/2) cos(p4r/2) cos{qiT/2) 



sm^{p3R/2) 
sin2(p3/2) 



(3.92) 



^ 7-(0), 

C2(G)C2(F) f d^p 



d^k 



BZ 



(27r)4 J (27r) i i 27r A;2 



c/g4 1 ^ + ^P4 + I'?4 - pI 



BZ 



{p + ql) 



sin2(/c4r/2) . 2 



sin2(A;4/2) 



sm^{k3R/2) cos{piT/2) cos{qiT/2)R. 



(3.93) 



For T — > 00, the factor '^"Z ^^^^2')^ gives a linear dependence on T, but the two 

cosines both give factors of T^^, so that in total all of these integrals go with 
T~^ in the limit and hence vanish. 

Then there are the graphs with two spatial and two temporal links where 
is inserted into the line connecting the two temporal links. The corre- 
sponding parts of a;^^^ are denoted by i^^^2'- 



< 



RT2 ^ conn 



MG)C2iF) f d^p 



{2ny 



d'^k 



dqi 1 



P 



27r A:2 ^2 



{p +Qa) 



, , (27r)4 

BZ BZ 



< 



sin^(fe3/2) 



tan(p4/2) tan(g4/2) 



UJ 



,{4B) 



RT2 > 

conn 



C2{G)C2{F) f d^p f d^k f dq^ 1 



BZ 



(27r)4 J (27r 



p 



BZ 



27T fe2^2(^ + ^2^ 



sin2(A;3i?/2) sin(p4f /2) sin(g4f /2) 
' sin2(/c3/2) tan(p4/2) tan(g4/2) 

^RX2 ^ conn 

C2{G)C2{F) f d'^p f d^k J dq4 1 



r d^p r d^k r 
J (2^ J (2^y 



p 



27r 



BZ BZ 
o,,^ ^,^,sin2(A;3i2/2) 



k'p'{if + qi) 



3sin^(/c4r/2) — 



sin^(A;3/2) 



cos((p + g)4/2)Sr(p4, 54) 



(3.94) 



(3.95) 
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sm2(/c3i?/2) 



sin(j;4f /2) sm{kif/2) 



-V(A;3/2) 

sm(j94f /2) sm(g4T/2) 



tan(p4/2) tan(q'4/2) 



+ sm^{p3R/2)ER{k3,-k3) 
C2{G)C2{F) 



6 



tan(p4/2) tan(gf4/2) 

TT 

(ig4 1 



(3.96) 



(Tk 



fs - 008(^4^)1 W{p3R/2) — 
V / si 



+ (l + 2cos(p3i?)) sm2(/c4f/2) 



sm2(A;3i?/2) sin(p4r/2) sin(g4r/2) 
siT?{k3/2) tan(p4/2) tan(g4/2) 

sm2(fc3^/2) sm(p4f /2) sm(g4f/2) 
sm^{k3/2) tan(p4/2) tan(g4/2) 



-3sin2(A;4r/2) 



sm 



,R/2) 



sm^(fc3/2) 



cos((p + q)i/2)T,T{pi, qi) 



-sin2(p3i?/2)Sij(fc3,-A;3) 



C2{G)C2{F) f d^p 



6 



sin(j34f/2) sin(g4f/2) 
tan(p4/2) tan(g4/2) 

7r 

(ig4 1 



(3.97) 



(27r)4 

BZ BZ 
,2 



(27r)4 7 27r fc2 .3 



3sin2(fc4T/2)^ 
+^sin2(p3^) 



sm\k3R/2) sm{{p + q)iT/2) 



sin2(A;3/2) tan((p + 9)4/2) 
sin(j)4T/2) s'm{q^T /2) 



tan(p4/2) tan(g4/2) 
Not all of these terms vanish seperately, but their sum gives: 



(3.98) 



C2{G)C2{F) f d'p 



d^k 



BZ 



{2^f£ (27r)4 J 27r k^ p2 ('^ 



dq4 1 



BZ 



sm^{k3R/2) sin(p4f /2) sin(g4f /2) 
' sin2(fc3/2) tan(p4/2) tan(g4/2) 



cos(A;4r) cos(p3A). 

(3.99) 



This expression vanishes in the limit of large T, again because of the fast oscil- 
lations of the cosine. 

The next large group consists of the graphs with three spatial links and one 
temporal link. Denoting the corresponding part with uj^rhjit, they give: 

^r(0),,(4A) 

^ ^ ^RRRT ^conn 



C2{G)C2{F) f d^p 



BZ 



d^k 



(27r)4 J (27r) 

BZ 



dq4 1 Pl + ql 
'^^f(p+qt) 
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o7 r^sm [kiT 2) — 7— ; — sm ip^R/z) (3.100) 

sm2(fc3/2) ^ 4 / y tan(g4/2) ^^■^ ^ ^ ^ ^ 

^ ^RRRT -^conn 

C2{G)C2{F) r d^p r d^k J dg^ l ^3^4 sin(g4T/2) 
2 (27r)4 (27r)4 27r p ^ + ^2) tan(g4/2) 



BZ 

sin- 
sin 



L^(p3^/2) sin(fe3^/2) . 

w — 7T\ • /o^ - 2cos(^33i2/2)OR(p3,g3) 

in(p3/2) sin(A;3/2) j 

\kifl2) cos(p4T/2) sin(p4/2)) 

(3.101) 



sin(A;3^/2) 

■ sin(A;3/2) ^ ^ , 

■ (cos^(A;4f /2) sin(p4f /2) cos(p4/2) 



— sin 



^ ^RRRT ^conn 

C2{G)C2{F) r d^p r d'^k J dg^ l sin(g4r/2) 
2 y (2^7 (2^7 2^fc2 .2/^ , tan(g4/2) 

■ (sin(i?4(r + l)/2) + sin(p4(T - l)/2)) sin(p3/2) sin(p3^/2) {^ir - ^2r) 

+ sin2(p3i?/2)sin(p4(T+l)/2)S^(fc3,-A:3) 

^ Tin{k3/2) '^^^3i?/2) sin(p3/2) sin(p4(r - l)/2) cos(/c4r)ER(p3, ^3) 
+ sin(p4r/2) sm^{p3R/2) cos(p4/2)Sr(A;3, -A;3) 

- sin(p4(r - i)/2) ^^'^(y/Jj) sin(p3/2) sin(p3i?/2)SR(fc3,^^3 
sin(fc3//j 



- sin(p4(T' - 1)/2)^^|^ sin(p3/2) cos{psR/2)ORiks,ps) 
+ sin^{p3R/2) (sm'^{k4f/2) sm{p4{f - l)/2) + sin(p4f'/2) cos{k^f)^ 

+2sin(p4(r - l)/2) sin^(fc4f/2) ''"^y/f cos(p3-R/2) sin(;,3/2)Oii(A;3,P3) 

sin(/c3/2) 

-2cos(p4r/2) sin(p4/2) sin^(fc4f/2) ^in\p3R/2) 

sm^(A:3/2) 

^ ^RRRT ^conn 

C2{G)C2{F) r d^p r d^k J dg^ 1 sin(g4T/2) 

12 J (27r)4 7 (27r)4 J 27r fc2 ^2 , ^2] tan(g4/2) 



(3.102) 



(sin(p4(T + l)/2) + sin 
• (El/? - T.2R + Si?(0, fcs) - Si?(p3 + A;3, -^3)) 

"'"'^f!^ sin'ipsm (sin(p4(f - l)/2) + 

sin^(A;3/2) ^ 



+ 



.(P4(T - l)/2)) sin(p3/2)sin(p3-R/2) 

.(p4(r + l)/2)cos(A;4f)) 
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- sin(p3/2) sm{p3R/2) sm{p4{f - l)/2) cos(fc4f) ''"^ff ^{^^ ^^(fcs, j^g) 

- sin2(p3i?/2) sin(p4(T + 1)/2)Ek(A:3, -^3) 

+2sm\p3R/2 f'!''i^fl^^ (sin(p4r/2) cos(p4/2) 

sm^(/C3/2) V 

- sin2(A;4T/2) sin(p4(f + l)/2)) 

+ sin2(p3^/2)^^^^%^(sin(p4f/2)cos(p4/2) 

sm (K3/2) ^ 

- sin2(/c4T/2) sin(p4(f - l)/2)) 

+ sin2(fc4f /2)^^^%^ sin(p4(T - l)/2) sin(p3/2) 

sm(A:3/2) 

• (sin(p3^/2)SR(p3, A;3) + cos(p3-R/2)Or(p3, ^3 

- sm\p3R/2) sm(p4T/2) cos(p4/2)Sfl(/c3, -/e3) 

+ ■ /ON ^™ ^^3-^/2 

sm (A;3/2j 

■ (sm^{k4f/2) sin(p4(f - l)/2) + cos(A;4T' ) sin(^?4f /2) cos(p4/2)) 

^r(0),,(4f?) 
^ -'^ ^RRRT ^conn 

C2{G)C2{F) r d^p r d'^k J dg^ l p4 sm{qif/2) 

12 J (2^7 (2^7 + a^) tan((/4/2) 

?{kzR/2) sin(p4f /2) cos (594/2)^ 

■(fc3i^/2) sin((p + k%R/2) . ^ . 

— 77 7^ — sin(r?3it/2) sin(»3/2) 

..11(^3/2) sm((p+fe)3/2) ^ ^ ^ 

sin2(fc4f /2) + cos(/c4T)) sm{p4{f - l)/2) 
+ (sm(p4(r + l)/2) + sin(p4(r + l)/2)) sm{p3R/2) sm{ps/2) 

■ (Sii(p3,0) - T.R{k3,P3 - ^3)) 

+ sin(p4(r + l)/2) sm^{p3R/2)R . 



(3.103) 



sm' 

sin 
+ . 

sm 



(3.104) 



All of these contributions go to zero in the limit of large T. 

What remains are the graphs with four temporal links and the graphs with 
three temporal and one spatial link. Only these can give non-vanishing contri- 
butions in the limit of large T; the first ones will be denoted by uj!^j, and the 
others with w^j^j^y "^^^ explicit results are: 
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= -2 



C2(G)C2(F) 



P 



2^ 1 (27r)4 J 2^fc^ 

BZ ' BZ -TT 



sin2(A;4/2) 

C2{G)C2{F) 



■V -u ujrpj, ^conn 



tan(p4/2) tan(g4/2) 



d^/c f dq4 1 



sin^(p3i?/2) (3.105) 



BZ 



(27r)4 y (27r) 



P 



BZ 



27r fc2^2 (^ + g2 



+4.sm{{k - p)3R/2) sm{{k + p)3R/2) 



sin(j34f /2) sm(^4f /2) sm'^{k4f/2) 
tan(p4/2) tan(g4/2) 8111^(^4/2) 

sin(g4f/2) sin2(/c4r/2) 



tan(g4/2) siii2(A;4/2) 



-4sin2((fc +p)3i?/2) cos((p + g)4/2)OT(P4, A;4)OT(g4, -^^4) 



Ot (^4,^4) 



(3.106) 



^ j-v LUj-irp -^conn 



C2{G)C2{F) f d^p f d^k f dqi 1 



I (2^0^ y (2^7 



+2 



-2 



tan(p4/2) 

sin(fe4f/2) 
sin(A;4/2) 
sin(p4f/2) sin(A;4f /2) 



(27r)4 {2-kY J_ 27rfc2^2(p-^ + ^2^ 
sin^(i>3i?/2) (Sit - S2r)L,=_,, cos(g4/2) 



^sm(p4T/2) ^. 2 



sm\k3R/2) (Est - S4T) cos((p + 9)4/2) 



Or(54, A;4)sin(g4/2) 



tan(p4/2) sin(A;4/2) 

(sin^(A;3^/2) cos{p3R) - s\n^{p3R/2) cos(A;3^)) 

sin(p4f/2) sin(g4f /2) a, . 

/ ; vy4 / ^sin2(p3i^/2)Sr((?4,-g4) 



tan(p4/2) tan(g4/2) 
sin2(A;4T'/2) 



sin^(/c4/2) 



■sm^(A;3i?/2) 



(St(p4, 94) cos((p + q)i/2) - Ot{p4, qi) sin((p + g)4/2)) 



sm(p4r/2) sin(/c4T/2) 
tan(j»4/2) sm(A;4/2) 
(ST(g4, ^4) cos(g4/2) - Or(g4, A;4) sin(g4/2)) 
sin(p4f /2) sin(A;4f /2) 
tan(p4/2) sin(fc4/2) 

S7'(fc4, (?.i) cos(g4/2) - OrikA, q^) sin(g4/2)) 



+ sm^{p3R/2) cosikaf) 
+ sm^{k3R/2) cosipsf) 



(3.107) 
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C2{G)C2{F) r I d^k J dqjl f 



+ 



+2 



(Sit - S2T)|p4=g4 

^;i^sin^(*3A/2)cos((p + ,)./2) 



sin(fc4/2) ™ vo^-/ -/ vvi- ' -i/"*/ -/ 

- T,4t) + St(P4 + q4, -k4) - St(p4 - k4, Qa)] 
Af/2) sin^rfc.f /2) 



■ [(^3T ~ ^41' J -r ^ryi^A t y4) —r^AJ — 
sin(p4f /2) sin(q'4f /2) sm^(k4f/2) 

tan(p4/2) tan(g4/2) s{n^{k4/2) 

■ {sm^{pzR/2) cos^{k^R/2) + sir? [k^R/ 2) cos^(p3^/2)) 

„sin(p4f/2)sin(A;4f/2)^ , . , 

tan(p4/2) sm(/c4/2) 

■ (s\n^{pzR/2) cos{k^R) - s\n^{kzR/2) cos{p^R)) 

1, A, X A^ sin(p4T/2) sin(fc4f /2) 

tan(p4/z) sm(K4/2) 

• i^T{k4-, Qa) cos(g4/2) - OrikA, Qa) sm(g4/2)) 

. 9/, A ,„x sin(A;4r/2) , sm(p4r/2) 
- sin^ fc3i^/2 . \^ ; / cos P3i^ , ^7 ' 
sin(fc4/2) tan(p4/2) 

• (St(A;4, 94) cos(g4/2) + OrikA, Qa) sin(g4/2)) 



2., A, xSm^(A;4f/2) 
■ (St(p4, 94) cos((p + g)4/2) - OriPA, Qa) sm{{p + q)A/2)) 



^ -L/ ujrprp -^conn 



C2{G)C2{F) r d^p r d^k J dqAi p 

W{p2R/2f^^^l^ (St(P4,0) - ^T{kA,PA - kA)) 



-{kAf/2) 



+ sm'{k3R/2 f''!^^fl^^ {T.T{kA,PA + Qa) - ^Hpa^Qa + kA)) 
sm(k4/2) 

s m(p4f /2)sin(g4f /2) ^ ^^^ ^ 



+ 



tan(p4/2) tan(g4/2) 



sin^(fe 



,.4T/2) sin((p + g)4f /2) . ^ . 
ikAl2) tan((p + 9)4/2) ^^^^/'^ 



(3.108) 
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sin(p4r/2) sin(fc4T'/2) sin((fc + q)if/2) 
tan(p4/2) sin(A;4/2) sin((A; + q)4/2) ^ 

• (sin^(p3^/2) cos{k3R) + sin^(A;3^/2) cosipsR 



(3.109) 



RTTT ^conn 

C2{G)C2{F) r d^p r d^k J dg^ 1 pI + qI 

sin^(fe4T/2) ,,^.(,^^^,) Sm(,4T/2) sm(,4T/2) ^.^ 
sin2(A;4/2) ^ * / ^ tan(p4/2) tan(g4/2) ^/-'^ / v j 

RTTT ^conn 

-C2{G)C2{F) J^J^J 

BZ BZ -TT y [P 

• (cos(p4f /2) cos(p4/2) cos(g4/2) - sin(p4(r - l)/2) sin(p4/2)) 

(3.111) 

RTTT '^conn 

C2{G)C2{F) r d'^p r d^k 1 dq^}_ l_ 

2 J (27r)4 y (27r)4 7 27r fc2 ^2 f ^ , 



+ «4 



2p4 cos(p4/2) sin(p4r/2) cos{qi/2) sm^{p^R/2) (Sit - ^2T)\p^=q, 
-2p4 cos(p4/2) sin(p4f/2) sin^(p3i?/2) cos(fc3^) 3111(94/2)0^(94, ^4) 
+P4 cos(p4/2) sin(p4T/2) sm^{p3R/2) ""'^^^^J^^ ^T{k4, -k^) 



-P4 cos(p4/2) sm(p4T/2) sm'{psR/2) cos(A;3i?)^^^%^ 

sin(K4/2) 

• (St(94, k4) cos(g4/2) - OT(g4, A;4) sm{q4/2)) 
+2p4 sin(p4(r + l)/2) sm''{psR/2f-^^^I^^T{k4, -k4) 

-2sin^(fe3^/2) ''^^f;7o? sm^(y3^/2)p4COs(p4/2)sin(p4r/2)ST(fc,g4) 

+ sm2(p3i?/2)^^^^J;^^4 sin(p4(r - l)/2)ET(fe4, -k4) 

, ,sm(A;4f/2) 
+ sin^ fc3i?/2 

sm(fc4/2) 
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• (g4 (sin(p4(T' - l)/2) sin(p4/2) - cos(p4r/2)) Orih, q^) 

-P4 sin(p4(T - l)/2) cos(g4/2)i;T(/c4, ^4)) 
+P4 (sin(p4(T - l)/2) - sm(p4(T + l)/2)) sin2(p3i?/2) 

81112(^4^/2) Sm(g4f/2) / ^ .2., 

9/ ; / — ; X (cos(A;3i? + sin2(A;3i2/2)) 

sin2(A:4/2) sin(g4/2) V ^ V ^ / 

+2sin2(A;3^/2)^%^ sin2(^,3^/2) 
sin(K4/2j 

•54 (sin(;j4(r - 1) sin(;j4/2) - cos(p4r/2)) Orih, qi) 

RTTT '^conn 

C2{G)C2{F) r (fp r d^k } dq^ l 1 
3 7 (27r)4 y (27r)4 7 27r ^2 .3 f ^ , 



2p4 cos(p4/2) sin(p4f/2) cos(g4/2) sm^{p^R/2) 

■ (Sir - S2T)|p4=g4 + St(0, 54) - St(A;4 + 94, -/^4) 

. . . 9/ A, , sm2(A;4T/2) sin((/4T/2) 

+2P4COS P4/2 sin^4r/2 siv?{jp^R 2) . j f / ^ ^'^^ > 

sm (fc4/2) tan(g4/2) 



■003^(^3^/2) 



(3.112) 



+2p4 cos(p4/2) sin(p4T/2) sui^{p^R/2f^^$^^^ cos{k^R) 
•sin(g4/2)Or(A;4,g4) 

-pi cos(p4/2) sm(p4f /2) sin2(p3^/2)^^^^^^^^ST(/i:4, -/i:4) 

^ ^ ^ sinf/c4T/2') 

-P4 cos(p4/2) sin(p4r/2) siT?{pzR/2) cos{k^R) . \, 

Sin(K4/Zj 

• (i;r(A:4, 94) cos(q'4/2) - Ot(^4, 94) sin(g4/2)) 

-2sin2(p3i2/2)sm2(fc3i?/2)^%^g4 

sin(A;4/2) 

• (cos(A;4f /2) - sin(p4(f - l)/2) sin(p4/2)) Or(fc4, ^4) 

+fo sinfe(f - l)/2) sm^(p3iJ/2) cos' feii/2) -'°'^^f /f f {f ' 

Sin (A;4/2) tan(g4/2j 



,V p4 cos P4/2 sm 094^/2 
sin A;4/2 



+2 sin2 (p3i2/2) sin^ ( A;3i?/2) 

sin(A;4/2) 
•cos(g4/2)ST(g4, A;4) 

-2sm\psR/2)p4 cos(p4/2) sin(p4T'/2)^^^^^i^ST(A;4, -^4) 
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+ sm\p3R/2) sm\q3R/2) 



sm{kif/2) 



sin(A;4/2) 

•94 (cos(p4r/2) - sin(p4(T - l)/2) sin(p4/2)) OT{q4, ^4) 
+ sm\psR/2) sm\q^R/2 f'!^''fl% ^ sm{pi{f - 1)/2)St(<?4, h) 

- sin2(p3i?/2)p4 cos(p4/2) sm{p^f/2f-^^I^^T{k4, -h) 

sm^{k4f/2) sm{qif/2) 
sin2(/c4/2) tan(g4/2) 



+Pi (sin(p4(T + l)/2) + sin(p4(r - l)/2) 



•sm2(p3i?/2) (sin2(A;3i?/2) +cos(A;3i? 



/ r(0), ,(4£:) 

RTTT ^conn 

C2{G)C2{F) f d^p 



BZ 



(27r)4 J (27r) 



(ig4 1 



27r fc2^2(^^ + ^2 



2sin(p4(r + l)/2) sin2(p3i?/2) 



sm{qiT /2) , 



tan((74/2) 

-2 sin(p4(f + l)/2) sm2(p3^/2) sm2(/i:3^/2) 
sin(A:4T/2) sin((g + k)if /2) 



sin(/c4/2) sm((g + A:)4/2) 
+ sin(p4(T - l)/2) sin2(p3i?/2) 



■ cos(g4/2) 
sin(g4T/2) , 



tan((74/2) 

- sin(p4(T - l)/2) sin2(p3^/2) sm2(A;3^/2) 
sin(fc4f /2) sin((g + fc)4T/2) 

■ sin(/fc4/2) sm((g + A:)4/2) ^ 
+2sin2(p3i?/2) sin(p4f /2) cos(p4/2) cos(g4/2) 

• {Tsriqi, 0) - St(A:4, Q4 - ^4)) 



+ sin2(p3^/2) sm(p4f /2) cos(p4/2; 



sm(g4T/2) - 



tan((74/2) 

+ sin^(p3i?/2) sm(p4T/2) cos(p4/2) cos{p^R) 
sm{kAf/2) sin((g + k)4f/2) 



sin(A:4/2) sm((g + A:)4/2) 



■ cos(g4/2) 



(3.113) 



(3.114) 



The functions T,it and S2T are obtained from Si and S2 by replacing R 
with T, with p4 and k^ by A;4 (see section r2.2.2|) . The new functions S3T 
and T,4T are the even parts of the following sums: 



T-3 T-2 T-1 
Y,^rp = gj(fc4-P4-z)(t-l)/2 ^ ^ ^ ^-ikih ^ipih ^iqih 

h=0 /2=/l+l«3=/2+l 
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S4T 



f-3 f-2 f-1 



f-3 f-2 f-1 

=4-P4-2)(r-i)/2 ^ ^ J2 t 

il=0/2=/l + li3=/2 + l 



3(^4 



The explicit results are: 

1 sin(g4f /2) sin(p4(f - 2)/2) sin(A:4(T - 2)/2) cos 

2^qnn = 

" -"-'54/2) sin(p4/2) sm{k4/2) 

i{ip-k)4f-2)/2) 
sin((p 



4 sin((;4/2) sin(^ 
1 sin(g4T'/2) sin(p4/2 — ^4) sin( 
4 sin(g4/2) sin(p4/2) 



+ 4 



1 cos((p + g)4 , 

sin(g4/2) sin((p + 9)4/2) 
1 cos(p4/2 — ki 



Diii|^^4/ oiiHV/^ ~ ^)4/2) 

(f - 2)/2 + 94/2) sin(/c4(f - 2)/2) cos( 
Hm((n-\- a)J2) sin(A:4/2) 
_ sin(p + q- k)i{f - 2)/2) 

4 sin(g4/2) sin((p + 9)4/2) sin((p + q- h)i/2)) 
.f/2) sin(A;4(f - 2)/2) sin(p4(f - 2)/2) cos(p4' 



S4r - - 



sin(/c4/2) 



sm 



l(P4/2) 



1 sin(g4 
4 sin(g4/2^ 

1 sin((j4f /2) sin(p4 - ^4/2) sin((p - k)i{f - 2)/2) 
"'" sin(A:4/2) sin((p — A;)4/2^ 



4 sin(g4/2^ 

i((g -fc)4(t-2)/2 + g4/2) sin 



(3.115) 



1 COS^^^ --y^v- -,i - . 1^ 

sin(g4/2) sin((g — fe)4/2 



.(p4(T-2)/2)cos(p4) 
sin(p4/2) 



1 cos(p4 - fc4/2) sin(p + q - k)^{f - 2)/2) 
4sin(g4/2)sin((g- A:)4/2) sin((p + g - A:)4/2)) ^ ' ' 



These results with L^^^ above have to be compared to the expressions one gets 
by inserting S^^^ into the graphs with oj^*^) in the limit of large T. The graphs 
contributing to < S^^^ oj'^'^'^ >o can be split in the same way as outlined above, 
and with similar arguments as before, it can be shown that the contributions 
coming from u;^]j, (jJ^^, and t^^/j/j^T vanish in the limit of large T. Hence only 
the results for the graphs with uj^)p and oj^^rpj,, which can be obtained from 
the sections 12.2.21 and 1,3. 2. 2[ are given here: 



d2 



-2- 



C2{G)C2{F) f d^p f d^k 1 k' ^ + k sin^(j34T/2) sin^(fc4r/2) 

^2^^ sin2(p4/2) sin2(A;4/2) 



(27r)4 J (27r)4 p2 
BZ BZ 



sm^{p3R/2)sin^{k3R/2) 
C2{G)C2{F) f d^p 



(3.117) 



d'^k 1 }p + k 



BZ 



(27r)4 J (27r)4 p^ 



BZ 



sin (P4T/2) sin2(A:4T/2) 
sin2(p4/2) sin2(A:4/2) 



+40|,(p4,A;4) 



sin2(p3i?/2) cos2(/c3^/2) + sin2(A;3i?/2) cos^(p3i?/2) 



(3.118) 
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and 



(Pk 1 k^ + k 



BZ 



(2^)%4 i^^rp' (k^ 



A sm(p4r/2) 
. sin(p4/2) 



+ 



sin 

sin2(p4/2) I sin2(A:4/2) 



T + 



sin(p4r/2) sin(A;4T/2) 
sin(p4/2) sin(A;4/2) 



cos(A;3i?) 



sin(j34r/2) sin(A;4r/2) _ sin((p + /c)3r/2) \ 
sm(p4/2) sin(A;4/2) sm((p + A;)3/2) / 



12 



(2vr) 



^2 

d^/c 1 A;2 + 



•sm2(p3i?/2) (3.119) 



,2/^ ^-'/o^ c;^2/ 



(27r)4 p2 



(k^y 



+2 



•^ sin^(p4T/2) ^ 
sin2(p4/2) 

sin(p4T/2) 



sin^(p4r/2) sin^(fc4r/2) ^ A/o\ma 

sm2(p4/2) sin2(A:4/2) ''''^ ^"^"^^ ^> ^ ^ sm{p,/2) 

sin(p^r/2) sin(A:^r/2) sm{{p^ + fe^)r/2) 



(2 - cos{k3R)) 



/ q(o), ,(4£;) 



sin(p,,/2) sm(fc,,/2) sin((pi, + ku)/2) 



(Sir — S2t) 

sin2(p3i?/2) 



(3.120) 



C2(G)C2(F) /• d*p f d'^k le + k 



2,2 r- 



BZ 



(27r)4 J (27r)4 p2 



BZ 



fc2 



sin2(p4r/2) 
sin2(p4/2) 



rsin2(p3i?/2) 



- sm2(p3^/2) sm2(A;3i?/2) 



sin(p4r/2) sin(A;4T/2) sin((p + fe)4r/2) 



sin(p4/2) sm(A;4/2) sin((p + A;)4/2) 



(3.121) 



RTTT ^conn 



<q(0)^(4B) . 

^ RTTT ^conn 



= -C2{G)C2{F) 



d'^p 



(3.122) 

d^A; 1 4A:4 siiP{k3R/2) sinip^f /2) 



BZ 



(27r)4^7^ (27r)4p2 ^^2^2 sin(fe3/2) sin(p4/2) 



■cos2(p3i?/2) 



sin2(/c4f /2) sin(j04f /2) 
sin(A;4/2) sin(p4/2) 



+ 2cos(A;4r/2)Or(fc4,J34) 



(3.123) 



RTTT -^conn 

C2{G)C2{F) f d*p 



d'^k 1 k^kii 



BZ 



{2.YJ^ (2vr)4p2 



2sin2(p3i?) 



sin(p4T/2) sin2(fc3i?/2) 
sin(p4/2) sin(/j3/2) 



sin(A;4i?/2)Sr(A:4,P4) 
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+2si..^tefi/2) "'''f ""'f'"/'^ cos(faf/2)0.fe. k.) 

sin(p4/2) sm(/c3/2) 

+ sm {p3R 2)^- — — . 

sin(p4/2) sm(A;3/2) 

• (sm{k4f/2)T.T{p4, k4) + cosik^f /2)Ot{pa, k^) 



(3.124) 



^ ^ RTTT -^conn 



C2{G)C2{F) r d^p r d^k l kk 

sin'fef/2) °'°fe^/;' ^'"^fef si„(faT/2)E.fe. h) 

^ sin^(p3/2) sm(A;3/2) 

^9«^n2r. 5/o^ ^i^'(:P4T/2) sin^(/c3i?/2) sm^{k4f/2) 
+Zsm [p,K/Z) ^.^2(^^/2) sin(A;3/2) sin(A;4/2) 

, . o, ,sm2(p4f/2) sin2(/i:3i?/2) ^, , , 

-4sm\psR 2) ■ 2, L ■ V ^"^^4^/2 QTP4,fc4 
sm^(p4/2) sm(A;3/2) 

+2sin2r„ p/g^ sm^(l^4T/2) sin2(fc3i?/2) 

+2sin (psR 2) _ — — ■ , s 

sin^(p4/2) sm(A;3/2) 

■ (^sm{k4T /2)T,t{p4, k^) — cos^k^T /2)OT{pi, k^ 

<c{0)^(4i?) 

^ RTTT -^conn 

C2{G)C2{F) f d^p f d^k 1 kh sin(p4r/2) . 2 ^ 

= ^ J W7 J (^^T^^Rw^^^^ (^^^/'^ 

BZ BZ \ j 

sin^M/2) . „^ ^,,, sin((p + fc)4f/2) 

• sin(fe3/2) ""^'^^/' ^in((p + fe)4/2) • ^^'^^^^ 

Now using the expressions (|3.1()5j) to (|3.126j) . it has to be shown that 

hm < ^ < ^(0)^(4) _ (3_;l27) 

f^oo f^oo T 

For the parts with uj^'^^\ this is easy to show; using 



(3.125) 



one indeed gets: 



and 



hm < L(0)4'?) >eon„= hm i < S^^^J^f^ >,onn (3.128) 
f-*oo f^oo T 



hm < L^^^'uJ^^^Ij, >conn= li™ i < S^^^UJ^^^Ij, >conn= 0. (3.129) 
f^oo f-*oo T 
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All the other contributions have to be checked explicitly by evaluating the eight- 
respectively nine-dimensional integrals using a numerical integration routine 
like Vegas fill. This work, which requires lots of computer time, is still in 
progress. 



Chapter 4 

Energy Sum Rule 



In deriving the energy sum rule in the form 



V{R, 13) = Um 



(4.1) 



essentiahy three steps were necessary: 



1. Taking the derivative H1.19() of the potential, calculated on an anisotropic 
lattice, with respect to the anisotropy parameter ^, and then returning to 
the isotropic lattice ^ = 1; this led to ()1.2U() . 

2. Introducing the abbreviations r]± and using the formula derived in 
for r/+; the result was (|1.23|) . 

3. Taking the limit of large T and thereby restricting the sum over all pla- 
quettes to one fixed time slice. 



There is no problem with the second step; the formula for rj^ used there 
was proven in analytically, as well as checked numerically. But the first and 
third step are questionable; its not entirely clear if in the continuum limit, the 
potential really becomes completely independent of ^ (hence it is not clear if 
the derivative (|1.19|) really is zero), and the behaviour for large T is not clear, 
too. 

In the first section of this chapter, some basic calculations will be done, 
and the energy sum rule will be cast into an equivalent form which is easier 
to check. The second section deals then with the explicit perturbative check. 
This will lead to a better understanding of the various terms which contribute 
to the potential energy (the energy in the electric and in the magnetic fields 
and the trace anomaly). The various contributions will be classified and their 
magnitudes will be compared in the third section. Finally, the limit of large T 
will be examined. 
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4.1 Preliminaries 

It will turn out that it will be helpful to write the difference between the spatial 
and temporal plaquettes in the first term of (|1.23j) in the following way: 



-Vt + Vs 
Then the potential is given by: 



-Vt-Vs + 2rs. 



(4.2) 



V = lim — 

f^oo T 
= lim — 

= lim — 



V- < -Vt -Vs + 2Vs >gg-0 + 



290 



<S> 



qq-O 



-J- < —J + -^Js >qq~0 H 7, < >qq-0 

(3 2,gQ 

r/_ < S* >qq-0 +—^ri_ < Ss >qg-0 H 7, < S > qq^Q 



2d{F) 



d{F) 



2go 



(4.3) 



where Sg denotes the part of the action which comes from summing only over 
the spatial plaquettes. 

Because the expectation value of the action appears in two places in the 
formula above, the identity ()3.5|) 



lim — < S >qq-0,subtr- 
T — ^oo 1 



(4.4) 



becomes very useful now. When using it, one has to pay attention that the 
self-energy contributions have to be subtracted. Inserting this identity above, 
one gets: 



V 



9t 



d 



V- 



2d{F) '-dg, 



V 



2 dg( 



V + lim 



1 gl 







f -^oo T d{F) 



??_ < Ss >qq-0 (4.5) 



Until this point, everything had been exact. Now the following expansions will 
be used, where the explicit expressions for Pt and Pt were taken from j42j : 



2[dC 



1 /d_ 
2[m 



e/3+2d(F)ect(o±^/5±^c.(o 



+ 0{gi) 



(3 + 2d(F)c,(0 + 2d{F)C^ T 1/3 T ^c.(e) ± + Oig'o) 

f3 + 2d{F)—^^p±2d{F) 



2d{F) 



2d{F) dcs\ 



dct I dcs 

P + d{F){^-^ 



d{F){ 



5=1 
5=1 




f3od{F) 
^ + cd{F) 



+ 0{gl) (4.6) 
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Additionally one needs the leading order of the expansion of the lattice beta 
function: 



?o « „3 , ^/„5 

with 



/3i(ffo) = -a^ = -PogI + 0{g'o) (4.7) 



. 11^^2(0) 



487r2 

Then one arrives at: 

d ' 1 4 5 - 1 ^ 4 (9 



+ lim — 

f^ooT 

dgo 



2<Ss >gg-0 +5oC < Ss >gg-o\ + 0(5o) 

5 



5=1 f^5o 



V 1 r 
+ lim — 



where 



2<S, >gg-o +glc < Ss >gg-o\ + 0{gl), (4.8) 

has been used. The formula (|4.8j) is equivalent to the energy sum rule up to 
next-to-leading order; it will be checked now in the following section. 

All that is needed in order to do this are the expansions of the potential and 
of < Ss >gg-o- The explicit form of the potential up to next-to-leading order 
can be taken from chapter 2; for the other term, one gets the following expansion 
up to next-to- leading order, which is similar to the ones for < SW >conn (|3-2U|) 
and < LW >conn (EISU): 

< Ss >gg-0 = -gi < Si'^^J^^ >,onn +gl < 5f 5(2)a;(2) >^„„„ 

-g'o < si'^lisWfJ'^ >_„ < 5(0)42)^(2) 

-g'o < SPJ'^ >conn +g'o < Si'^S^'^J'^ >,onn 
-So" < 5(0)u;(4) -^4 < 5(1)^(3) 

-g^ < J2) >^ _^Q^g6y (4 9) 

In contrast to (|3.2flj) and (|3.52|) . here a disconnected term (the last one) appears, 
because of the definition ()1.8|) of < O >gg-o. 

For simplicity, but without loss of generality (because of the symmetry of the 
lattice), in the following the Wilson loop will be chosen to lie in the 3-4-plane. 

4.2 Pert ur bat ive check 
4.2.1 Leading order 

In leading order, (|4.8j) reduces to the following simple formula: 

y = 9o^V - 2gl lim i < Si^^uj(^^ >,onn (4.10) 

Cgo T^ooT 
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where use has been made of the fact that both the expansion for < Ss >qq-o 
and the expansion for the potential start with gg. But the latter one leads also 
to the following identity: 



(4.11) 



hence in leading order, the energy sum rule is equivalent to the requirement 

(4.12) 



lim i < 5(°)a;(2) 
r-»oo T 



i. e., in leading order, the expectation value of the spatial plaquettes has to 
vanish in the limit of large temporal extent of the Wilson loop. 

An explicit calculation for this correlator can be done by first looking at the 
insertion of si^^ into an arbitrary gluon line: 



< ^^ >conn= 6^^{27r)'6{p + 9)5^^^^^^^^<5^ 



with 



1 M<4 
/x = 4 



(4.13) 



(4.14) 



Then the result for the correlator with the Wilson loop is: 

d'^p cos(p4f) sin2(p3i?/2) ft^ 

BZ 



< 



^conn,subtr~ 



Taking the limit, one obviously obtains: 

lim i < 5(0)a;(2) > 



conn,subtr~ 



0; 



(4.15) 



(4.16) 



the expectation value of the spatial plaquettes (the energy in the magnetic 
fields) does indeed vanish in leading order, and therefore the energy sum rule 
is true in leading order. 

Looking at the relevant Feynman diagrams makes this a bit clearer. On the 
one hand, there are the possibilities where both gluon lines end on spatial lines 
of the Wilson loop: 




Here a sum over all spatial plaquettes is to be understood (which is denoted by 
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the perspective view of the plaquette), and, as usual, the end points of the gluon 
lines can be attached to all four links of the plaquette. Obviously, such diagrams 
can give no contribution proportional to T; therefore their contributions to 
conn,subtr Vanish when dividing by T and then taking the limit. On 
the other hand, the sum over all of the diagrams of the following form: 



would give a contribution proportional to T. But such diagrams can not appear, 
because the propagator in Feynman gauge connects only links which are parallel 
to each other. Only spatial plaquettes appear here, hence no gluon line can 
connect the plaquettes to temporal line of the Wilson loop. 



4.2.2 Next-to-leading order 

Using the result derived above 



lim I < SPJ^^ >conn,suUr= 0, (4.17) 
T^oo T 



one sees that only the higher terms in the expansion of < Ss >qq-o can give non- 

>qq-Q 



vanishing contributions in the limit of large T. Hence the term Qq < Ss >q 



in 1)4. 8|) is of order and does not contribute in next-to- leading order: 



T> 2 "9 T> , 9ct 
dgo 



9o^y + 2 lim ^<Ss >qq-o +0{gl). (4.18) 

5=1 ogQ r-»oo T 



The remaining graphs which contribute to < Sg >qq^o can, as usual, be 
divided into the three groups 

• vacuum polarization graphs 

• spider graphs 

• graphs with two independent gluon lines 

And again as usual, these three types of graphs will be examined separately. 



The vacuum polarization graphs 

The relevant graphs are shown in figure ETTl Using the same argument as for the 
one diagram in leading order, one sees easily that most graphs do not contribute 
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Figure 4.1: Vacuum polarization graphs contributing to the expectation value 
of the spatial plaquettes 



for T — > oo; only the second and the fourth graph survive: 

^ conn.ext 

iim ■ — = iim ^ 

f-^oo T f^oo T 

^ c(2), ,(2) ^ ^ o(0) c{2) , ,(2) ^ 

< ' >conn T < iJppl^^ ' > conn 

= nm ^ = nm 

f^oo T f-^oo T 

^ o{0) o(2) , ,(2) ^ ^ o(l), ,(2) ^ 

(because of the insertion S^s^\ respectively sf'\ only graphs can contribute in 
which at least one of the external gluon lines ends on a spatial line of the Wilson 
loop — but the sum over such graphs does not give a contribution proportional 
to T and hence vanishes in the limit). The subscript "ext" in the first two terms 
means that there, only the insertions of S^^^ into external lines are considered. 
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The remaining graphs (number two and four) will be studied later in more 
detail. 

The spider graphs 

As in section 13.4.21 there are two spider graphs which have to be taken into 
account here: 



The first graph behaves for T — > oo like < 

1)^(3) > 

conni the ordinary spider 

graph without plaquette insertions: 

lim i < 5(0)5(i)a;(3) ^-^ I < ^(1)^(3) q (4 20) 

f^oo T f-*oo T 

Treating the second graph is even more simple; it does not contribute at all, 
even for finite T: 

< 5(i)a;(3) 0^ (4 21) 

The reason for this is that si^^ contains a three-gluon vertex. Because this 
vertex only incorporates spatial plaquettes and the gluon lines are in Feynman 
gauge, all three gluon lines can only be connected to the spatial lines in the 
Wilson loop. But then the polarizations of all three gluons which meet at the 
vertex are identical — and for this configuration, the vertex vanishes. 

Therefore both spider graphs give no contributions. The only graphs which 
are left now, beside the two vacuum polarization graphs mentioned above, are 

The graphs with two independent gluon hnes 

Here one has to consider two types of contributions. First, as usual, there are the 
contributions coming from uj^^\ but, as already mentioned, here additionally a 
disconnected contribution appears. This is depicted in the following two graphs: 



78 



CHAPTER 4. ENERGY SUM RULE 



The first graph gives contributions proportional to (C2(-F))^ and to C2{F) ■ 
C2{G); the product of the two loops in the second graph gives a color factor 
(C2(-F))^. An explicit calculation yields: 

^ ijg iu ^conn 



BZ 



(27r)4 p2 lsin2(p3/2) sin2(p4/2) 



d'^p cos(p4f'/2) sin2(p3^/2) 



{p -Pt) 



and 



(27r)4 (p2)2 sin2(p3/2) 

tiZi 

+C2(F)C2(G) ■(...) (4.22) 



< >o< 5f ^(2) >eonn 

sin2(p3i?/2)sin2(p4f/2) / 1 1 



= -0(r (FW^ ( P sin-(p3i^/2)sin-(p4r/2) 1 

^ (27r)4 ^2 Uin2(|,3/2) ^sin2(p4/2 



J (^ f^2)2 ,,,2. -^^) ' 



BZ 



hence the contributions proportional to (C2(-F))^ cancel exactly, and the dis- 
connected part vanishes. 

But there are still contributions left from cv^^^ which are proportional to 
C2{F) ■ C2{G). These graphs can be classified into three categories: 

1. Graphs in which only temporal links of the Wilson loop appear; these 
links can not be connected with the spatial plaquettes using the gluon 
propagators in Feynman gauge, and therefore their contributions vanish. 

2. Graphs in which only spatial links of the Wilson loop appear; these do 
not give contributions proportional to T and therefore in the limit of large 
T, they vanish. 

3. Graphs in which temporal as well as spatial links appear - these have to 
be treated explicitly. 

In order to examine the contributions of the graphs of the third category, 
split uj^^'' again into the parts A to and then look only at the graphs with 
spatial as well as temporal links; write for these wj^/^'' to The results for 

the various contributions are: 



C2{F)C2{G) f d^p f d^k p -pI 1 sm^{p3R/2)sm^{hT/2) 



6 J (27r)4 J (27r)4 (^2)2 p sin2(p3/2) sin2(fc4/2) 

[1 — cos{k^R) — cos{piT) + cos{piT) cos{k^R)\ 

+ {C2{F)f .{...) (4.24) 



4.2. PERTURBATIVE CHECK 



79 



/ <?(0), ,(4B) 



C2{F)C2{G) r d^p f d^k V -pj l siii^{p^R/2) siii^{k4T/2) 



BZ 



BZ 



^ c(0), 

C2{F)C2{G) f d^p 



BZ BZ 
sin2(p3i?/2) / sm2(/c4r/2) 



d^/^ ^ — P3 1 
(27r)4 J (27r)4 (^2)2 ^ 



sin2(p3/2 

,2 



sm2(A:4/2) 



sin^(p3/2) 



\^ sin^(A;4/2) 



sin2(A;4r/2) ( sin^{p^R/2) 
siv?{p^/2) 



sin2(A;4/2) \ 

siv?{p^R/2) 
sin2(p3/2) 



- i? (1 - cosiksR)) 



-2 



su^{kiT/2} ■ 

COS(p4i ) — ;— 27 COS(/C3itj 



conn 

C2{F)C2{G) f d^p 



r d'^p r d'^i 

J (2^7 (2^ 



sm^{ki/2) 
d^k p -pI 1 



(27r)4 y (27r)4 (^2)2 ^2 

BZ BZ ' I 



24 

sin2(p3i?/2) sm2(A;4r/2) 



-3 



sin^(p3/2) sin^(fc4/2) 

sin^(p3i^/2) / sin^(fc4r/2) 
sin2(p3/2) \ sin2(A;4/2) 

sin2(A;4T/2) / sm2(p3i?/2) 



(8-2cos(j94T)cos(fc3i2)) 



sin2(A;4/2) 

sin^(p3i^/2) 
sin2(p3/2) 



sin2(p3/2) 



COs(]94r) 



sm2(A;4T/2) 
sin2(A:4/2) 



+ 3T 



sin^(fc4r/2) / sm^(ff3i?/2) ^ 



sin2(fc4/2) 

< 'Ujjj, >c 

C2{F)C2{G) 



12 



BZ 



d'^p 



d'^k p 



BZ 



^ sm^{pzR/2) 
sin2(p3/2) 



T {I - cos{paT)) + 



(27r)4 (p2)2 p 

s\n^{kiT/2) 



(4.25) 



(4.26) 



(4.27) 



sin2(A;4/2) 



R{l-cos{k^R)) 



.(4.28) 



And finally, because of the symmetry properties in the colour indices one gets 
as usual: 

< Sf'>U^iP >conn= 0. (4.29) 
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Adding everything, the result is: 

(0) (4) ^ _ C2{F)C2{G) f f d'^k f^-Pl 1 

i^RT>conn;C,{F)C,(G) " ^ ] J T^^^f. 



(27r)4 J (27r)4 fp2f ^2 
BZ BZ ^ ' 

sin^(p3^/2) , ^, sin^(fc4r/2) 

■ 2/ /o^ cos(p4T) . a cos(fc3i?) 
sm (^3/2) sm (A;4/2) 

(4.30) 



Divide by T and take the hmit: 



c(0) ('4") c(0) (4) 

< ^ >eonn;C2(i^)C2(G) _ -^^^ < ^flr >conn;C2(i^)C2(G) _ g 

f^oo T f^oo r 

(4.31) 

Hence nearly all contributions vanish, and only the vacuum polarization graphs 
contribute to < 5^ >gg-o in next-to-leading order! 

Check in next-to-leading order 

Making use of the fact that only the vacuum polarization graphs contribute to 
the expectation value of the spatial plaquettes, the result is now simply: 

5(0)i(5(l))2u;(2) 

>rnr 



lim i f< 5(0)5(2)^(2) >_„ - < 5(0)i(5(i))2.;(2) 



^z?^ 4 / d^p cos(p3fi) (^) 
= C'2(F)5o / (27r)3 /.j^2 nV(p,0). (4.32) 

^ ^ (f ) 

As usual, the self energy contributions have been subtracted here. The vac- 
uum polarization tensor 11^^^ (with the insertion of the sum over the spatial 
plaquettes) is given by the sum of the following two contributions (given for an 
arbitrary dimension d; see appendix El): 



ytLoop,, 



BZ *'J = 



d 1 ^2 

d k ^ — ^ k kj^kj 



6ij{p + 2k)^{p + 2k)y cos{pi/2) cos{pj/2) 

+6^,^{2f+k)i{2f+k)j cos{k^j2) cos{k^/2) 

{p - kf cos((p + k)i/2) cos((p + k)j/2) 

+ \^5^i (^2f+k)j{2f+k)^ cos{ki/2) cos{p^/2) 
+ ip^k)j{p + 2k)^, cos((p + k)i/2) cos{pj/2) 
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+ {p - k)^{2p + k)j cos{{p + k)i/2) cos{k^/2) + 



(4.33) 



T-iTadpole,s 



d 



(2d-5)(d-2) 7 

j:^! " M + 



2,d{d - 1 

( {2d-l){d-2) ^ . . 
I 4d(d-l) ^° - id ' ^ 



(5(i-2)(d-2) 



6d(d - 1) 



Ao + 



1\ 



(i-2 
Qd{d-l) 



a2 



d- 2 



6(d-l 



+ 



(4d-l)(d-2) 
12d(d - 1) 



+6] 



1 

-( 

6 

A -2d 



+ ( C2{F) - -C2{G) ) { 5^,^^^^Aop 



{d-2f 



Ao + ^ b + 



d d 
{d~2)\ . . 
didTT) ''''''' 



^1 + K 



d{d - 1) 
4 -2d 



AoP^ 



^ Ao + - ) Pt^Pu^l^K 



with 



(4.34) 
(4.35) 



and the usual abbreviation 



BZ 



d'^k 1 



(see appendix [nj. 

As already explained above, vacuum polarization graphs in which the sum 
over the spatial plaquettes is inserted into an external line can not contribute; 
therefore the ghost graphs and the graph with the insertion of the integration 
measure do not appear here, because they do not have internal gluon lines. 

Only the p, = v = d = A components contribute in the limit of large T, and 
from these, only the part with p4 = 0. Then the formulas above reduce to: 



n 



1 -2 1 



C72(G)Ao [-:P -^]+ C2{F)-Aop 



n 



44 



C2(G) 



12 

d^k 



1 



-J 



(4.36) 



1 



BZ 
i.2 



(27r)^ 



kA ( {p + kf + kl 



{2k)ik 3--jf +fc [l--ki\[2p + k 
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- 1 



(4.37) 



There is a source for a potential infrared divergence in the p-integral in H4.32() : 
one has to check that the vacuum polarization tensor goes to zero for p ^ 0. 
Inserting p = above, the results are: 

and therefore 

ng (0) = (0) + u^r'" (0) = 0. (4.38) 

Hence there is no infrared divergence, and the integral can be done without 
problems. 

Inserting the formula 1)4. 32() for the expectation value of the spatial plaque- 
ttes, (j4.18|) becomes: 



T> 2 d ^y^dct 



,t±,V^2C,iF)4 j |^j£!^nlJ(p-0) (4.39) 

=1 % BZ^ ' [P 



Then the expansion for the potential 

V = glV^^^+9tV^'^^+0{gl) (4.40) 

yields 

d^p cos{pzK) ^ 

BZ 



y(4) ^ 2V^A) dc^ 



V^^^ + 2C.iF)J^^-^uihm. (4.41) 



The explicit expression for the potential can be taken from chapter [3 Addi- 
tionally, use 



dct 



= -7C2(F) • 0.586844 + AC2{G) ■ 0.005306, (4.42) 



derived in [15], to obtain finally the following expression, which is equivalent to 
the energy sum rule in next-to-leading order: 

d^p cos{psR) (s) 



BZ 



C'2(G') 5 , nnr'T.nnA , C'2(i^), 



f 0.023220 \nR + 0.057400) + ^Mo.051644. (4.43) 

Now it is convenient to split up y\^'^<^p°^^ into two parts, corresponding to 
the group theoretical factors: 

Y^Tadpole ^ (j^(^Q-^Y\^adpole,l _^ ^^-j jjradpo«e,2 . (4.44) 
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Figure 4.2: Comparison of the calculated values to the fitted curve; see text 



then (|4.43|) splits also into two parts: 
(I) 



BZ 



(2vr) 



2 

P 



AttR 



-^0.051644 



(4.45) 



/rn f COs(p3i?) /„radpo«e,s,l/^ n\ , TjLoop,s.l ^ ^ p, 

^ ' / (27r)3 /-J^\2 V (P>0) +n44 ^' (p,0 



f 0.023220 In + 0.057400) 



(4.46) 



In the first part ()4.45p . insert the explicit expression 11^'^^°''^'^ (p, 0) 



3A0P ; the remaining integral gives in the continuum limit (see section r2. 2. 1|) : 

d?p cos{p^R) 1 



BZ 



(27r)' 



-a 
P 



AirR 



(4.47) 



Hence the first part of ()4.43p requires simply the following formula to hold: 



-Ao = 0.051644. 

3 



(4.48) 



Inserting for Aq the numerically calculated value (see appendix C), one sees 
that this formula is indeed satisfied — to very good accuracy! 

The second part ()4.46|) is a little bit more complicated. The integral on 
the left hand side can be only evaluated numerically; this has been done for 
R between 1 and 10 using again the routine Vegas from [IJ. The resulting 
numbers can be compared to the function on the right hand side. 
A fit of the calculated values to a function of the form 



alnR + b] , 



AirR 



using Mathematica, gives: 



a = 0.02268 ± 0.00078; 



b = 0.05792 ± 0.00132 



(4.49) 
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If the energy sum rule is valid, the results should be a = 0.2322 and b = 0.05740. 
The numbers from the fit agree very well with this in the range of their numerical 
errors. Hence the second part of (|4.43|) is also valid. 

For illustrative purposes, figure gives a plot of the calculated values (left 
hand side of (|4.46l) '). as well as of the curve on which they should lie if the 
energy sum rule is valid (right hand side of 1)4. 46() ). The very good agreement 
between the points and the curve is obvious. 

Summarizing: Both parts, ()4.45|) and 1)4. 46() . of the formula 1)4. 43|) . which 
is equivalent to the energy sum rule in next-to-leading order, are satisfied with 
good numerical accuracy, or in other words: 

The energy sum rule is valid up to next-to-leading order with good numerical 

accuracy. 

4.3 Contributions to the potential 

Using the validity of the energy sum rule and the results obtained for the spatial 
plaquettes, one can now analyze the various contributions to the potential: 

V = electric field energy -|- magnetic field energy 
-|-trace anomaly 

= -r/_ hm <Vt >qq-0,subtr 1™ ^ < >qg-0,subtr 

T-»oo T T-*oo T 

1 1 

+7 Jim ^ < ^T^^(x) >qg-o,subtr (4.50) 

4 T^oo T 

First, look at the contribution of the magnetic field energy. 

V— ^ 'Ps ^qq—0,subtr ~ ^ Ss ^qq—O.subtr 

— < >qq—0,subtr ~^'^^9o ^ ^qq—0,subtri 

where the constant c is again given by 

_ / dct dcs 

the functions q and Cg can be found in 42j. 

As already pointed out in the last section, in the limit of large T, the leading 
order of < Sg >qq-o,subtr vanishes, and the lowest non-vanishing contribution 
is of order g^; hence up to next-to-leading order, the second term does not 
contribute here. The first term is given by the two vacuum polarization graphs 
where the sum over the spatial plaquettes has been inserted into an internal 
line. According to the results obtained in the last section, one gets then for the 
contribution of the magnetic field energy: 




magnetic field energy 
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^ <'^s >qq-0,subtr 
T^oo T 

C2{F) (llC2{G)^^ 
AttR 
F 



^^©5o fC2(G) fo.023 In ^ + 0.058) + 0.052 C72(F)) . (4.51) 



For the electric field energy, use 

^ T^t ^ qq—0,subtr — S ^qq—0,subtr ^ 5*5 ^qq~0,subtr^ 

and 

lim ;^ < O >gq-0,SM6tr— "a"- 

T^oo T f^/j 

Inserting the explicit form of the potential and of < Sg >qq-o,subtri one then 
obtains for the contribution of the electric field energy: 

electric field energy 



-f]- ^ <^t >qq-0,subtr 

T^oo T 

C2{F) ( 2 ^ llg2(G') 4. A 1 dCs 

— 9n H n — Qn In 7.501it — 

C2(F) 



5o' + |c2(F)5o^ 
5=1 8 ^ 



, „ U +5oC^2(G)(o.070 In i? + 0.141) +0.302 ff^C2(F)) . (4.52) 
47rit ^ ^ ^ ^ 

Finally, for the trace anomaly one can use 
and again 

r 1 c _ 

lim < J >qq-0,subtr— H^T^- 
T-^oo T op 

Then one sees that the contribution of the trace anomaly to the potential is 
given by: 

trace anomaly 
1 1 

= T Jim < "^F^fii^) >qq-0,subtr 
4 T^oo T 

11C2{F)C2{G) 

967r2 4vri? 

« -0.011C2(F)C2(G)^. (4.53) 

4ttR 

Looking at the three contributions 1)4. 52() . 1)4. 51() and 1)4. 53p . one notices 
several things: 
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Figure 4.3: Contributions to the potential in NLO, proportional to (C2(-F))^, 
divided by g^: electric field energy (dashed), magnetic field energy (dot-dashed), 
potential(solid) 
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Figure 4.4: Contributions to the potential in NLO, proportional to C2{F)C2{G), 
divided by g^: electric field energy (dashed), magnetic field energy (dot-dashed), 
trace anomaly (dotted), potential (solid) 



• The contribution from the electric field energy is negative and appears 
already in leading order. Both group theoretical factors contribute in 
next-to- leading order. 

• The contribution from the magnetic field energy is positive and appears 
first in next-to-leading order; but even there, only the vacuum polarization 
graphs give a non-vanishing contribution. Again, both group theoretical 
factors contribute. 

• The contribution from the trace anomaly is negative and appears first in 
next-to- leading order; it is proportional to C2{F) ■ C2{G). 

Comparing the explicit numerical factors which appear in the three con- 
tributions, or even better, looking at the plots where the magnitudes of these 
terms are displayed (figures 1^31 and for the special case of SU{2>)^ see figure 
14. 5|) . one additionally sees that the energy in the magnetic fields is smaller than 
the one in the electric fields, so that only a part of the electric field energy is 
cancelled and the overall contribution of the energy in the fields is negative. 
The trace anomaly gives an additional negative contribution, which for large R 
is small compared to the energy in the fields. This is in strong contrast with the 
situation for a confining potential, discussed in section ll.3.3( where the trace 
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anomaly and the energy in the fields both gave exactly equal contributions to 
the potential. 

4.4 Restriction to a fixed time slice 

What remains to be shown is that one can restrict the sum over the plaquettes 
to a fixed time slice: 



1 ^ 

hm — <Vt >qq-0 

T-*oo T 

Hm i < "P, > 

f-.oo T 



jim < V[{t) >gg_0 
jim < r',{t) >gg-0, 

T->oo 



(4.54) 
(4.55) 



where V^it) respectively Pg(t) denotes the sum over all plaquettes on the time 
slice t. Using 

(4.56) 



Vs + Vt = Is 



P 



and assuming that 



>gg-0= lim < L{t) >qq-0 
T^oo T T^oo 



lim — < S >r 



(4.57) 



(see section I4.4jl , it suffices to show that the restriction to a fixed time slice 
works for the spatial plaquettes; this is equivalent to showing that 



1 



^qq~Q= lim < Ls{t) > qq-Q 
T^oo T T->oo 



Um — < S. >r 



(4.58) 



The left hand side of this equation has already been discussed in section 14. 2| 
for the right hand side, one gets an analogous expansion: 



< L, > 



s -^qq—O 



.4 




Figure 4.5: All contributions to the potential in NLO for SU(3), divided by 
gQ: electric field energy (dashed), magnetic field energy (dot-dashed), trace 
anomaly (dotted), potential (solid) 
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4.4.1 Leading order 

For calculating the correlator < Li'^^w^^^ >conn, it is convenient, as usual, first 
to look at the insertion of L^^ into an arbitrary line. The result is very similar 
to the one obtained when inserting si^^ (j^EJ; the only difference is that the 
fourth component of the momentum is not conserved: 

< Ai{p)A^{q)Lf'> >_„= 6^^ {2nf6{p + q) E ^'"Ij' ( ^t- (4-60) 

With this, one gets the following simple result for the correlator with the Wilson 
loop in leading order: 



conn 



BZ 



(27r)4y 27r^2f„-^ + «^2^ sin2(p3/2) 



■ sin(p4r/2) sin(g4T/2) cos((p + q)i{nc,i - t)), (4.61) 

where again a Wilson loop lying in the 3-4-plane has been used and nc,4 is 
the fourth coordinate of the center of the loop. For T — > oo, this correlator 
vanishes; for the special choice nc,4 = t the correlator vanishes even for every 
finite T, because then the function under the integral is odd. 

Thus one sees that in leading order, the restriction to one time slice indeed 
works: 

lim ^ <Ss >qq-o= lim < Ls{t) >gg-_o +0(5-0)- (4-62) 

f ^00 T f ^00 



4.4.2 Next-to-leading order 

The relevant graphs for the next order can be found essentially already in 
section [4.2.21 the only crucial difference is that in that section, the plaquettes 
were summed over all space-time, whereas here they are summed only over 
the fixed time slice t. Now, as usual, it is convenient to look at the three 
different types of graphs seperately. For simplicity, again only the special case 
^c,4 = t = will be treated. 



The vacuum polarization graphs 

Here one can use arguments very similar to the ones in the sections 14.2.21 and 
l,'-{.4.2l First look again at the graphs where L^^^ is inserted into an external line. 
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They give a contribution proportional to 



{2-kY J 27r sin(p3/2) 



sin(p4T/2) sin(g4T/2) 



^ \ sin(p3) sin(p4) J p2rf j^q2^ 



BZ 



Obviously the function under the integral is odd in 54 and therefore the integral 
vanishes even for every finite T. This is in accordance with the results of section 
14.2. 2[ where it was shown that in the limit of large T, the vacuum polarisation 
graphs with an insertion of S^^ into an external line do not contribute. 

Next consider an insertion of llf'^ into an internal line and the two graphs 
coming from the vertex insertions L^p and (corresponding to the first two 
graphs in figure \?>.?>\ but here only with spatial plaquettes) . These give the 
following contribution: 



d^p f dqi siiP{p3R/2) sm{pif/2) sin{qif/2) 



^Sa — 6^13 — 6413 

sin(p«/2) sin(g^/2) 



n^«(^'P4,<?4), (4.64) 



where H^" represents the vacuum polarization tensor with L^^ inserted into one 
of its internal lines respectively the two contributions with the vertex insertions 
or . In the limit T 00, this reduces to: 

^g^(^)E/||r ™'/g,f'' nti(ao.o). (4.65) 



Now first look at the gluon tadpole graph. The contribution coming from 
this graph is proportional to: 

P^BZ ^ ^ BZ ^ ' ^ + 



5'p[5pxr'-fprx)5l (4.66) 

where again the Delta-function coming from the four-gluon vertex has been 
extraced explicitly from F. Now split this four-dimensional Delta- function up 
into a three-dimensional spatial Delta-function and another for the temporal 
components. Additionally, as explained above, it suffices to look at the special 
case p4 = ^4 = 0: 



BZ BZ I [I ^ 

E ^i.fx^'iiP^ 0), {q, 0), r, s)6; {6,xf' - f,r,) 6i 
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Carrying out the four s-integrations, using the second and third Delta-function, 
gives: 

E / lPr,i^^nP + ^^t^pTiim, iQ,0),r, -r)6; ^'''^l~;'^' 6i. (4.67) 

If one compares this with the result one would have obtained if one would have 
inserted the sum over all spatial plaquettes si^\ using (|4.1H|) . and then looking 
again only at p4 = ^4 = (the only important part for T — > 00), one sees that 
one gets exactly the same result, and therefore: 

hm < 40)a;(2)5(2) lim I < 5(0)^(2)5(2) (^^gg) 

f^oo f^oo T 

holds, if only the insertion into the internal line is considered. But above it 
already had been shown that an insertion of into an external line gives for 
T ^ 00 the same result as an insertion of si^^ — therefore the result is finally 
that an insertion of L^^ into any line of the gluon tadpole graph is equivalent 
to an insertion of Sg^^ in the limit of large T. 

Exactly the same arguments can be made for the graph where L^^ is inserted 
into any line of the gluon loop graph. There one obtains: 



lim < l(0)u;(2) f5(i)f hm I < S^'^J'^ (s^'A' > 



conn - 



(4.69) 



So, summarizing, for the vacuum polarization graphs with an insertion, the 
result is: 



hm - < 40)^(2)5(2) >_„ + < lW.;(2)1 (5(1))' > 

_ < r(0)^(2)a(2) - < L(0)a,(2) q(2) > ^ 



conn 



U-< 5fa;(2)5(2) + < 5fa;(2)i [s^'))' 



lim 4 I - < S^'^LO'-^'S^'' >ronn + < S\'"uj'-''^ ( S"-'' ]^ >, 
— < 5^ )a-'( ^ S^p^p >conn ~ ^ ^ ^rneas ^conn ^ j (4.70) 

where now insertions into all gluon lines are allowed. 

What remains are the two additional vacuum polarization graphs, incorpo- 
rating L^p and , analogous to the first two graphs in figure IIOI but with 
the sum over the plaquettes restricted to a fixed time slice. The operators 

(2) 

and Ls both only contain spatial links, but in the limit of large T, only the 
4 — 4— component of the vacuum polarization tensor 11 contributes - hence the 
contributions of these graphs vanish in the limit: 

lim < l(i)5(i)u;(2) lim < if^J^^ >conn= 0. (4.71) 
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Again that is identical to the results one obtains if one had inserted Sg instead 
of Lg, divided by T and taken the limit: 

^ q(1)q{1), ,{2) ^ / cCS), ,(2) ^ 

< Js -J^ 'iJ^^ ' >conn ,. < ' >conn „ /. „o\ 

iim = hm = 0. (4.72) 

f^oo T f^oo T 



Therefore one sees that for all vacuum polarization graphs, in the limit of 
large T, the contributions coming from an insertion of Ss (and dividing by T) 
are exactly the same as the contributions coming from an insertion of Ls, so 
for these graphs, the restriction to one fixed time slice works. 



The spider graphs 

In section l4.2.2l it was seen that the two spider graphs do not contribute in the 
limit of large T. Therefore if one calculates them with Ls instead of Sg, their 
contributions should also go to zero. 

For the first spider graph, corresponding to < L^s^^to^'^'^ >conn, one can use 
the same argument as in section l4.2.2l L^^^ contains a (slightly modified) three- 
gluon vertex, which connects only gluons with a spatial polarization with each 
other, and does not conserve the fourth component of the momentum. Using 
Feynman gauge, the polarizations of the gluons at the vertex are the same as 
the one on the Wilson loop; and because there is only one spatial direction 
available on the Wilson loop, all three gluons meeting at the vertex have the 
same polarization. But for three gluons with the same polarization, the three- 
gluon vertex vanishes - hence the result is simply: 

- and that is equal to the result for < S's^-'u;*^^'' >conn- Thus obviously one has 
lim I < Si^^J"^^ >conn= lim < L«c^(3) _ (4 74) 

T^oo T f^oo 

The treatment of the second spider graph is more complicated; it has to be 
calculated explicitly. Fortunately one can use the results of section 13.4.21 for 
this; using that here only spatial plaquettes are inserted, the formula obtained 
there simplifies considerably, so that in the end one gets: 

BZ BZ -TT 

\ Xp + 2k)^cos{p3/2)Ej^^{pTk)j . ^ ^^^^^ ^ 

2i / sm((/c + q)4T/2)OR{p3 + ks, -ks) 
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{p + q)4 cos((p + k)3/2) [k -kU +{2p + k)^ cos{qi/2)k 



+2i 



p^'k"^ [k (p + k\ + (p^q) 



T.ji'^P + k)jkjk3 cos{q4,/2) 
fp'k'^ [k + qi\ i(p+k\ + [p^qj 



\ 



J 



sin(fc3i?/2) sm((p - q)iT/2) 



cos((p + k)3R/2) cos{kif/2) 



+2i 



sin(A;3/2) sm{{p - (7)4/2) 
/ 

{2f+ k)3 cos{k3/2)ki - j:j{^P + k)jkjk cos(A;3/2) 



V 



p^k'^ (k + 94) ( (p + k^ + (p^q) 



{p^2q)iCos{p3/2)(k -kl 



\ 



p^k^ (k + qjj \^ \^p + kj +{p-q 
■ sm{{k + q- p)Af /2)Or{p3 + k^, -k^ 

-2i 



{k^q)ACOs{{p + k)3/2) 


f - pI^ 


sm{k3R/2) sin((/c - g)4T/2) 


p^k? + qj 




[ 


p + k^ 






sin(A:3/2) sin((A; - (7)4/2) 



• cos((p + k)3R/2) cos{kif/2) 



(4.75) 



A careful analysis of the individual terms reveals that for large T, all of them 
go to zero, so that one indeed gets: 

hm 1 < 5W5(i)a;(3) >_„= lim < if^S^'^J^ >eo„„= 0. (4.76) 
The graphs with two independent gluon Unes 

There are again two graphs to consider here: first, the disconnected one, cor- 
responding to < L^s^uj^"^^ >conn< ^^"^^ >0) ^'^^ sccond the one containing 
a;*^^^. But using the results obtained in leading order, where it was shown 
that < L^'^^tJ^^) >conn= 0, the disconnected contribution vanishes. In the same 
way, one can show that 



< 



0. 



For the remaining graphs, it is convenient to distinguish between the following 
three classes again, as already in section 14.2.21 The argumentation is also very 
similar: 
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1. Graphs in which only temporal links of the Wilson loop appear; these 
links can not be connected with the spatial plaquettes using the gluon 
propagators in Feynman gauge, and therefore their contributions vanish. 

2. Graphs in which only spatial links of the Wilson loop appear; these will 



be discussed below. The relevant parts of a;^^^ will be denoted by 



(4) 
RR- 



3. Graphs in which temporal as well as spatial links appear - these are the 
most complicated ones and will be also discussed below. The relevant 



parts of w^^) will be denoted by 



CO 



(4) 
RT- 



As usual, is split into the parts B to F (part A has already been treated 
above), and as usual, lo^^^^ gives no contribution at all. 

Here again the results from section 13.4.21 can be used, but because only 
spatial plaquettes are inserted, the formulas simplify considerably. The graphs 
of the second category give the following result: 



< L^S' > 



C2iG)C2{F) f d^p 



BZ 



(2vr)^ 



dq4 1 



BZ 



d'^k 



(2^)4 J 2tt y ^2 



{p + ql) 



sin((/c - p)4T/2) sin((/c + q)4T/2) 



sin2(p3i?/2) sin2(A;3^/2) 



sin2(p3/2) sin2(A:3/2) 
4sin((p + k)4f/2) sin((g - fc)4f /2)0|(p3, h) 



(4.77) 



{0),„(4C) 



RR 



> 



C2{F)C2{G) f d^p f d^k f dqi 1 



p -pi 



BZ 



{27r)\l (27r)4 7 27r ^2 ^2 + ^2 j 



BZ 



P [P + Qt 



2n ^/o^sin(fc3i?/2) 



■sin^(/c4r/2)- 



2 (Si 



sin(A:3/2) 



cos(p4T/2) cos(g4r/2) 



+ 



< ^RR > 



C2{F)C2{G) 
6 



sin(/c3i?/2) 
sin(fc3/2) 

d'^p f d^k 



BZ 



(27r)4 J (27r)4 J 2tt ^2 ^2 f ^ 



Sfl(p3,-P3) + 



dq4 1 



sm{p3R/2) 
sin(p3/2) 



(4.78) 

^R{P3,h) 



pI 



BZ 



P" {p + Qa) 



■sin^(/c4r/2)- 



cos(p4T/2) cos(g4T/2) 



+2 



sin(A:3/2) 

sin(/c3i?/2) sin2(p3i?/2) sin(p3i?/2) 



S2) |p3„fc3 + 2 S,j(0, -^3) - 2 Sr(p3 - A:3, -^3) 



sin(A;3/2) sin2(p3/2) sin(p3/2) 



Si?(p3,^3) 
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sm{k3R/2) 



sin(A:3/2) 

^RR 



^r{P3, -Ps] 



(4.79) 



< L(PU^^ > 



C2{F)C2{G) f d^p f d^k f dqi 1 



r d'^p r d'^k r 

J (2^ J J2^ J 



p -Pi 



6 J (27r)4 J (27r) 

~ — TT 

• cos(p4T/2) cos(g4T/2) sin2(A;4T'/2) 



sm(fe3ii/2) 



27r k2 ^ ^2) 

sm(/c3i?/2) 



sin(A;3/2) 

sin(A;3fl:/2) sin(p3fl:/2) sin((p + k)zR/2) 



sin(A;3/2) 



sin(A;3/2) sin(p3/2) sin((p + A;)3/2) 



(4.80) 



Because of the fast oscillations of the two cosines for T ^ oo, the contributions 
of these integrals vanish in the limit of large T. On the other hand, the graphs 
of the third category give: 



< 4°)4t ^ > 



C2{G)C2{F) f d^p f d'^k f dq4 1 



r d^p r d"k r 

J (2^ J (2^7 



P 



pI 



cos{p4T/2) cos{q4T/2) 



2 J (27r) 

sin2(p3i?/2) sin2(A;4f/2) 
' sin2(p3/2) sin2(/c4/2) 

< Lf'Ji^' > 

C2(G)C2(F) / d*p f d^k f dq4 1 



27r fc2^2(^^ + ^2 



(4.81) 



r d^p r d^k r 

J '(2^ J (2^7 



p -pi 



{2tiY J 27Tk^p2('^ 



{p + ^1) 



+ 



sin2(/c4r/2) ^. 
sin2(/c4/2) 

sin2(fc4r/2) 



BZ BZ 

siT?{kiR/2) cos{p4f/2) cos(g4r/2)Sij(p3, -ps) 



sm\k4/2) 

< Lf^-t? > 
C2{G)C2{F) 



sui^{kzR/2) cos(p4r/2) cos(54r/2) 



sin^(p3i^/2) 
sin2(p3/2) 



(4.82) 



d'^p 



d^k 



BZ 



(27r)4 J (27r) 



dq4 1 



^-^ "2 



27r fc2^2(p-? + ^2) 



sin2(A;4f/2) . ^ 



+2 



sin2(A;4/2 
sin2(A;4f/2) . ^ 



sm\ksR/2) cos(p4r/2) cos (^4^/2)2^(^3, -^3) 

sin2(p3i?/2)" 



sin2(fe4/2) 



sin^(A;3i?/2) cos(p4r/2) cos(g4r/2)- 



sin2(p3/2) 



(4.83) 



< L^I'^^P > 
C2iG)C2{F) 



d'^p 



d^k 



BZ 



(27r)4 J (2n) 



dq4 1 



"-^ "2 

P -Ps 



BZ 



27r fc2^2(p-^ + ^2j 
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• '''''in7o? sin2(fc3i?/2) cos(p4T/2) cos{q^f/2)R. (4.84) 
sm (A:4/2) 

As already pointed out in section r-{.4.2[ for T ^ oo, the factor ^™ iit^ln\ gives 

a linear dependence on T, but the two cosines both give factors of T~^, so that 
in total all of these integrals go with T^^ in the limit and hence vanish. Thus 
one gets: 

hm < 4°)u;W >= hm ^ ^ = 0. (4.85) 

f^oo f^oo T 



Taking together all of the results obtained in the last subsections, one sees 
that indeed 

lim < Ls >qq-o= hm ^ (4.86) 

f^oo f^oo T 

is satisfied up to next-to-leading order — up to that order, the restriction to one 
fixed time slice works for the spatial plaquettes. Thus the check of the energy 
sum rule is now completed. 



Chapter 5 



Summary 



In this work, it was shown that in lattice perturbation theory, both sum rules, 
H1.12() and H1.23() . hold up to next-to-leading order in the coupling constant. 
Additionally, the possibility to restrict the expectation value of the action and 
the expectation value of the magnetic field energy to the sum of the plaquettes 
on a fixed time slice has been investigated. Two spin-offs of the check were 
a proof of the transversality of the gluonic vacuum polarization on the lattice 
in leading order and a proof of the gauge invariance of the expectation value 
of the Wilson loop up to next-to- leading order. This is not completely obvi- 
ous, since gauge transformations on the lattice are implemented via unitary 
transformations of the link variables. 

The scaling behaviour of the potential which is used in the derivation of the 
action sum rule was checked explicitly by using known results for the potential 
(HOlEnilSI- The crucial part of this sum rule is the identity (|1.7j) respectively 
(|^-i.5jl . The perturbative examination of this identity yielded methods and valu- 
able results which could be used in the examination of the energy sum rule. 

Additionally, it opened up a way to proof the gauge invariance of the ex- 
pectation value of the Wilson loop perturbatively up to next-to-leading order. 
By expressing this expectation value as a polynomial with respect to the gauge 
parameter, it was possible to show that its dependance on the gauge parameter 
vanishes in all orders of the gauge parameter and up to next-to-leading order of 
the coupling constant. For eliminating the dependance on the gauge parameter 
in first order, the action sum rule could be used. 

The examination of the restriction of the expectation value of the action 
to one fixed time slice (the possibility to replace this expectation value by T 
times the expectation value of the Lagrangian, the sum of all plaquettes on a 
fixed time slice) turned out to be much more difficult. It was first shown that 
this is possible in leading order. For the graphs in next-to-leading order, this 
was accomplished for the vacuum polarization graphs, the spider graphs and 
some of the graphs with two independent gluon lines; work on the rest is still 
in progress. The eight- and nine-dimensional integrals have to be evaluated 
numerically, and that requires lots of computer time. 

On the other hand, all parts of the energy sum rule were shown to be true up 
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to next-to-leading order: first it was checked without the restriction to a fixed 
time shce. It was shown that both the expectation value of the magnetic field 
energy and the contribution from the trace anomaly vanish in leading order, so 
that the only contribution to the potential in that order stems from the energy 
in the electric fields. 

In next-to-leading order it turned out that the expectation value of the 
(euclidean) energy in the magnetic fields receives its only contributions from 
the graphs in which the sum over the spatial plaquettes is inserted into internal 
lines of the gluonic vacuum polarization. By computing the contributions of 
these graphs numerically and examining the two parts of the potential with 
different group theoretical factors separately, the validity of the energy sum 
rule was confirmed with good numerical accuracy. 

For the expectation value of the magnetic energy it was possible to show 
that it can be restricted to one fixed time slice. The leading order was again 
relatively easy and could be demonstrated and explained explicitly. In the next- 
to-leading order, similar arguments as before concerning the expectation value 
of the action could be used to take care of the vacuum polarization graphs. The 
contributions from the spider graphs and from the graphs with two independent 
gluon lines were calculated explicitly and shown to vanish, as expected. 

Now the questions which were posed in the introduction can be addressed 
using the results obtained in the checks. Unfortunately the conclusions which 
were obtained here are not very helpful because they apply only if the coupling 
constant is small enough to allow perturbation theory to give sensible results. 
In contrast, in real physical systems like mesons, the coupling constant is large 
and perturbation theory breaks down. But the proven validity of the sum 
rules in the small coupling regime suggests that they are true also in the non- 
perturbative region, hence it it still possible to draw some sensible conclusions 
from the results. 

The first interesting point to notice is that there are cancellations between 
the expectation values of the energy in the electric and in the magnetic fields: 
the (euclidean) electric field energy is always negative, the magnetic field energy 
is positive. This is true for both contributions (with different group theoretical 
factors) to the potential. But here the magnitude of the magnetic field energy 
is much smaller than the magnitude of the electric field energy, in contrast to 
Monte Carlo simulations in the non-perturbative regime where they both have 
comparable sizes. Hence the problem of large cancellations between these two 
contributions which appears in these simulations does not appear in the weak 
coupling limit. 

Another result is that the trace anomaly contributes only a very small part 
to the potential for large quark-antiquark separations; this is in strong contrast 
to the case of a confining potential, where it contributes exactly one half of the 
potential energy! 

To close, I will give an outlook to the open problems, possible extensions and 
future projects. Now that the sum rules have been checked perturbatively in 
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the regime of a small coupling constant, one should use Monte Carlo simulations 
in order to perform a check in the physical regime of large coupling constant. 
There the same questions will be interesting as investigated here: the problem 
of cancellations between the contributions from the magnetic and electric field 
energies and the contribution of the trace anomaly, which is expected to have the 
same magnitude as the energy in the fields in the confining region. Additionally, 
a look at the sum rule for the glueball mass would be interesting. 

Another promising approach is using a nonperturbative model like the sto- 
chastic vacuum model of Dosch and Sominov 43 . A recent calculation of the 
quark-antiquark potential in this model and subsequent comparison with the 
predictions of the lattice sum rules gave good consistency j44j . 

On the other hand, one could use the sum rules to check the consistency of 
the results obtained in Monte Carlo simulations. One first step in this direction 
was done already in j29j . where a (corrected) version of Michael's action sum 
rule was used. 

A possible extension is the incorporation of dynamical fermions into the 
sum rules, so that one could study the effects of string breaking more closely. 
Alternatively the sum rules could be investigated for finite temperature, which 
should shed some light on the phase transition to the quark-gluon plasma. 

Hence there is much potential for future work on the lattice sum rules, and 
lots of additional interesting results can be expected. 



Appendix A 

General SU(N) formulas 



A.l Basics 

The Lie group G = SU{N) has N'^ — 1 generators. In the fundamental repre- 
sentation of the group SU{N), these generators, denoted by T^, are given by 
hermitian, tracclcss, complex N x matrices. They obey the following basic 
commutation and anticommutation relations: 

[T^,T^] = iY^fABcT^ (A.1) 

c 

{T^,T^} = -^^S^^'ida^K^+Y.dABcT^. (A.2) 

where d{F) = N is the dimension of the fundamental representation F and 
idfK^p^ denotes the <i(F)-dimensional identity matrix. The real numbers /abc 
are called structure constants. 

In the adjoint representation, the generators are denoted by t"^ and are given 
by complex [N'^ — 1) x (A^^ — 1) matrices, whose elements are: 

{t^)^^ = -ifABC. (A.3) 

A.2 Traces 

For the generators in any representation R, one always has: 

TV(T^r^) = T{R)6^^ (A.4) 

with a constant T(R) depending on the representation. For the fundamental 
representation, T{R) is simply i; for the adjoint representation, it is N. Using 
this and the commutators and anticommutators given above, one gets: 

Tr(T^T^) = ^S^'^ (A.5) 

IV(r^T^r^) = ^{dABC + ifABc) (A.6) 

r^pnArpBrpCj^D-^ ^ SaB^CD " ^ y^XfABEfcDE — dABEdcDE) 
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+ o y^XfABEdcDE + dABEfcDE) (A- 7) 

E 

Tr(t^t^) = N5ab. (A.8) 

A. 3 Sums 

For any representation i?, the quadratic Casimir operator C2{R) is defined by 

J2T^T^=:C2{R)ida^Ry (A.9) 

A 

In the special case of the fundamental representation F, one gets 

iV2 - 1 

C2{F) = (A.IO) 
and for the adjoint representation G: 

C2{G) = N. {All) 
Additionally, the following relation holds for every representation R: 

r^n^^^^. (A.12, 

where d{G) is the dimension of the adjoint representation, which is equal to the 
order of the group, i. e. the number of its generators. Hence for SU{N), the 
formula gives: 

TiR) ^ (A.I3) 

This agrees with the results for T{R) given above for the fundamental as well 
as the adjoint representation. 

Using the formulas given above and the symmetry respectively anti-sym- 
metry of the structure constants and the dABC^ the following sums can be 
evaluated: 

Y.<^AAB = (A. 14) 

A 

Y^fADEfcDE = C2{G)5ac = N5ac (A.15) 

D,E 

Y^dADEfcDE = (A.16) 
D,E 

J^dADEdcDE = i^^C2{F) - - C2{G)j Sac = ^^Sac (A.17) 

E f ABC f ABC = C2{G)d{G) = N{N^ - 1) (A.18) 
A,B,C 

E dABcfABC = (A.19) 



A,B,C 



E dABcdABC = d{G) (aC2{F) - J- - C2{G)) = 



2 (A^'-l)(iV^-4) 



A,B,C 



N 

(A.20) 



Appendix B 

Sums along the Wilson Loop 



Some sums along the Wilson loop appear so often that it is convenient to 
summarize them here. 

B.l Unrestricted sums 

The easiest sum is the one which contains only one vector potential: 

I 

R-l f-1 R-1 

= ^^("-0 + lfi) + J2 ^^("-0 + Rfi + lC')-Y^ A^ino + Rfi + fi>- Ifi) 

1=0 1=0 1=0 

f-1 

-J2Mno + fi>-li>) (B.l) 

1=0 

Insert the Fourier representation for A: 

Aa{x) = J ■^A4p)e'P-+'P-/^. (B.2) 

BZ 

This yields: 



,4 / R-l T-1 

dp ' 



R-l ^ ^ f-1 

1=0 1=0 
^4 / pip^R _ 1 pipvf _ 1 



^ J (27r)4^'^^^^' p'^2isin(p^/2) + '^°^2zsin(p./2)' 



^pipy.R+ip^T _ r _^ ^pip,.T 

'"'^ -2i sin(p^/2) '"^ -2i sin(p,/2) , 
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E 



BZ 



(27r)4 
""^ sin(p,,/2) 



sin(p^/2) I' J 



BZ 



(B.3) 



Using this, one immediately gets: 
Y,<AfA^{x) >o 



(B.4) 



-2i 



f d^P sm{p^R/2) sin(p^f/2) ip(n„+[,R/2+uf/2-x) ^i^P - 
J (27r)4 p2 



= +A5 



BZ 
h,l2 

AB f d^P sin2(p^^/2)sin2(p^f/2) 



sin(p/3/2) 



(B.5) 



(27r) 



+ 



sm.\pj2) sm.\p,/2) 



The last sum is the main ingredient of < w^^^ >o and is also crucial for calcu- 
lating the vacuum polarization graphs. 



B.2 Restricted sums 

The following sum appears in w^^-* as well as several times in u;*^^) and thus is 
needed often: 



E 

h<l2 



d^k 



(27r)4 J (27r) 



[Ap{q), A^{k)] e^(«+fc)(no+AA/2+Pf/2) 



BZ BZ 

I - 2.5^^5^^ sin((A; + q\f/2)OT{k^, q^) 

sm(kf,R/2) sm{q^R/2) 



sin(fe^/2) sin(g^/2) 
sin(fc^f /2) sin(g^ j;/2) 



sin((A; - q)uT/2) 



sm{ku/2) sin(g^/2) 
(cos(g^f /2 - k^,R/2) +ism{q^f/2 + k^R/2 



-{n^iy,R^T\ 



(B.6) 



Appendix C 

Some common integrals 



In several calculations, integrals of the form 

d'^p p^pv . . . 



I 



BZ 



(C.l) 



appear. All of them can be evaluated exactly, leaving only one constant which 
has to be determined numerically up to n = 2, and an additional constant for 
n = 3. For convenience, all results here will be given for arbitrary dimension 
d as well as for the relevant case d = A. Hence the denominator is in general 
given by: 

f = Epl- (C.2) 

/i=0 

Now first, the three most elementary integrals are: 



/ 

BZ 



BZ 



I 

BZ 



d'^p 


Pi 


{2'kY 




d'^p 


1 


{2'kY 




d^p 




(27r)'^ 


p2 



d ^^-'^ 

Ao (C.4) 
0, (C.5) 



the first following from the symmetry of the lattice, the third from the fact 
that the integrand is odd, and the second simply being a definition of the one 
remaining constant Aq. This constant can be calculated more easily by noting 
that the corresponding integral can be rewritten in the following way: 

d'^p 1 dt f ^'^^ e~*P^ 



(27r)<^p2 _/q J r2TT¥ 

BZ BZ 



BZ \ 
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dtlj ^ exp (-4t sm2(p/2) 



oo 



2 

where it has been used that the integral 



1 /■'^ 

-J^ dte-'Xit), (C.6) 



gives 27r times the Bessel function Iq (t) . The remaining one-dimensional integral 
can be evaluated using numerical integration routines much more easily than 
the original d-dimensional one. In the relevant case d = 4, the explicit result is: 

Ao ^ 0.154933. (C.7) 

From these three elementary integrals, one immediately gets: 



d> Pi _ 1 



BZ 



{27tY (p2)2 d 



-Ao. (C.8) 



The general idea to treat the more complicated integrals can be found in ap- 
pendix B of ISTI; they can be evaluated using partial integration. For example: 



= 

BZ 



d'^p d sin(p^) 

d''-p f cos(p^) sin(p^) • 4sin(p^/2) cos(p^/2) 



BZ 



(2vr)d V p2 (^2)2 

d'^p ( l-2sin(p^/2)2 _ sin2(p^/2)cos2(p^,/2) \ 

(27r)rf \ p^ (p2)2 j 

BZ 

2 (^2\2 ^'^■^> 



(2tt¥ \ p^ (p^ 

BZ ^ ' ' 

and therefore, using the elementary integrals mentioned above: 

d'^p pf, 4-2(i^ ,1 



(27r)rf(p2)2 d " d 

BZ 

Now make again use of the symmetry of the lattice: 

d'^p {ff 



^0 + -. (CIO) 



1 = 

BZ 



(2'kY ij>^? 

d f J^JL + fd'.d) fJ^^ 

J {2ttY{p^Y^^ 'J (27r)rf(p2)2 

BZ BZ 
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with fi ^ u in the second term. This gives: 

d^P pIpI 2d -a 



f dy p;pt 2d-4 



(2^)"! (]52)2 ci(d - 1) 

Using the same methods, ah integrals of the type HC.1|1 can be evaluated. 
The results up to n = 2 are given here: 

/ d'^p pIpI ^ rin^A ^^^n A A 2^-4 A 

= (5^^(-Ao + ^) + (l- V)^Ao 
A-2d^ 1,, 2(i-4 ^ 

= (-1ao + ^)V + ^Ao (C.12) 

BZ 

For n = 3, one needs an additional constant, which is defined in the following 
way: 

Ai:=(d-4) / Z"'^. ^ , (C.14) 



BZ 



Then one gets the following results: 

;Ai (C.15) 



d'^p pI 1 



BZ 

, 2 



{2t:Y {p^f 2d 

BZ 

d'^P pIpI _ 1 



(27r)'^(p2)3 2d{d-l) " d{d-l) 

BZ 

d'^p {ptf 3-2d^ 1 4 



Ao + — -Ai for fi^u (C.17) 



{27r)d (p2)3 ' d d 

BZ 

, 2 



A0 + ---A1 (C.18) 



- ' A,+ 'a, for,^. (C.19) 



BZ 



{2TT)d (p2)3 d{d-l) " d{d-l)' 



pIpIpI _ 2d -6 



-Ao - -— — -73 — -Ai (C.20) 



BZ 



(2vr)d (p2)3 l)(d-2) " d(d-l)(d-2) 

for /i 7^ 7^ A 
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For d —> 4, the additional constant Ai can be calculated exactly by extracting 
the infrared divergence: The integral 



f 

J (27r 



<fp 1 



is divergent for p ^ Q. But in that limit, can be replaced by p^: 

d'^p 1 

I 

BZ 



lim Ai = lim (ci- 4) / / ^ , ^ ^ . 



Additionally, the region of integration can be restricted to a (small) sphere 
around the origin: 

lim Ai = lim(d - 4) / d^a /"^/"'ttt^A = Hm -^R'^-\ 

where O,^ is the surface of the d-dimensional sphere. Taking the limit and 
inserting the explicit value $^4 = ^(27r)^, one finally gets: 



Appendix D 

Fourier transform of the 
potential 



In section 12.11 the following expression for the quark-antiquark potential in 
momentum space was given: 



alia) 



C2{F) 



l + 9o(«) 



/3o In 



vr 



7 H 

' 33 



31 \ A(1,0,A^) 



47r 



+ R{N) 



(D.l) 



In order to get the dependance of the potential on R, one has to perform 
a Fourier transformation. In the literature usually only the result is given; 
here I explain explicitly how the calculation can be done. The first term is 
easy — simply use the residual theorem: 



lim 



(fq 1 



,iq-R 



JqR 



lim 



oo 

hm , \ [ ^ ^ ^ e^'^^dq = Mm — , 
a^o (2'KYiR J + a^o AttR 



q dq dx 



^iqRx 



-1 



q"^ + 



AttR 



(D.2) 



In the second term, additional to the pole, there is a branch cut on the 
positive imaginary axis because of the logarithm. Thus one has to use a more 
complicated path of integration: 



d^q 



(2^)2 



= lim , , ^ 

a,l3-*0j (27r)3 Q^ + a^ 
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1 ^ 
= lim - — — / 

a,/3-^o (2Tr)^iR J 



-2 + 



lim - — — 
a,/3^0 (27r)2ii? 



^2 + a2 



7^ + q2 



In 



-|- 



(D.3) 



The first term can again be evaluated using the residual theorem, the second 
vanishes, and the third gives a contribution from the discontinuity across the 
branch cut: 



lim — — In f ^ I — lim - — „ 
/3^o AttR I /x2 I ,^^^0 (27r)2ii? 



^2ilmlnf2!^) 



i/3-e 



(2 \ 
^ 1 + lim ^ / -e"^^dq 
Ai2 y ^^0 2TrR J q ^ 

lim -J-^ In 1^1 + lim / -e'^^dx 



/-- 

J X 

PR 



(D.4) 



Now the second term gives a modified Gamma function, which can be evaluated 
exactly in the limit /3 — >^ 0: 



1 



7?^ 



1 



= hm ^— In I ^ 1 + hm r(0, (3R) 

1 //3^\ 1 
hm In -77 - hm (ln(/3i?) -|- 7) 

ln(//2i^2^ ^ 27 
4^rR ■ 



(D.5) 



Putting everything together, the potential in coordinate space is given by: 



AttR 
l + 5o(a) 



■ / 7r2i?2 3i\ ^(l,0,iV) 



33 



47r 



(D.6) 
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